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Complex  physical  quantum  systems  represent  a  diverse  variety  of  problems  in 
theoretical  physics,  ranging  from  mesoscopic  physics  to  quantum  chaos.  The  challenge 
consists  in  extracting  similar  and  universal  features  from  apparently  different  systems. 

The  energy  spectrum  is  such  a  common  characteristic,  revealing  ultimately  impor- 
tant clues  about  the  nature  of  these  systems.  Random  matrix  theory  is  an  important 
tool  to  identify  the  universal  part  and  also  to  provide  a  semiphenomenologic  theory 
for  the  system-specific  behaviour.  Standard  random  matrix  theory  is  very  success- 
ful in  explaining  many  features  of  the  energy  spectra,  but  is  restricted  to  systems 
exhibiting  chaotic  spectra. 

It  is  a  natural  approach  to  extend  this  theory  to  include  systems  that  show  a 
transition  from  the  highly  correlated  chaotic  towards  an  uncorrelated  integrable  spec- 
tral behaviour.  Of  special  interest  is  the  metal-insulator  transition  in  the  disordered 
conductor  problem  where  the  metallic  regime  is  chaotic  while  the  insulating  regime 
is  integrable. 

The  approach  used  here  defines  a  generalized  random  matrix  ensemble,  based  on 

viii 


numerical  studies  of  microscopic  models  of  the  metal-insulator  transition  It  introduces 
an  internal  parameter,  which  can  be  considered  as  the  measure  of  chaos  in  the  spectra 
or  the  measure  of  disorder  in  the  system.  The  significance  of  this  ensemble  is  that  it 
is  exactly  solvable  and,  therefore,  permits  on  a  reliable  basis  the  in-depth  study  of  the 
corresponding  physical  properties.  The  generalized  ensemble  belongs  to  a  different 
universality  class  and  is  therefore  already  interesting  in  its  own  right. 

The  general  features  of  this  new  ensemble  are  discussed,  and  a  detailed  statistical 
analysis  is  performed,  allowing  the  comparison  with  other  investigations  of  either 
numerical  or  experimental,  available  in  the  literature. 

Furthermore,  the  model  is  expanded  to  include  effects  of  an  external  parameter 
like  an  magnetic  field.  This  requires  a  dynamic  formulation  of  the  random  matrix 
approach,  which  was  pioneered  by  Dyson.  The  expanded  model  for  the  generalized 
ensemble  is  compared  to  the  standard  case  and  predictions  for  physical  systems  ex- 
hibiting transitional  statistics  are  made. 
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CHAPTER  1 
INTRODUCTION 


In  the  area  of  quantum  physics  complex  systems  have  remained  a  challenging  and 
rewarding  subject  to  investigate.  The  challenge  lies  in  the  absence  of  exact  solutions 
and  the  ineffectivity  of  perturbation  approaches  due  to  nonlinear  resonances.  How- 
ever, random  matrix  theory  has  proven  to  be  quite  useful  for  a  partial  description 
of  the  system,  mostly  in  the  chaotic  region.  Until  recently  a  random  matrix  model, 
which  can  also  capture  the  physical  aspects  of  systems  whose  characteristics  lie  be- 
tween regular  and  chaotic  behavior,  was  not  available. 

The  main  purpose  of  this  thesis  is  to  motivate,  introduce  and  investigate  a  general- 
ized random  matrix  model,  which  can  describe  certain  systems  in  the  crossover  region 
between  chaotic  and  regular  behavior.  Before  detailing  the  main  model,  it  is  useful 
to  review  the  general  methods  necessary  for  the  investigation  of  the  approach. 

1.1     Random  Matrix  Theory 

Introduced  first  by  Wigner  in  the  1950s  [1,  2,  3],  random  matrix  theory  (RMT) 
is  a  versatile  tool  for  analysing  a  variety  of  complex  systems.  Originally  envisioned 
to  be  applied  to  the  energy  eigenvalue  spectrum  of  nuclei,  the  usefulness  of  RMT 
has  expanded  into  solid  state  physics  and  into  nonhnear  and  chaotic  systems.  The 
pioneer  work  by  Wigner  has  been  refined  as  well  as  expanded  by  Dyson  and  Mehta 
[4,  5,  6,  7]  and  developed  into  a  sound  mathematical  theory. 

The  principal  idea  is  as  simple  as  it  is  elegant.  Physical  problems  are  often  solved 
by  finding  the  eigenvalues  of  a  matrix.  For  problems  that  are  'sufficiently'  complex 
it  is  reasonable  to  assume  that  the  matrix  elements  are  random  variables.  Analysing 
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an  ensemble  of  matrices  leads  then  to  a  distribution  of  the  respective  eigenvalues  and 
therefore  to  interesting  physical  properties.  Wigner  considered  the  choice  of  random 
variables  as  matrix  elements  as  a  way  to  represent  our  lack  of  knowledge  about  a 
highly  complex  system. 

However,  physical  problems  usually  exhibit  global  symmetries,  which  have  to  be 
imposed  on  the  matrix.  These  global  symmetries  restrict  the  configuration  space  of 
the  eigenvalues  leading  to  important  changes  in  the  distribution  of  the  eigenvalues. 
To  be  more  specific  let  us  discuss  a  more  explicit  example  from  here  on  (another 
important  example  is  discussed  in  the  second  chapter).  The  random  matrix  often 
represents  the  Hamiltonian  H  of  a  system,  while  its  eigenvalues  are  the  energy  levels. 
As  symmetry  we  can  choose  time  reversal  symmetry,  which  leads  then,  as  well  known 
from  quantum  mechanics,  to  a  real,  symmetric  Hamiltonian.  This  reduces  the  actual 
number  of  independent  variables  and  therefore  independent  matrix  elements  from 
A''^  to  N(N+l)/2,  whereas  the  configuration  space  of  the  eigenvalues  is  apparently 
only  N-dimensional.  The  mapping  from  the  space  of  independent  matrix  variables 
to  the  space  of  eigenvalues  is  therefore  contractive  and  associated  with  a  nontrivial 
Jacobian.  In  case  that  the  matrix  elements  or  random  variables  obey  a  standard 
gaussian  distribution,  one  can  calculate  this  Jacobian  explicitly.  This  is  instructive 
and  will  be  outHned  herein  for  the  case  of  N=2.  A  general  symmetric,  real  Hamiltonian 
takes  then  the  form 

"-(lit)  (") 

Diagonalized  the  matrix  can  be  parametrized  as 

'^=i    0      E,)  C'^) 

where  Ei,  E2  represent  the  eigenvalues  of  the  matrix  H.  Apparently  a  simple  rotation 
can  transform  the  Hamiltonian  into  its  diagonal  form  //  =  DED"^ ,  where  D  can  be 


parametrized  by  a  single  'angle'  6. 

D=(    ^-°<®)      ''<^^]  (13) 

1^  -sin(0)     cos(0)  J  ^^■'^^ 

With  these  definitions  one  can  easily  evaluate  the  functional  dependence  of  the  matrix 
elements  on  the  eigenvalues,  which  gives  the  Jacobian  and  consequently  the  invariant 
volume  dhiidhi2dh22 


dhiidhi2dh22  =  det 
resulting  in 


dh]  1  dh-\  I  dh\  i 

dEi  dEi  dEi 

dE2  dE2  dE2 

dh]  1  dhi  1  dh\  i 

ae  d@  90 


dEidE2de  (1.4) 


dhndhudh22  =  |^i  -  E2\dEidE2d<d  (1.5) 

This  argument  can  be  extended  to  N  eigenvalues  [4,  8]  leading  to  N(N-l)/2  'rotation 
angles'  Q^. 

dhu..dhii\[dh2i..dh2N--dhNi..dhiyfM  —  (1.6) 

N 
Fi(0i,..,0;V(iV+l)/2)      n      \E^-  E,\dEi..dEMde,..dQN(N+,y2  (1-7) 

The  function  Fp  (here  P  =  I)  depends  only  on  the  'rotation  angles'  0,-  and  has 
therefore  no  influence  on  the  distribution  of  the  eigenvalues  Ei.  However,  the  factor 
ni<j,,^j  \E^  —  Ej\,  which  is  often  called  the  Jacobian  factor,  depends  on  the  eigenvalues 
and  will  consequently  affect  the  probability  distribution  of  the  eigenvalues.  The  most 
obvious  feature  of  the  Jacobian  factor  is  that  it  vanishes  at  E,  =  Ej  for  I  <  i,j  <  N. 
The  probability  for  that  to  occur  would  therefore  be  of  course  zero,  and  the  likehhood 
of  two  eigenvalues  being  very  close  to  each  other  would  obviously  be  quite  small.  This 
fundamental  result  is  called  level  repulsion  and  has  a  profound  impact  on  the  behavior 
of  the  eigenvalues. 

Note  that  the  derivation  has  only  assumed  a  matrix  with  random  elements  and 
a  symmetry  condition  imposed  on  the  matrix.  The  geometrical  restriction  from  the 
space  of  independent  parameters  to  the  space  of  eigenvalues  originates  already  from 
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the  appearance  of  the  Jacobian  factor  and  therefore  the  level  repulsion.    This  very 
general  derivation  shows  the  robustness  and  generality  of  the  level  repulsion. 
The  Jacobian  factor,  however,  will  be  altered  by  the  type  of  symmetry  imposed  on 
the  matrix.  In  the  example  above  a  real  symmetric  matrix  has  been  chosen,  while  in 
physics  two  more  symmetries  are  relevant. 

In  general  a  matrix  can  also  exhibit  complex  numbers  as  matrix  elements  with  the 
hermiticity  condition  H^j  =  z/f^,  being  imposed  on  the  matrix  enforcing  real  eigenval- 
ues. Then  a  unitary  transformation  U  can  be  used  to  diagonalize  the  Hamiltonian. 
A  hermitian  matrix  has  A^  +  2(|)(A^^  -  N)  =  N"^  independent  parameters,  which  are 
mapped  onto  the  N-dimensional  eigenvalue  space.  Consequently,  the  same  rationale 
as  above  can  be  applied  to  obtain  the  invariant  volume  element 

dhii..dhiMdh2\..dh2N--dhNi--dhNN  =  (1-8) 

N 

F2(ei,..,0Ar(/v-i))    n    \E^-  E,\^dE,..dENdQi..dQN^r^_^^  (1.9) 

The  Jacobian  factor  shows  in  case  of  unitary  symmetry  Y[\Ei  —  Ejp  with  strong 
similarities  to  the  case  of  orthorgonal  matrices.  Only  the  power  of  the  difference 
coordinate  is  changed.  It  should  not  surprise  us  that  the  level  repulsion  is  stronger 
for  hermitian  matrices  because  of  the  stronger  restriction  of  the  mapping  from  A'"^  to 
A^  dimensions  compared  to  A'^(A^  +  l)/2  to  N . 

This  result  also  appears  in  the  last  relevant  symmetry,  where  the  matrix  cannot  be 
represented  by  complex  or  real  numbers  but  by  quaternions.  Physical  examples  will 
be  outhned  later  in  this  chapter,  however,  it  should  be  mentioned  here  that  systems 
with  spin  ~  and  no  rotational  symmetry  will  be  represented  by  this  type  of  matrix. 
Even  though  the  derivation  for  the  invariant  volume  element  is  more  involved,  it  is 
nevertheless  straightforward  and  can  be  found  in  the  literature  [9,  4,  5]  with  the  result 


dhii..dhiP4dh2i..dh2N--dhf^i..dhj^j^  =  (1-10) 


N 


F4(0i,..,e^(27v-i)    n    \E,-  E,\'dE,..dENdQ,..dQj^(2N-i)  (l.li; 


>N(2N-l) 


5 
Again  the  Jacobian  factor  is  only  changed  in  the  power  of  the  difference  variable 
which  can  be  represented  by  the  symmetry  parameter  /3.  By  use  of  group  theoretical 
arguments  it  can  be  shown  that  /3  can  only  take  the  values  1,2,4  depending  on  the 
symmetry  case  [4,  5].  Usually  the  three  cases  are  called  after  the  type  of  diagonaHzing 
transformation  necessary  and  the  distribution  of  the  random  variables. 

For  real  symmetric  matrices  the  ensemble  of  matrices  is  therefore  named  Gaussian 
Orthogonal  Ensemble  (GOE),  for  hermitian  matrices  the  ensemble  is  called  Gaus- 
sian Unitary  Ensemble  (GUE)  and  for  quaternion  matrices  the  ensemble  is  known  as 
Gaussian  Symplectic  Ensemble  (GSE). 

We  would  like  to  discuss  the  physical  implications  of  these  three  basic  ensembles. 
The  one  which  is  usually  the  easiest  to  treat  i=!  the  gaussian  unitary  ensemble  in  which 
no  symmetry  except  for  the  hermiticity  of  the  Hamiltonian  matrix  is  present.  This 
situation  is  suitable  for  systems  which  are  measured  with  a  large  applied  magnetic 
field  in  order  to  break  the  time  reversal  symmetry.  In  case  the  time  reversal  symmetry 
is  present  in  the  system  the  analysis  has  to  be  further  subdivided  into  the  cases  of 
integer  and  half-integer  total  spin.  It  is  relatively  easy  to  show  that  in  the  case  of 
integer  spin  a  basis  can  be  chosen  in  which  the  matrix  elements  are  real  numbers  [4]. 
The  random  matrix  ensemble  which  is  describing  a  real  symmetric  matrix  is  the  GOE. 
This  applies  usually  to  systems  which  possess  an  integer  spin  and  which  are  measured 
without  a  magnetic  field  present.  The  last  case  is  intact  time  reversal  symmetry  and 
half-integer  spin.  The  finding  of  an  appropriate  basis  is  slightly  more  involved  and 
leads  to  2x2  matrices  as  basic  'numbers',  usually  called  quaternions.  The  resulting 
matrix  is  then  a  symplectic  unitary  matrix  [4,  5],  leading  to  the  GSE. 

To  discuss  the  aspects  of  random  matrix  theory  further  we  have  to  introduce  a 
distribution  function  of  the  eigenvalues,  which  is  laid  out  next. 
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1.1.1     The  Joint  Probability  Distribution  Function 

The  central  quantity  to  completely  describe  a  random  matrix  ensemble  is  the  Joint 
Probability  Distribution  function  (JPD).  An  elegant  derivation  of  the  JPD  for  ran- 
dom matrices  can  be  realized  by  defining  an  Information  Entropy  analogously  to  the 
usual  entropy  in  statistical  mechanics.  Maximization  of  the  entropy  with  some  given 
constraints  results  then  in  the  JPD.  This  argumentation  was  first  used  by  Balian  [10] 
and  Porter  [9].  The  information  entropy  is 

I{P{H))  =  j  dVP[H)\nP{H)  (1.12) 

Here  P(H)  represents  the  probability  distribution  for  the  matrix  H,  while  dV  is  the 
invariant  volume  element.  The  information  entropy  is  also  subject  to  the  condition 
that  the  probability  is  normalized  to  one. 


jdVP{H)  =  l  (1.13) 


One  can  easily  see  that  this  single  constraint  is  not  sufficient  to  define  a  suitable  prob- 
ability distribution.  The  level  repulsion  from  the  Jacobian  would  lead  to  a  nonzero 
probability  that  eigenvalues  are  at  infinity  and,  therefore,  the  distribution  cannot  be 
normalized.  Hence  another  constraint  is  necessary  to  confine  the  eigenvalues  around 
the  origin.  The  usual  choice  of  constraint  is  to  postulate  the  following  trace  condition 

I  dVTr[HH+]P{H)  =  Const.  (1.14) 

This  condition  prevents  very  large  eigenvalues  in  the  random  matrix  restricted  by  the 
finite  constant  given.  The  choice  of  the  trace  of  the  matrix  is  a  natural  one,  as  the 
trace  is  invariant  under  all  types  of  transformations  in  consideration,  i.e.  orthorgonal, 
unitary  and  symplectic.  It  can  be  even  shown  that  each  invariant  of  a  NxN  matrix 
under  a  nonsingular  similarity  transformation  can  be  expressed  in  terms  of  the  first 
N  powers  of  the  trace  [4,  11].  This  assumption,  however,  is  not  the  most  general  one 
in  order  to  construct  a  random  matrix  ensemble.  That  point  will  be  important  in  the 
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discussion  of  a  generalized  model  in  the  second  chapter. 

The  information  entropy  and  the  two  constraints  lead  then  to  a  maximization  problem 
with  two  Lagrange  multiphers  A,/i 

I[P{H),  X,i2)  =  I  dV{P{H)  HP{H))  +  XP{H)  +  ^Tr{HH+)P{H)}        (1.15) 

The  standard  procedure  is  to  set  the  first  variation  in  P(H)  as  a  function  of  A  and  /i 
equal  to  zero,  resulting  in 

5I{P{H),X,  fi)  =  I  dVSP{H){ln{P{H))  +  1  +  X  +  fj,Tr{HH+}  =  0  (1.16) 

The  choice  of  5P{H)  is  arbitrary,  and  therefore  the  integrand  has  to  be  zero  finally 
giving 

P{H)  =  exp{-{l  +  X  +  ixTr{HH^))  (1.17) 

By  expressing  now  the  JPD  in  the  diagonal  representation  with  the  Jacobian  factor, 
noting  that  Tr{HH^)  =  Y.i  ^h  ^^^  rescaling  the  energy  variable  E^  =  y/JIei,  we 
arrive  at 

P{E.)  =  Cm,01[  I^^-  -  E,f^xp{-J2Ef)  (1.18) 

Here  C^^p  represents  a  normalization  constant  in  order  to  satisfy  Equation  1.13. 
Equation  1.18  is  a  fundamental  result  derived  only  with  the  assumptions  of  global 
matrix  symmetry  and  statistical  independence  of  linear  independent  components. 
It  is  common  to  rewrite  Equation  1.18  in  the  more  suggestive  form 

P{E.)  ^  C\,pexp{-(3J2H\E^  -  E,\)  -  E^')  (1-19) 

The  surprising  fact  is  that  the  exponent  exhibits  the  same  behavior  as  the  2D-electron 
gas  in  a  harmonic  oscillator  potential.  This  so-called  Coulomb  gas  analogy  is  not  only 
important  for  depicting  a  particular  system,  but  also  for  applying  known  approxima- 
tion techniques  from  the  2D-electron  gas.  Furthermore,  it  allows  to  speak  of  the 
eigenvalues  as  particles  (electrons  or  charges)  which  is  a  useful  and  common  analogy. 
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From  equation  1.18  we  can  also  introduce  the  terminology  of  the  two  body  inter- 
action or  logarithmic  repulsion  represented  by  the  first  term  in  the  exponent  and  the 
one  body  potential  or  confining  potential  which  is  the  second  term  in  the  exponent. 
So  far,  the  discussion  has  been  limited  to  the  distribution  of  eigenvalues,  which  is 
sufficient  to  find  the  probability  distribution  and  its  moments.  However,  the  question 
can  be  extended  to  the  eigenvectors  of  the  system  which  contain  information  not  avail- 
able in  the  eigenvalue  distribution  alone.  The  eigenvector  distribution  is  in  particular 
important  in  cases  where  eigenlevel  transition  are  triggered  by  external  influence.  An 
example  is  the  description  of  nulear  reaction  width.  Also  the  description  of  localized 
or  extended  states  in  the  theory  of  localization  needs  to  know  complete  states.  Porter 
[9]  was  the  first  to  realize  that  it  is  possible  to  determine  the  distribution  function  for 
the  nth  component  (1  <§:  n  -C  N).  of  an  eigenvector.  There  result  was 

P{^)  =  \h ^exp{~x/2)dx  (1.20) 

V  Zn  -y/X 

where  a]  =  Nx  with  a,-  represents  the  nth  component  of  the  eigenvector.  Naturally  it 
is  called  Porter-Thomas  distribution.  One  important  quantity  which  uses  the  eigen- 
vector distributions  is  the  inverse  participation  ratio.  It  measures  how  extended  the 
state  is  and  presents  a  helpful  tool  in  the  theory  of  disordered  conductors. 

1.2     Correlation  Functions  and  Spectral  Statistics 
1.2.1     Correlation  Functions 

The  Joint  Probability  Distribution  function  contains  all  information  about  an  en- 
semble. However,  usually  it  is  not  necessary  to  know  the  entire  distribution  (which 
would  be  also  quite  cumbersome),  since  it  is  sufficient  to  concentrate  on  the  lower 
moments  of  the  distribution  (or  correlation  functions),  which  contain  the  major  phys- 
ical aspects.  A  natural  choice  are  the  n-level  correlation  functions.  These  are  defined 


&s 


N\ 


iw        f     f 
/?Ar(xi,..,x,/v)  =  ,  ,,  y  ■■J  dxr,+i..dxMP{xi,..,XN)  (1.2i; 


9 

The  function  Rn  describes  the  probability  density  that  one  eigenvalue  is  at  xj,  one  at 
X2,  Up  to  xat.  The  1-level  correlation  function  (also  equivalent  to  the  first  moment) 
is  the  level  density. 

/oo  /-oo 

••/       P{x,X2,..,XN)dx2-  -dXM  (1.22) 

-oo         J  —oo 


Even  though  the  n-level  correlation  functions  have  a  direct  probability  meaning,  for 
computational  purposes  it  is  useful  to  define  separat  correlation  functions  which  are 
only  nonzero  when  the  n  different  levels  are  correlated,  thereby  eliminating  level  self- 
correlations.  Such  a  set  of  functions  is  called  Cluster  functions  and  is  usually  defined 
as  the  logarithmic  derivative  of  the  n-level  correlation  function.  We  have  first  to 
define  the  generating  function  for  the  /?„ 

oo       i       .  .  i 

Roc{t)  =  J2  -\j  ■■  J  M^u  -,  Xn)  n  2/;^^.]  (1-23) 

The  generating  function  for  the  Cluster  functions  is  similarly  defined. 

oo       I  '„  ^  i 

Tooit)  =  E  ^[/  -  /  U^u  -  ^n)  n  yjdx,]  (1.24) 

1=1    •  J-1 

Now  the  Cluster  functions  are  derived  from  the  following  relation 

i?oo(0  =  exp(r^(t))  (1.25) 

resulting  in 

X        1  51n/?oo(l) 

The  1-level  and  2-level  correlation  function  in  terms  of  Cluster  functions  are 

Ri{x^)  =  T^{x,)  (1.27) 

Rl{x^,X2)   =  Ti{x,)T,{x2)  +  T2{XuX2)  (1.28) 
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1.2.2     Statistical  Analysis 

The  goal  of  statistical  analysis  is  to  provide  a  tool  for  comparing  random  matrix 
theories,  analytical  models,  numerical  calculations  and  experimental  data  on  the  sta- 
tistical properties  of  a  set  of  eigenvalues.  Even  though  statistical  analysis  can  usually 
not  detect  any  specific  physical  aspects,  it  is  quite  useful  to  understand  broad  physi- 
cal aspects  of  certain  systems  and  to  reveal  universalities  and  quahtative  features  in 
complex  systems. 

Comparing  the  statistical  properties  of  different  systems  requires  to  'unfold'  the 
spectrum.  This  procedure  has  to  be  followed  by  both  the  random  matrix  theory 
and  the  experimentally  and  numerically  observed  spectra  to  obtain  comparable  data. 
The  goal  is  to  fix  the  first  moment  properties  of  the  spectrum,  thereby  attributing 
the  se-cond  moment  properties  a  universally  comparable  meaning.  In  the  case  of  ex- 
perimental or  numerical  data  a  level  sequence  {£^„}i<„<Ar  is  given.  Then  one  analyses 
the  staircase  function  N{E)  of  that  set,  where  the  staircase  function  is  defined  as  the 
number  of  levels  with  an  energy  smaller  than  E.  This  function  can  be  extended  to  a 
smooth  approximating  curve  N{E)  (for  example  by  spline  interpolation).  Mapping 
the  spectrum  {En}  onto  a  new  set  of  values  {E'^}  (the  unfolded  set  of  eigenvalues) 
by  means  of 

E'n  =  A[N{En)  -  B]  (1.29) 

(A,B  are  two  arbitrary  constants)  results  then  in  a  constant  level  density,  which 
becomes  quite  obvious  from  the  straight  hne  increase  in  N{E').  The  most  common 
choice  is  to  set  the  average  spacing  to  1.  Now  the  fluctuations  given  by  some  second 
moment  expressions  provide  the  leading  behavior. 

The  theoretical  approach  will  usually  also  exhibit  a  nonuniform  density.  By  means  of 
the  mapping 

a{E)dE  =  a{E')dE'  =  dE'  (1.30) 
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which  requires  the  invariant  volume  to  remain  constant,  one  arrives  by  definition  at 

a{E')^l. 

The  nearest  neighbour  spacing  distribution  is  defined  by  the  probabiHty  to  find  the 

spacing  AE  between  any  two  eigenvalues.   As  such  it  is  a  close  range  characteristic 

of  the  eigenvalue  distribution.  For  a  random  matrix  distribution  the  Wigner  surmise 

for  real  symmetric  matrices  holds  true. 


nX  ,        TTX^ 


^(^)  =  ^exp(-^  (1.31) 

D  represents  here  the  average  spacing.  Wigners  argumentation  [2]  was  based  on  the 
solvable  2x2  matrix  system  and  then  generalized  for  larger  matrices.  It  should  be 
emphasized  that,  despite  the  fact  that  1.31  is  not  exact,  it  matches  the  exact  function 
surprisingly  well  [4].  For  all  three  ensemble  types  the  spacing  distribution  reads 

Pq[x)  =  Apx^exp{-Bpx^)  (1.32) 

with  Ap  and  Bp  being  constants.  The  slope  at  the  zero  argument  is  exact  even  for 
larger  numbers  of  eigenvalues  and,  therefore,  represents  a  possible  way  to  distinguish 
the  ensembles  as  well  as  distributions  which  differ  from  the  three  standard  types,  i.e. 
orthorgonal,  unitary  and  symplectic. 

In  considerable  contrast  to  the  Wigner  surmise  is  the  distribution  of  completely  ran- 
dom levels,  the  Poisson  statistics 

1  X 

P(x)  =  -exp(--)  (1.33) 

The  two  distributions  are  plotted  in  figure  1.2.2  demonstrating  the  numerical  ease  to 
differentiate  statistical  signatures.  This  as  well  as  the  simple  evaluation  for  numerical 
eigenvalues  makes  the  nearest  neighbour  spacing  distribution  a  popular  tool  for  the 
statistical  analysis. 

The  number  variance  is  another  important  quantity  for  statistical  analysis.  It  is 
defined  as  [4] 

<n^>=f     I    {R^{x)5{x-x')  +  R2{x,x'))dxdx'  (1.34) 
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Figure  1.1:  The  nearest  neighbour  spacing  distribution  for  Poisson  and  Wigner-Dyson 
ensembles  in  terms  of  the  variable  s  =  x/D,  where  x  is  the  coordinate  and  D  represents 
the  mean  level  spacing.  The  dotted  curve  is  /?  =  2  and  the  dashed  one  for  /?  =  1. 


13 

where  L  is  the  range  (support)  on  which  eigenvalues  are  present.  In  terms  of  a  fixed 
mean  level  spacing  D,  after  the  spectrum  has  been  unfolded,  the  number  variance 
reads 

/■-A/'/2     /■- ./V/2 

<n^>-<n>^=/  /  (Six- x) -T2(x,x'))dxdx'  (1.35) 

Jn/2     Jn/2  ^        ' 

It  measures  the  statistical  deviation  of  the  eigenvalue  system  and  involves  as  such 
only  information  of  the  first  and  second  moment  of  the  distribution.  In  contrast,  the 
nearest  neighbour  spacing  distribution  contains  information  of  all  moments. 

Another  integral  long-range  measure  of  eigenvalue  correlations  is  called  As-statistics. 
The  As-statistics  was  introduced  by  Dyson  and  Mehta  [7]  as  their  third  A-statistics  to 
measure  the  long-range  statistical  properties  of  eigenvalues  for  a  finite  set  of  eigenval- 
ues. Only  the  As-statistics  proved  to  be  a  reasonable  measure  for  eigenvalue  rigidity. 
It  measures  the  deviation  from  a  totally  rigid  spectrum,  that  is,  the  deviation  from 
an  equally  spaced  eigenspectrum.  More  precisely,  it  measures  the  variations  of  the 
staircase  function  around  a  best  fit  line. 

We  will  discuss  in  the  following  the  derivation  of  a  general  expression  for  the  A3- 
statistics,  since  this  measure  is  not  so  widely  known.  Then  we  will  calculate  the 
Aa-statistics  for  the  limiting  cases  of  Poisson  and  Wigner-Dyson  statistics. 

Assume  we  have  a  finite  sequence  of  energy  levels  (they  can  come  from  experimental 
or  numerical  simulations).  We  can  scale  the  eigenvalues  onto  an  energy  level  sequence 
with  constant  density  p  —  Const.  As  discussed  above  with  the  appropriate  rescaling, 
one  can  always  achieve  a  unit  density.  The  set  of  levels  ranges  without  restriction 
from  0  to  Etot-  The  staircase  function  is  defined  as 

N{E)  ^  l'^  p{E)dE'  (1.36) 

The  variations  around  a  straight  line  are  given  by 


A3     -     min{N{E)-aE-b)^  (1.37) 

=    min-  r  dE\N{E)-aE-bf  (1.38) 
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The  average  denoted  by  the  solid  hne  takes  a  sufficient  range  of  energy  levels  into 

account.  The  parameters  a,  6  have  to  be  chosen  so  that  the  extremal  conditions  are 


met.  The  conditions  ^  =  0  =  ^  lead  to  a  =  ^EN{E)  and  b  =  N{E).  Then 


A3  =  A^2  _  yv    -  — £;A^(£;)  (1.39) 

To  obtain  a  statistically  relevant  measure  we  have  to  take  a  statistical  average  over 
an  ensemble  of  eigenvalue  configurations.  Denoting  <>  as  the  ensemble  average  we 
have 

A3=<A3>  (1.40) 

In  terms  of  the  density-density  correlation  function  this  can  be  expressed  as 

A3  =  -  TdEf    (IE  I    dE'  <piE)p{E')>  (1.41) 

n  Jo         Jo  Jo 

\  I       rn       ^     rn     rE  pE 


-—  r  dE  r  dE  f    dE  f    dE'  <p{E)p{E')  >  (1.42) 

n"'  Jo  Jo         Jo  Jo 

r  dE  r  d'E  [    dE  f    dE'  <p{E)p{E')>  (1.43) 

JO  Jo  Jo  Jo 


A'n?  Jo  Jo 

The  density-density  correlation  function  can  be  calculated  from  the  first  two  correla- 
tion functions.  Using 

<  p{E)p{E')  >=  5{E  -  E')p{E)  +  p{E)p{E')  +  Y^iE,  E')  (1.44) 

and  observing  that  Y2  —  T2/ D'^  one  can  analytically  calculate  the  As-statistics.   For 
the  Poisson  case  (1^2  =  0)  we  easily  get 

Id  ' 

For  the  GUE  {Y2  =  ^^j^^Ti)  this  can  also  be  integrated  in  the  large  N  limit  (see  [5]) 

A3,Wigner- Dyson   =   — [ln(27rn)  "t"  CoUSt.]  (1-46) 

We  can  already  see  the  drastic  difference  between  the  linear  increase  in  the  Poisson 
case  and  the  slow  logarithmic  rise  for  the  Wigner-Dyson  model,  which  so  apparently 
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indicates  the  difference  between  a  chaotic  and  a  regular  spectrum.   The  large  N  be- 
havior of  the  As-statistics  is  quite  similar  to  the  number  variance,  only  at  small  N  do 
these  two  statistics  differ  qualitatively. 

1.3     Application  to  Physical  Systems 

The  wide  physical  application  of  random  matrix  theory  makes  it  almost  impossible 
to  provide  a  comprehensive  coverage.  Therefore,  we  will  restrict  ourselves  to  several 
important  cases.  More  examples,  pertaining  to  the  transitional  model  at  the  core  of 
this  work,  will  be  discussed  later. 

1.3.1     Nuclear  Physics 

As  already  mentioned,  the  first  application  of  random  matrices  in  physical  systems 
was  in  nuclear  physics.  It  was  of  particular  importance  to  understand  the  excitation 
spectrum  of  nuclei,  typically  determined  by  slow  neutron  resonance  measurements 
[12].  The  lower  lying  energy  levels  are  quite  well  understood  by  an  independent  par- 
ticle model,  where  the  nucleons  move  independently  in  an  average  potential  [13]  This 
approximation  breaks  down  for  higher  lying  excitations,  and  the  subject  becomes  a 
full  interacting  N-body  problem,  where  N  is  the  sum  of  protons  and  neutrons.  The 
high  level  of  interaction  makes  the  exact  prediction  of  energy  levels  almost  impossible. 
This  was  the  reason  for  Wigner's  proposal  in  1955  [2]  to  accept  the  lack  of  information 
as  an  advantage  and  to  assume  the  energy  level  distribution  (at  least  locally)  to  be 
equal  to  the  distribution  of  eigenvalues  of  a  random  matrix. 

It  is  obviously  quite  difficult  to  prove  the  exactness  of  a  random  matrix  approach, 
because  specific  models  are  usually  neither  solvable  nor  as  general  in  nature  as  the 
random  matrix  approach.  Therefore  one  has  to  resolve  to  numerical  comparison  and 
statistical  analysis  in  order  to  show  the  accuracy  of  Wigner's  assumptions.  In  the 
case  of  nucleonic  behavior  the  two  level  properties  (as  for  example  number  variance, 
Aa-statistics)  were  in  good  agreement  with  the  predictions  from  RMT  (a  good  com- 
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prehensive  reference  is  the  book  by  Porter  [9]).  In  contrast,  the  one  level  quantities 
like  the  constant  level  density  showed  less  satisfactory  agreement  [5].  This  represents 
one  example  where  the  statistical  behavior  demonstrates  good  agreement,  whereas 
direct  correlation  functions  still  contain  specific  physical  properties  complicating  a 
comparison  with  RMT. 

1.3.2     Quantum  Chaos 

The  first  indication  that  random  matrix  theory  is  relevant  for  describing  the  spec- 
tral properties  of  chaotic  systems  resulted  from  studying  the  quantum  stadium  bil- 
liard [14].  A  stadium  billiard  system  is  either  integrable  or  chaotic.  The  authors 
showed  that  in  the  chaotic  case  (noncircular  stadium)  the  level  spacing  was  following 
a  Wigner-Dyson  distribution  opposed  to  a  Poisson  distribution  in  the  integrable  case 
(circular  stadium).  Today  we  know  that  the  dynamics  of  a  classical  system  and  the 
spectrum  of  its  quantum  analogue  are  tightly  related.  If  the  spectrum  exhibits  level 
repulsion,  the  classical  system  is  chaotic.  Due  to  this  phenomenon  random  matrix 
theory  plays  an  integral  role  in  the  understanding  of  chaos.  The  better  understand- 
ing of  the  transition  between  these  two  limits  has  remained  an  outstanding  topic. 
In  classical  systems  constants  of  motion  are  destroyed  by  nonlinear  resonances.  The 
quantum  counterpart  is  the  destruction  of  good  quantum  numbers.  As  long  as  con- 
stants of  motion  (good  quantum  numbers)  are  present,  degeneracies  in  the  spectrum 
exist,  which  imply  at  least  some  Poissonian  behavior.  Classically  chaotic  systems  can 
also  exhibit  many  local  features  in  their  phase  space  (scars,  KAM  surfaces)  which 
make  the  transition  to  a  regular  spectrum  to  some  degree  system  dependent.  Despite 
these  possible  objections,  the  search  for  a  transitional  theory  represents  an  important 
topic  in  quantum  chaos.  An  early  phenomenological  distribution  by  Brody  [15]  has 
often  been  used  to  fit  numerical  data  of  transitional  systems.  Another  approach  was 
based  on  the  division  of  the  phase  space  into  chaotic  and  integrable  regions  [16],  where 
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the  distribution  function  was  superposed  from  the  (known)  distributions  in  chaotic 
and  integrable  portions  of  the  phase  space.  However,  a  large  class  of  transitional 
systems  seems  to  exist,  which  do  not  follow  the  characteristics  of  this  transitional 
model.  Further  relevant  and  more  recent  models  are  described  later  in  the  discussion 
including  the  model  presented  herein. 

1.3.3     Mesoscopic  Physics 

The  rising  interest  in  mesoscopic  physics  is  another  field,  where  random  matrix  the- 
ory can  be  applied  to.  A  mesoscopic  system  is  usually  considered  to  range  between 
macroscopic  physics,  where  systems  can  be  treated  as  self  averaging,  and  the  micro- 
scopic level,  where  fluctuations  can  be  very  large.  In  other  words,  the  system  is  still 
large  enough  (>  10~^m)  to  exhibit  a  bulk  statistical  behavior  but  small  enough  for 
fluctuations  effects  to  occur  (<  10~^m).  The  usual  quantity  for  defining  a  mesoscopic 
system  size  is  the  electron  phase  coherence  length  /g. 

One  particularly  interesting  area  in  mesoscopic  physics  are  disordered  conductors, 
which  exhibit  a  transition  from  metallic  conductors  to  insulators  depending  on  the 
disorder.  The  basic  assumption  here  is  that  the  quantum  coherence  length  is  larger 
than  the  mean  free  path  of  electrons,  therefore  the  physical  transport  properties 
should  be  insensitive  to  microscopic  details.  The  main  quantity  to  be  determined  in 
the  disordered  conductor  problem  represents  the  conductance  itself,  which  is  defined 
by  a  whole  conductance  distribution  P(g).  This  can  be  seen  from  the  fact  that  the 
conductance  has  to  be  a  random  variable,  which  fluctuates  with  each  realization  of 
the  disorder  potential.  The  usual  transport  theory  is  based  on  the  Landauer  formula, 
which  describes  the  electrons  to  be  in  momentum  channels,  flowing  from  reservoir 
with  chemical  potential  /ii  to  a  reservoir  with  chemical  potential  ^2-  The  current 
depends  then  in  the  following  way  on  the  voltage 

t^[jj:T.,]V  (1.47) 


.2 
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where  T  represents  the  transmission  matrix  and  y  the  conductance  quantum.  Then 
we  have  for  the  conductance 

oo 

5  -  E  T^J  =  TrTT^  (1.48) 

This  simplified  approach  was  pioneered  by  Landauer  [17,  18]  and  refined  later  in 
[19,  20].  The  formula  can  be  used  to  establish  a  scattering  theory  approach  which 
relates  ingoing  flux  into  an  interaction  area  to  the  outgoing  flux.    The  scattering 

I  matrix  is  usually  replaced  by  the  transfer  matrix  M  which  relates  fluxes  on  the  right 

•1 
1  to  the  fluxes  on  the  left. 

!  *^(o)  =  (o')  ('-«' 

The  current  conservation  law  requires  then  the  relation 

I  |/|2_|0|2  =  10'|-|/'r  (1.50) 

An  appropriate  choice  of  parametrization  and  the  relation  which  connects  the  transfer 
and  scattering  matrix  leads  then  by  use  of  the  Landauer  formula  to  the  conductance 
[21,  22,  23] 

2=1  '  ' 

where  the  A,-  are  eigenparameters  of  the  parametrization  chosen.  This  relation  implies 
that  the  conductance  is  now  given  by  an  entire  distribution,  and  mean  conductance 
and  conductance  fluctuations  involve  now  an  average  procedure  over  the  distribution. 
The  most  striking  feature  in  theses  systems  is  the  appearance  of  universal  conduc- 
tance fluctuations  ([24]  and  experimentally  [25]),  which  results  in  the  prediction  of 
^  in  terras  of  the  dimensionless  conductance.  This  quite  fundamental  phenomenon 
actually  prompted  the  development  of  the  theory  mentioned  above. 

Quite  active  is  also  the  area  of  ballistic  quantum  junctions.  Advances  in  construct- 
ing microstructures  make  it  possible  to  experimentally  measure  mesoscopic  systems 


19 

in  the  ballistic  regime.  In  these  systems  the  scattering  approach  depicted  before  is 
usually  applied  in  form  of  the  scattering  matrix.  The  junction  is  described  by  a 
Hamiltonian  which  is  strongly  coupled  to  the  leads  [26].  The  appropriate  ensemble 
for  this  configuration  is  Dysons's  circular  ensemble  [27,  28].  It  allows  to  calculate 
the  variance  of  the  conductance  fluctuations  (resulting  in  ^)  as  well  as  to  predict  the 
hneshape  of  the  weak  localization  corrections  [29].  It  should  be  noted  that  the  present 
random  matrix  theory  can  only  account  for  chaotic  junctions,  whereas  the  nonchaotic 
junction  in  a  random  matrix  formulation  is  still  an  open  question.  The  difference 
in  ensembles  (from  gaussian  unitary/orthogonal  to  circular)  originates  from  the  fact 
that  in  disordered  conductors  the  quasi  one  dimensionality  is  present  in  the  multi- 
plicative law  2.7,  whereas  ballistic  micro  junctions  show  no  dimensional  properties 
and  are  therefore  sometimes  regarded  as  'zero-dimensional'  systems. 

This  is  quite  similar  to  the  ballistic  quantum  dot.  Quantum  dots  separate  the 
leads  by  tunnel  junctions.  The  aspect  of  Coulomb  charging  effects  is  present  here, 
which  determines  the  spacing  of  the  conductance  resonances  [30],  but  the  resonance 
amphtude  can  still  be  calculated  from  random  matrix  theory  [31]. 

Another  interesting  study  area  in  mesoscopic  physics  are  small  metaUic  particles. 
These  systems  consist  of  an  entire  ensemble  of  approximately  equal  sized  particles, 
usually  mounted  on  a  (insulating)  matrix.  This  system  can  then  be  probed  with 
external  fields,  and  the  response  can  be  measured.  The  basic  assumption  here  is 
that  we  practically  have  no  knowledge  about  the  particles,  which  makes  the  maxi- 
mum randomness  (and  therefore  random  matrix  theory)  applicable.  This  field  was 
pioneered  by  Gor'kov  and  Eliashberg  [32]  and  has  remained  an  active  area  (see  [33]). 
Of  interest  here  is  that  the  experiments  can  also  be  repeated  with  more  disordered 
metaUic  particles,  which  can  also  crossover  to  insulating  behavior. 


CHAPTER  2 
THE  TRANSITION  MODEL 


In  the  last  chapter  we  discussed  properties  of  standard  random  matrix  theories, 
which  are  now  an  accepted  method  for  determining  features  of  chaotic  systems  and 
the  metalhc  behavior  of  disordered  conductors.  This  chapter  introduces  a  generahzed 
ensemble  of  random  matrices,  based  on  a  maximum  entropy  ansatz,  which  is  suitable 
for  describing  the  transitional  region  of  chaotic  systems  and  disordered  conductors. 

2.1     Maximum  Entropy  Ansatz 

We  consider  in  the  following  a  real  symmetric  random  matrix  H  with  matrix  el- 
ements H,j.  In  order  to  determine  a  probabiUty  density  P{H)  over  the  space  of 
matrices  H,  it  is  necessary  to  'measure'  the  distance  between  matrices,  that  is,  a 
metric  has  to  be  defined.  The  natural  choice  is 

ds'  =  Tr{dHf  =  Y:^{dHuf  +  2^idHi,y  (2.1) 

The  volume  element  is  then  given  by 

dH  =  l[dHulldH,,  (2.2) 

The  metric  and  the  volume  element  are  invariant  under  real  orthogonal  transforma- 
tions and  the  addition  of  a  constant  matrix. 

The  importance  of  these  invariances  can  be  exemplified  by  the  similarity  to  the 
approach  taken  in  differential  geometry  which,  given  a  metric  c/^,^,  defines  a  length 
element  by 

^^^  =  Y^gutM{x)dxudx^j,  (2.3) 
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and  volume  element  by 

dV^\detg{x)\'''lldx^  (2.4) 

Both  measures  are  now  invariant  under  the  less  abstract  mappings  of  rotation  and 
translation  of  the  space  {x^}. 

The  following  discussion  is  restricted  to  the  orthogonal  case  (/?  =  1),  but  the 
generalization  to  the  other  cases  (unitary  and  symplectic)  are  straightforward,  as 
they  only  change  the  dimensional  arguments.  With  the  matrix  measures  we  can  write 
the  diiferential  probability  as  a  product  of  the  invariant  measure  and  the  probability 
density. 

dP{{H,,})  =  p{{H,,})dH  (2.5) 

The  probabihty  density  expresses  the  probability  that  the  set  of  matrix  elements  Hij 
hes  in  the  interval  [H,j,Hij  +  dHij].  The  parameter  space  is  apparently  A^(yV  -  1)- 
dimensional  (orthogonal  example). 

The  probability  density  is  now  transformed  to  a  N-dimensional  parameter  basis 
to  span  the  space  by  independent  components.  As  most  obvious  and  usual  choice  of 
parameters  serves  the  eigenvalue  basis,  while  the  radial  eigenparameters  of  the  transfer 
matrix  approach  to  conductivity  are  another  example.  Lets  call  this  N-dimensional 
parameter  set  {xi}.  The  probability  density  can  then  be  expressed  as 

dP{{x,})  =  J{{x,})w{x,)  n  dx,  (2.6) 

The  probability  density  w{xi)  Hi  dxi  gives  the  measure  of  the  .r,  in  the  N-dimensional 
parameter  space,  while  the  Jacobian  J{x,)  gives  the  measure  of  the  x,  in  the  N(N-l)- 
dimensional  matrix  space. 

The  non-compactness  of  the  matrix  space  (which  is  simple  to  see  by  construct- 
ing a  sequence  of  matrices  which  exhibits  in  the  limit  diverging  matrix  elements  or 
eigenvalues)  will  make  the  integration  over  that  space  diverge. 


22 

This  argument  delineates  wiiy  the  weight  function  if(xj)  has  to  decay  for  large  x,. 
Without  it,  the  first  principle  f  dP  =  I  cannot  be  obeyed.  However,  the  weight  factor 
w{xi)  is  not  specified  by  the  approach  so  far  and  has  to  be  chosen  appropriately.  In 
contrast,  given  a  certain  symmetry  of  the  matrix  H  and  a  unique  set  of  parameters, 
the  Jacobian  factor  J(xj)  can  uniquely  be  determined  (except  for  a  multiplicative  con- 
stant). Therefore,  the  statistical  properties  have  to  be  fixed  by  defining  an  additional 
condition  on  the  probability  density,  which  leads  to  the  determination  of  w{x,). 

In  the  derivation  of  the  probability  density  in  standard  random  matrix  theory,  the 
condition  was  to  restrict  the  average  trace  of  the  matrix  H.  This  definition  is  equivalent 
with  the  (rather  natural)  assumption  that  each  linearly  independent  component  is  also 
statistically  independent.  The  random  matrix  can  then  be  imagined  as  containing 
matrix  elements  which  obey  a  gaussian  distribution  and  which  lack  any  correlation 
between  each  other  (except,  of  course,  the  symmetry  requirement). 

Furthermore,  the  trace  depends  only  on  the  N-eigenparameters.  This  property 
ensures  the  zero  parameter  dependence  of  standard  random  matrix  theory  (after  the 
average  eigenparameter  separation  has  been  scaled  out).  Certainly,  this  assumption 
has  already  physical  relevance,  but  it  is  restrictive  for  other  physically  motivated 
matrix  element  correlations  [34,  35]. 

Therefore,  the  goal  is  to  find  a  more  general  way  of  fixing  the  statistical  properties 
of  the  random  matrix  approach.  The  Jacobian  or  interaction  term  J(x,)  provides 
already  the  level  repulsion,  but  leaves  the  eigenparameter  density  cr{x)  unspecified. 
This  leads  to  the  approach  to  fix  the  weigth  factor  iy(x;)  by  defining  the  density  a{x). 

This  idea  is  supported  by  the  results  from  experimental  and  numerical  studies. 
Often  the  density  (as  a  first  moment  of  the  distribution  P)  shows  system  dependency 
[36,  86,  37],  whereas  the  variance  (or  second  moment  of  P)  displays  universal  behavior. 

Since  a  general  theory  like  the  random  matrix  approach  is  most  suited  for  de- 
scribing universal  behavior,  it  is  plausible  to  take  the  system-dependent  density  as 
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an  input  parameter  to  determine  the  weight  factor  w{xi)  (and  therefore  the  confining 
potential).  Due  to  the  global  character  of  the  density  this  ansatz  is  often  called  global 
maximum  entropy  approach  (see  [21,  22]) 

For  transfer  matrices  M  a  different  kind  of  approach  is  possible.  It  originates  from 
the  multiphcity  law  transfer  matrices  obey.  Consider  a  system  (a  mesoscopic  wire  for 
example)  of  length  L  with  the  goal  to  add  a  small  piece  6L,  which  is  still  of  mesoscopic 
size.  The  associated  transfer  matrices  are  called  Ml  and  Msl-  Then  the  resulting 
transfer  matrix  is 

M  =  Ml+sl  =  MslMl  (2.7) 

If  the  two  pieces  are  statistically  independent,  the  probabihty  density  for  the  combined 
block  can  be  written 

PL+5l{M)  =  jpL{MslM)psL{M5L)dMsL  (2.8) 

The  measure  dM  is  the  invariant  measure  associated  with  the  transfer  matrix  ensem- 
ble. A  similar  integral  equation  appears  also  in  Markovian  processes  and  allows  the 
conversion  into  a  Fokker-Planck  type  differential  equation.  The  physical  restriction 
necessary  to  fix  the  statistical  properties  is  usually  chosen  to  assign  a  mean  free  path 
value  1  [39,  40]  This  restriction  leads  to  an  isotropic  probability  density,  which  is  also 
preserved  in  equation  2.8.  Then  one  obtains  for  5L  — >■  0  [] 

where  s  =  j  and  the  initial  condition  is  given  by 


w 


iP) 


S{X)  (2.10) 


This  equation  has  recently  been  exactly  solved  [41].  The  solution  shows  a  smaU  con- 
tribution from  a  second  two  body  interaction  term,  which  makes  the  global  approach 
for  transfer  matrices  just  an  approximation  to  the  problem,  albeit  an  excellent  one. 
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as  a  second  two  body  interaction  term,  is  not  present  in  the  giobal  formulation.  The 
same  problem  does  not  arise  for  the  eigenvalues  of  a  Hamiltonian  matrix,  since  the 
energy  levels  are  directly  the  eigenvalues  of  the  Hamiltonian. 

However,  it  is  important  to  mention  that  the  solution  of  2.9  has  a  potential  which 
displays  a  lognormal  dependence  for  large  arguments.  This  is  an  important  indication 
on  how  a  model  in  global  maximum  entropy  framework  should  be  constructed. 

2.2     Formulation 
2.2.1     Generalized  One  Body  Potential 

The  maximum  entropy  ansatz  determines  the  weight  factor  from  the  density.  Fur- 
thermore, the  weight  factor  can  be  determined  by  the  Lagrange  multiplier  method 
already  used  in  1.15.  In  general  one  has  the  functional  integral 

S[w{{x,})]  =  -  I  w{{x,})  \n{w{{x,}))d{x,}  (2.11) 

subject  to  the  constraining  conditions 

<  fn{{xi})  >-- an  (2.12) 

Here  i  runs  from  zero  to  the  number  of  constraints  and  <>  represents  the  expectation 
value 

^fa.TW    /nf^x,/({x.})P^({xJ) 

<  f{{^^\)   >=   r^N,       pu      ^^ (2-13 

J  n,    dx,Pp{{x,}) 
The  comparison  with  the  given  density 


iV 


(r{x)=<^5{x,-x)>  (2.14) 

gives  immediately  an  infinite  number  of  constraints  defined  by  a  function  n  — )■  x,  a„  — > 
V{x).  Lagrange's  procedure  results  then  in  the  variation 

5S  =-  f  dx6p{x)ln[p[x]]  +  1  -  /  V{x)  J2  ^{^  -  ^r)dx  =  0  (2.15) 
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obtaining 

N 
u;(xO-exp(-^y(xO)  (2.16) 

!  =  1 

Lets  now  concentrate  on  the  generic  form  of  the  Jacobian  in  equations  1.6,1.8,1.10, 
where  we  can  write  (again  including  the  symmetry  factor  (3) 

P{xi)    =    C0,Nexp{-l3H) 

N  N 

Cp^Nexp{^PJ2H\x,  -  Xj\)  +  J2V{x,))  (2.17) 

This  is  the  fundamental  result  for  the  JPD  from  the  maximum  entropy  ansatz.  The 
power  law  dependence  of  the  Wigner-Dyson  ensemble  is  replaced  by  an  arbitrary 
function  V,  which  is  determined  by  the  den^^ity.  The  generalization  of  Wigner's  re- 
sult 1.18  is  obvious.  In  addition,  the  potential  V  allows  for  parameter  dependence,  a 
possibility  which  is  not  present  in  the  original  standard  model. 

The  potential  can  be  approximated  in  the  continuous,  large  N  limit  by  a  mean  field  ar- 
gument, which  was  originally  used  by  Wigner  to  calculate  the  density  of  the  standard 
random  matrix  model  (the  famous  semi-circle  law). 

The  most  probable  configuration  of  the  eigenvalues  {.x,}  can  be  found  by  calculating 
1^  which  results  in 

dV      ^        I  ,        , 

-^  +  E|    ._     .1  =0  (2-18) 

By  ignoring  the  'self-energy'  of  each  eigenvalue  (the  eigenvalue  self-correlation),  and 
assuming  the  average  and  most  probable  configuration  to  be  sufficiently  close,  the 
potential  can  be  estimated  in  the  mean  field  sense  as: 


V{xi)=      a{x')ln{\x-x'\)dx'  (2.19) 

By  using  the  functional  derivative 

<  a{x)a{x')  >  -  <  a(x)  X  a{x)  >=  -^^||^  (2-20) 
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it  can  be  elegantly  demonstrated  that  the  variance  of  any  linear  statistic 


F  =  J  f{x)a{x)dx  (2.21) 

q  {^{^)  =<  Eit  K"^t  —  2;)  >  is  the  microscopic  density)  is  independent  of  the  potential  V 

;i  and  therefore  universal  [42].  This  underlines  the  success  of  the  random  matrix  theory 

■  constructed  by  means  of  the  global  maximum  entropy  approach. 

It  is  important  to  recognize  the  limitations  of  the  result  above.   A  necessary  con- 
dition has  been  the  ability  to  go  to  the  continuous  limit  in  the  case  N  -^  00.  This, 
however,  depends  on  the  strength  of  the  confining  potential.    A  'strongly'  confining 
potential  will  hmit  the  eigenvalues  close  to  the  origin,  and  therefore  in  the  A^  — )■  00 
j  limit  the  eigenvalues  will  move  very  close  together  allowing  a  continuum  approxima- 

',  tion.   In  case  of  a  'weak'  potential,  where  the  eigenvalues  move  towards  infinity,  the 

continuum  reasoning  does  not  apply.  The  role  of  weak  potentials  represents  a  central 

I  topic  in  this  investigation  and  will  be  pursued  below. 

I 

i 

2.2.2     The  Method  of  Orthorgonal  Polynomials 

,j  The  equilibrium  probability  density  function  2.17  requires  apparently  a  more  gen- 

j  eral  formalism  than  the  continuum  limit  approximation,  which  allows  to  cover  arbi- 

trary potentials  and  a  finite  number  of  eigenvalues.  This  formalism  is  the  method  of 
orthogonal  polynomials  and  was  first  used  by  Mehta  [4].  The  key  observation  consists 
of  the  fact  that  the  Jacobian  factor  Ot-^j  kz  -2;^^  (/5  =  2  for  simplicity)  can  be  written 

in  terms  of  a  Vandermonde  determinant 

2 


TV 


11  P'        -^J 


h3 


1  1         ..         1 

Xi        X2        ■■      Xn 


(2.22) 


By  using  properties  of  determinants  one  can  multiply  the  (j-l)th  row  with  an  arbitrary 
factor  and  add  it  to  the  jth  row.  Continuing  this  procedure  to  the  Nth  row  leads  to  a 
polynomial  of  the  jth  degree  in  the  jth  row.  Multiplying  each  row  with  the  potential 
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term  results  in 

N 


=    Adet{k,{x,)exp{--V{x,))  (2.23) 

where  kj{xi)  represents  a  polynomial  of  degree  j.  The  advantage  of  choosing  the  or- 
thogonal polynomial  set  ^^(xj)  which  corresponds  to  the  integration  measure  exp{  —  ^V{x))da 
for  the  polynomial  kj{x^)  is  obvious,  as  it  allows  to  evaluate  n-point  correlation  func- 
tions by  using  the  orthogonality  of  the  polynomials 

/•CO  I 

/       ri{x)rj{x)exp{--V{x))dx  -  SnSij  (2.24) 

J-oo  2 

[sn  is  the  normalization  of  the  polynomials  Tj).    The  density  (one  point  correlation 
function)  for  example  is  given  by 


/CO 
dx2..dxMP{xi) 
-oo 

N-l 

^     Y:rUx)exp{-V{x))  (2.25) 

n=0 

and  the  two  point  correlation  function  is 

/oo  /•GO 

..    /         dX3dxMP{xi,X2,..,XM) 
-oo       J —oo 


N  -,  N  1 

^     j:r]{x)exp{--V{x))Y:r',ix')exp{--V{x'))- 

i-O  ^  j=0  "^ 

(E  r.{x)r,{x')  exp(-iy(x)  -  ^^(x')))^  (2.26) 

As  a  simple  example,  the  standard  random  matrix  model  has  the  confining  poten- 
tial V{x)  —  x^,  which  leads  naturally  to  the  Hermite  polynomials  r,(x)  =  //,(x)  as 
orthogonal  polynomial  set. 

The  importance  of  this  method  lies  in  the  generality  in  which  a  potential  V(x)  can 
be  treated.  For  each  integration  measure  w{x)dx  an  orthogonal  polynomial  set  exists, 
making  the  solution  of  the  joint  probability  density  2.17  for  each  V(x)  in  principle 
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accessible.  However,  for  practical  purposes  it  is  better  to  use  a  potential  which  has  a 
well-known  and  recorded  associated  polynomial  set. 

2.3     Formulation  of  a  Spectral  Transition  Model 
2.3.1     Strongly  and  Weakly  Confining  Potentials 

The  maximum  entropy  ansatz  leads  to  the  relatively  independent  choice  of  the 
confining  potential,  which  is  only  fixed  by  the  density  configuration.  Remember  that 
this  additional  freedom  was  not  present  in  standard  random  matrix  theory.  The 
limitation  of  the  trace  of  a  random  matrix  is  straightforward  leading  to  a  harmonic 
oscillator  confining  potential. 

However,  the  trace  restriction  is  implied  by  statistical  independence  of  the  ma- 
trix elements.  They  obey  the  symmetry  condition,  but  otherwise  they  are  chosen  as 
randomly  as  possible.  This  'maximum  randomness'  approach  is  quite  useful  for  de- 
scribing physical  systems  from  which  we  have  minimal  information  such  as  completely 
chaotic  systems. 

However,  there  are  many  systems  which  do  not  belong  to  this  class.  The  matrix 
elements  usually  contain  further  correlations,  rendering  their  evaluation  more  difficult. 
The  global  maximum  entropy  ansatz  allows  to  use  a  measurable  quantity  like  the 
density  and  to  construct  then  a  reasonable  random  matrix  theory. 

As  already  mentioned  before,  the  confining  potential  can  fall  into  two  classes: 
weakly  and  strongly  confining.  By  imagining  the  eigenvalues  of  the  JPD  in  2.17 
as  two  dimensional  electrons  or  charges  which  repel  each  other  with  the  logarith- 
mic two  body  potential,  we  can  define  the  difference  between  these  two  confining 
potential  types.  In  the  strongly  confined  case  the  potential  is  approximately  given 
by  V{x)  =  f  a{x')ln{\x  —  x'\)dx'.  The  electron  density  a  has  to  be  neutralized  by 
the  potential  V,  otherwise  the  electrons  would  move.  Therefore  the  interpretation  of 
the  potential  as  to  arise  from  a  neutralizing  background  grid  of  positive  charges  is 
feasible.  This  grid  also  allows  the  electron  to  define  a  continuous  density,  even  though 
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the  electrons  repel  each  other.  In  the  case  that  an  additional  electron  is  added  to  a 
system  of  N  electrons,  an  additional  positive  charge  has  to  be  added,  too.  This  charge 
permits  the  electron  system  to  be  compact  and  not  to  expand  the  energy  coordinates 
to  infinity. 

The  situation  is  different  in  the  case  of  a  weakly  confining  potential.  It  cannot  be 
interpreted  as  a  neutralizing  background  charge,  because  the  electron  system  will  be- 
come larger  in  the  sense  that  the  difference  between  the  electrons  with  the  smallest 
and  the  largest  coordinate  grows  by  a  fixed  (N-independent)  amount.  For  infinitely 
many  electrons  the  position  coordinate  will  also  go  towards  infinity.  This  rationale 
allows  to  define  weak  and  strong  potentials,  depending  on  whether 

lim   <  {Max{{x,}))  >  (2.27) 

is  a  finite  limit  or  not.  A  local  version  of  this  condition  is 

f{xi)  =   lim   <  (x,  -  x,+i)  >  (2.28) 

going  to  zero  or  not.  The  finite  distance  between  charges  makes  a  continuous  de- 
scription, as  in  2.20,  not  exact  for  N  ^  oo.  This  is  why  the  orthorgonal  polynomial 
formalism,  which  is  explicitly  N-dependent,  is  more  appropriate  for  the  desired  tran- 
sition model. 

It  is  also  this  property  which  will  restrict  the  validity  of  the  universality  of  any  hnear 
statistics  in  equation  2.21  to  cases  where  a  continuous  limit  can  be  taken  in  an  appro- 
priate manner.  The  functional  derivative  reasoning  will  only  apply  for  a  continuous 
spectrum. 

There  have  been  already  investigations  concerning  the  question  how  strong  the  long 
distance  decay  behavior  of  the  potential  has  to  be  in  order  to  exhibit  transitional 
properties.  They  indicate  that  all  potentials  decaying  slower  than  V{x)  oc  x^^^  have 
a  different  qualitative  behavior. 
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2.3.2     Requirements 

The  method  of  orthorgonal  polynomials  enables  to  solve  for  the  solution  of  equation 

2.17  for  all  potentials  V(x),  which  does  not  violate  the  normalization  requirement 

/  P({x,})(i{x,}  =  finite.  To  construct  a  model  which  displays  transitional  properties 

in  its  spectral  statistics  therefore  means  to  choose  the  potential  V  appropriately. 

For  simplicity,  the  model  should  depend  solely  on  one  parameter.  This  is  an  obvious 

choice,  since  the  main  objective  of  the  discussed  model  is  to  capture  the  qualitative 

(!  features  of  a  transition  model.  For  physical  purposes  this  parameter  can  be  regarded 

I  as  a  measure  of  disorder  or  the  transition  parameter  for  chaotic  systems.   However, 

one  has  to  realize  that  an  exact  identification  of  this  model  parameter  with  specific 

l|  theoretical  models  is  difficult  due  to  the  semiphenomenological  nature  of  this  model. 

Even  more  important,  the  desired  model  should  pass  over  to  the  two  known  limits,  the 

Wigner-Dyson  random  matrix  model  and  the  Poisson  statistical  model.  This  requires 

I  the  model  to  have  a  weakly  confining  potential  at  large  arguments.    A  particularly 

:,i  good  choice  seems  to  be  the  lognormal  potential  K  a;  In  (x).    This  potential  in  the 

disordered  problem  occurred  first  in  the  one  dimensional  transfer  matrix  approach 

by  Melnikov  [43,  44].    Even  though  this  result  is  not  directly  comparable  with  the 

description  of  disorder  phenomena  via  energy  level  correlations,  it  still  represents 

a  cause  to  investigate  the  corresponding  transfer  matrix  related  q-Laguerre  model 

[45,  46].  This  model,  which  apparently  belongs  to  a  new  universality  class,  confirms 

the  unique  nature  of  the  weakly  confining  potential  including  translational  symmetry 

breaking. 

As  this  might  already  be  considered  sufficient  motivation  to  investigate  an  energy  level 
model  containing  a  weakly  confining  potential,  even  stronger  evidence  of  its  relevance 
comes  from  several  numerical  studies  [22,  47,  48,  49],  which  compare  a  similar  model 
with  specific  disordered  or  chaotic  systems.  Here  the  3~dimensional  tight  binding 
Anderson  model  should  in  particular  be  mentioned.    Canali  [47]  shows  quantitative 
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agreement  between  a  model  with  weakly  confining  potential  and  the  3D  Anderson 
model  at  the  critical  point.  For  example  the  linear  part  of  the  number  variance  and 
the  nearest  neighbour  spacing  distribution  show  good  agreement. 

The  standard  random  matrix  models  incorporate  a  strongly  confining  potential. 
The  transition  potential  should  have  a  similar  behavior  for  a  certain  value  of  the 
parameter.  While  the  weakly  confining  part  has  to  occur  at  large  arguments  x,  the 
Wigner-Dyson-hke  strongly  confining  part  of  the  potential  arises  at  small  arguments 
X.  Appropriate  for  describing  energy  levels,  the  limiting  potential  has  a  square  de- 
pendence. This  type  of  model  can  be  formulated  in  terms  of  orthogonal  polynomials 
with  the  Hermite  ensemble.  Therefore  the  other  limit  of  the  desired  potential  is 

YimV{x^{))-^x^  (2.29) 

It  is  obvious  that  this  choice  has  to  depend  on  a  parameter  in  order  to  display  a 
change  in  statistical  behavior.  The  potential  furthermore  requires  the  recovery  of  the 
Wigner-Dyson  limit,  for  example  a  square  dependence  for  the  Hamiltonian  model. 

2.4     Transition  Ensembles 

After  having  formulated  the  conditions  on  a  transition  model,  an  appropriate  type 
of  potential  has  to  be  chosen.  The  problem  consists  in  finding  a  potential  which  corre- 
sponds to  a  well  studied  polynomial  set  allowing  the  use  of  the  method  of  orthorgonal 
polynomials.  Furthermore,  this  potential  should  depend  on  one  parameter.  The  first 
approach  in  this  respect  was  proposed  in  [50]  with  the  potential  being 

V{x)^  -a\n\l  +  bx)  (2.30) 

It  should  be  noted  that  this  ansatz  has  two  parameters,  whereas  the  other  two  con- 
ditions are  certainly  satisfied.  This  ansatz  was  justified  by  conducting  a  numerical 
analysis  to  the  tight  binding  Anderson  model.  The  result  was  a  quite  reasonable 
agreement  in  the  correlation  functions  [50].    From  this  perspective  it  was  natural  to 
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investigate  the  q-Laguerre  ensemble  [45,  46]  as  an  attempt  to  model  the  transitional 
properties  of  transfer  matrices,  which  constitutes  a  direct  model  for  the  transport 
properties  of  disordered  conductors.  As  mentioned  before,  this  is  an  exact  solution 
which  allows  precise  predictions  from  the  formulation.  The  one  body  potential  is 
linear  in  x  for  small  x  and  shows  a  lognormal  dependence  for  large  x. 

The  properties  of  this  model  exhibited  a  decay  behavior  in  the  density  for  a  larger 
disorder  parameter  and  a  rapidly  changing  density  close  to  the  origin.  In  terms  of  the 
conductance,  the  model  was  found  to  be  in  qualitative  agreement  with  the  disordered 
conductor  physics,  i.e.  the  dependence  on  the  disorder  parameter  exhibited  both 
metallic  and  insulating  behavior. 

However,  the  small  discrepancy  between  this  model  formulated  in  the  global  entropy 
assumption  and  the  results  from  the  local  entropy  ansatz  (this  becomes  most  clear 
from  the  occurence  of  an  additional,  albeit  small,  two  body  interaction  term)  makes 
a  further  quantitative  agreement  disputable. 

Another  aspect  of  transitional  models  is  the  breaking  of  the  invariance  under  uni- 
tary transformations  or  U(N)-symmetry.  This  symmetry  is  inherent  to  standard  ran- 
dom matrices,  but  any  transitional  statistics  will  have  to  break  the  unitarity.  This 
is  immediately  obvious  from  the  fact  that  the  standard  random  matrices  have  no 
inherent  scale  other  than  the  number  of  eigenvalues.  In  disordered  conductors,  for 
example,  the  standard  random  matrix  theory  describes  well  the  eigenstates  of  the 
metallic  phase,  which  are  all  extended.  In  this  phase  one  set  of  eigenstates  describes 
the  whole  system  completely,  i.e.  every  other  eigenstate  set  can  be  found  by  linear 
combination.  In  the  insulating  phase  the  eigenstates  are  localized  and  therefore  of 
completely  different  nature.  Apparently,  the  'featureless'  standard  random  matrix 
theory  cannot  decribe  this  fundamental  transition. 

This  aspect  was  pointed  out  in  [53,  54]  along  with  the  argument  that  this  symmetry 
breaking  is  spontaneous  in  weakly  confining  potentials.   Other  approaches  recognize 
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the  symmetry  breaking  directly,  by  defining  a  random  matrix  theory,  whose  U(N)- 
symmetry  is  broken  by  construction.  This  approach  was  favored  in  reference  [55], 
where  JPD  had  the  following  form: 

P{H)  =  Cn  exp{-TrH^  -  h^N^Tr{[A,  H][A,  H]^))  (2.31) 

The  model  introduces  explicitly  an  unitarity  breaking  term,  which  depends  on  one 
parameter  h,  measuring  how  strongly  the  symmetry  is  broken.  The  matrix  A,  which 
itself  is  unitary,  indicates  a  preferred  direction  in  matrix  space  for  the  Hamiltonian. 
The  idea  behind  this  random  matrix  model  with  preferential  basis  is  to  average  over 
the  matrix  A.  The  model  has  the  corresponding  Wigner-Dyson  limit  for  h  =  0.  The 
interesting  result  from  this  theory  is  that  many  features  of  the  preferential  basis  model 
are  very  similar  to  the  weak  confinement  model  treated  in  here.  In  particular,  the 
two-level  correlation  function  (for  hN"^  —  const.)  in  the  critical  regime  is  identical 
to  the  one  in  the  transitional  model  treated  in  here  and  derived  in  the  next  chapter. 
This  relation  will  be  evaluated  in  more  detail  in  the  discussion. 

Another  way  of  breaking  the  symmetry  is  to  introduce  banded  random  matrices. 
The  physical  intuition  behind  this  model  is  the  notion  that  the  further  the  matrix 
element  is  away  from  the  main  diagonal,  the  smaller  it  becomes.  A  matrix  element 
on  the  main  diagonal  can  be  interpreted  in  terms  of  the  degree  of  the  overlap  between 
the  state,  in  this  case,  with  itself.  These  models  have  also  one  parameter,  commonly 
called  matrix  bandwidth  b,  which  measures  how  fast  the  matrix  elements  decay  with 
increasing  distance  from  the  main  diagonal.  The  transitional  parameter  is  usually 
cited  as  b/N,  with  N  being  the  number  of  eigenvalues  [35,  56,  58,  61].  So  far,  most 
investigations  point  into  the  direction  of  this  model  being  most  feasible  in  quasi-lD 
[59,  60]. 
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2.5     Applicability  to  Physical  Systems 

One  possible  application  of  a  transition  random  matrix  theory  is  to  understand  the 
energy  level  behavior  in  disordered  systems.  An  example  system  are  small  metalhc 
mesoscopic  particles. 

It  would  represent  a  particularly  important  step,  if  the  metal-insulator  transition 
was  better  understood  in  terms  of  a  random  matrix  formulation.  In  the  metalhc 
regime  the  Wigner-Dyson  random  matrix  theory  applies.  The  system  is  in  its  diffu- 
sive state,  where  electrons  lose  their  memory  on  the  scale  of  the  mean  free  path.  As 
long  as  the  mean  free  path  is  smaller  than  the  sample  size  L,  the  electrons  are  diffusing 
over  time  through  the  entire  sample.  This  'ergodicity'  defines  a  sample  characteristic 
time  and  energy  scale,  tc  and  Ec,  where  the  energy  scale  Ec  is  usually  called  Thouless 
energy.  The  Thouless  picture  connects  the  static  energy  spectrum  of  the  object  with  a 
transport  property.  However,  this  picture  breaks  down,  when  the  second  energy  scale, 
the  meanlevel  spacing,  becomes  comparable  to  the  Thouless  energy.  When  the  diffu- 
sive regime  is  crossing  over  to  the  ballistic  regime,  the  electrons  will  also  experience 
the  sample  edges  (boundary)  and  therefore  will  know  about  the  dimensionality.  The 
usual  standard  random  matrix  ensembles  do  not  provide  any  information  about  the 
dimension.  This  makes  them  only  suitable  for  describing  the  ergodic  metalhc  regime 
[57].  From  this  discussion  it  becomes  apparent  that  a  new  type  of  random  matrix 
theory  is  required  for  describing  the  statistical  properties  of  disordered  systems.  The 
interest  has  increased  recently,  as  it  was  suspected  that  at  the  critical  point  in  3D  the 
statistical  fluctuations  are  scale  invariant.  This  led  to  speculations  that  the  behavior 
at  the  critical  point  is  governed  by  a  unique  statistics,  i.e.  neither  Wigner-Dyson  nor 
Poisson  [63,  64].  This  renders  the  study  of  the  transition  model  very  interesting. 

Chaotic  systems  represent  another  area  of  physical  interest.  The  term  chaotic 
systems  should  be  defined  more  precisely.  We  are  usually  interested  in  quantum 
systems  whose  classical  counterpart  is  chaotic.     This  correspondence  is  necessary. 
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because  the  definition  of  cliaos  in  quantum  systems  is  not  straightforward,  whereas 
the  Ljapunov  exponents  give  a  precise  meaning  to  classical  chaos.  A  regular  classical 
system  describes  the  particles  as  moving  in  phase  space  on  invariant  tori  [65]  in  closed 
orbits.  Tiie  transition  to  chaos  is  then  seen  by  the  breakup  of  the  invariant  tori  by 
nonlinear  resonances  [8].  This  transition  as  a  function  of  a  parameter  will  initially 
start  only  in  parts  of  the  phase  space  and  will  gradually  proliferate  over  the  entire 
phase  space  to  complete  chaos.  On  this  basis  we  can  understand  the  terms  soft  and 
hard  (complete)  chaos. 

The  quantum  picture  is  more  complex,  as  particles  can  tunnel  from  their  present 
state  into  another  state  rendering  the  picture  above  inaccurate.  A  pragmatic,  albeit 
somehow  nonquantitative,  approach  to  quantum  chaos  uses  the  statistical  properties 
as  a  measure  of  how  chaotic  the  system  is.  As  already  mentioned,  the  statistical 
properties  of  regular  and  chaotic  systems  show  drastically  different  behavior.  The 
statistical  properties  are  also  easily  available  for  both  numerical  and  experimental 
sources. 

Standard  random  matrix  theory  has  been  very  successful  in  explaining  the  statis- 
tical properties  of  chaotic  systems  [66,  8],  while  directly  measurable  quantities  have 
been  more  difficult  to  bring  into  accordance  with  random  matrix  theory.  An  exam- 
ple is  the  energy  level  density  in  heavy  nuclei,  which  is  not  as  flat  as  the  gaussian 
orthorgonal  predicts  [9]. 

The  large  number  of  transitional  systems  makes  a  theory,  which  is  capable  of  pre- 
dicting the  transition  from  regular  to  chaotic  behavior,  desirable.  Some  systems,  like 
billiards  for  example,  are  either  integrable  (regular)  or  nonintegrable  (chaotic).  The 
transition  between  these  two  limits  is  usually  abrupt,  and  a  parameter  which  de- 
termines this  transition  is  not  always  obvious.  Furthermore,  billiards  usually  exhibit 
one  or  more  spatial  characteristic  length  scales,  the  shortest  periodic  orbit  being  most 
important. 
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Such  peculiarities  in  a  system  are  usually  difficult  to  explain  by  random  matrix 

theories.  But  others,  for  example  the  hydrogen  atom  in  a  strong  magnetic  field, 
show  transitional  behavior  as  function  of  a  parameter  (  the  magnetic  field).  It  can 
experimentally  be  generated  with  relative  ease,  and  the  Hamiltonian  equations  can 
be  analyzed  by  numerical  means.  Its  statistical  behavior  is  very  similar  to  the  weakly 
confining  potential  theory  [67].  For  this  class  of  systems  it  is  quite  important  to 
develop  a  theory  which  can  characterize  the  transitional  features  of  chaotic  systems. 
The  complexity  of  analytic  results  of  chaotic  systems  makes  the  comparison  of  the 
statistical  behavior  an  essential  part. 

Another  class  of  systems,  which  could  be  better  understood  by  looking  at  the  tran- 
sitional model,  are  quantum  systems  depending  on  external  parameters  [68,  69,  70, 
71,  72,  73].  The  external  parameter  is  physically  quite  different  from  the  transitional 
parameter  discussed  before.  The  transitional  parameter  is  connected  to  the  disorder 
in  the  system  and  represents  an  internal  physical  quantity,  which  is  fixed  for  a  given 
system.  One  needs  an  entire  ensemble  of  systems  in  order  to  get  meaningful  (ensem- 
ble averaged  or  disordered  averaged)  results.  In  contrast,  the  external  parameter  can 
by  varied  from  the  outside  and  for  each  given  system  independently.  There  are  quite 
many  systems,  which  fall  into  this  category.  Just  a  few  examples  are  mesoscopic  rings 
with  internal  magnetic  fiux,  mesoscopic  systems  with  external  electric  or  magnetic 
field  and  also  quantum  dots  not  in  complete  chaos  with  the  presence  of  an  external 
field. 


CHAPTER  3 
THE  DEFORMED  GAUSSIAN  ENSEMBLE 


As  outlined  before,  this  chapter  investigates  the  model  based  on  the  ensemble  of 
q-deformed  Hermite  polynomials  or,  in  short,  q-Hermite  polynomials.  This  model  is 
particularly  interesting,  since  it  provides  additional  insights  for  the  behavior  of  energy 
levels  in  disordered  and  chaotic  systems. 

3.1      Formulation 

The  general  theory,  described  in  the  second  chapter,  allows  to  define  a  unique 
random  matrix  model  by  specifying  an  orthorgonal  polynomiial  system.  Given  a 
complete  polynomial  system,  a  weight  function  (or  invariant  measure)  can  be  found 
which  corresponds  to  the  one  body,  confining  potential.  In  case  of  the  q-Hermite 
ensemble  the  potential  is  given  by 

1  1 
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^ai,a2,..,a„;(?)cx,  =  {auq)^{a2;  q)^  ■  ■(a„;g)co  (3.4) 


The  formulation  in  this  chapter  is  expressed  in  terms  of  the  general  variable  x.   For 
physical  interpretation  the  variable  x  can  be  replaced  by  the  energy  coordinate  E. 
In  order  to  ensure  the  suitability  of  this  potential  for  the  description  of  transition 
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energy  level  statistics,  three  conditions  have  to  be  met.  The  first  one,  the  dependence 
on  one  parameter,  is  immediately  obvious. 

However,  the  other  conditions  have  to  be  checked  explicitly,  with  regard  to  a  stan- 
dard square  dependence  of  the  potential  at  small  coordinate  x  and  a  lognormal  depen- 
dence of  the  potential  at  large  coordinate.  Jacobi's  triple  product  formula  provides 
one  way  to  verify  this. 

{q'-,q')oo{^zq-q')^[~'^--q')^=     Y,    q-\-  (3.5) 
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n=— oo 


Choosing  z  =  e'^^  we  arrive  at 
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The  sum  can  be  approximated  by  an  integral.  Completing  the  square  in  the  exponent 
results  in 

1  rco  , 

— —     =     /      e-^"+'""  (3.7) 
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where  7  is  related  to  q  by  <?  =  e"'^.  Scaling  the  coordinate  (x  -^  y^x)  is  eliminating 
the  problematic  7  — )•  0  limes.  The  joint  probability  distribution  function  is  then  given 
by 

P{{xn})  =  \xi  -  Xj\'^w{x)  (3.9) 

A  transformation  to  variable  sinh  x  —  x  gives 


zi;(x)  =  J(x,x)V27r7e-"^"'"''W  (3.10) 

where  the  Jacobian  J(x,  x)  has  to  be  independent  of  the  potential.  Now  the  expansion 
of  the  arcsinh  function  for  small  and  large  argument  x  shows  the  desired  requirements. 
Therefore,  the  q-Hermite  ensemble  is  a  suitable  choice  for  a  transition  statistics  model. 
The  general  theory  allows  to  write  the  JPD  as 

P{x,)  =  det\[Uxk)]'^[Mxj)]\  (3.11) 
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with  the  4^i{xj)  representing  the  eigenfunctions  ('single  level  wave  functions')  of  the 
model  given  by 

4>n{Xm)   =    —hn{Xm]q)\/w{x)  (3.12) 

where  the  polynomials  hn{xjn)  are  defined  by  the  recursion  relation 

1  -  g" 
hn+\{x]  q)  =  2xhn{x]  q) — hn-i{x;  q)  (3.13) 

and  the  relations 

ho{x-q)     -     1  (3.14) 

hi{x]q)     =    2x  (3.15) 

The  kn  are  the  normalization  factors  given  by  the  orthogonality  relation  [76] 

/oo 
hn{x]q)hmix;q)w{x)dx     =     7(5;  9)oo(<?;  <?)ng~"^"~^^'^<^nm  (3.16) 
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These  equations  define  the  q-Hermite  random  matrix  ensemble. 

In  the  mathematical  literature  these  polynomials  are  called  the  ^"^-Hermite  polyno- 
mials which  were  widely  investigated  [74,  75].  There  are  actually  infinite  possibilities 
to  generalize  the  q-Hermite  polynomials.  However,  it  is  important  to  note  that  even 
though  the  generalization  is  not  unique,  the  displayed  qualitative  behaviour  of  each 
ensemble  after  rcscaling  is  the  same.  This  rescaling  with  the  density  is  necessary,  as 
can  be  seen  from  the  the  difference  in  the  level  density  between  q-Hermite  and  scaled 
q-Hermite  polynomials,  which  we  will  discuss  later.   . 

3.2     Correlation  Functions 
3.2.1     Level  Density 

The  eigenvalue  level  density  is  then  easily  expressed  in  the  form 

Af-l 

aN{x;q)     =     ^  <^'(^)  (3-18) 

n=0 
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In  the  case  N   ^   oo  the  density  can  be  calculated  in  closed  form  (compare  the 
calculation  of  the  kernel)  resulting  in 

a^{x;q)  =  —  (3.20) 

27 

Here 

7  =  ln(l/<?)  (3.21) 

is  a  convenient  redefinition  of  the  deformation  value  q.  For  finite  N  the  density  has 
to  obey  the  sum  rule  J  aj^dx  =  A^,  and  therefore  the  density  has  to  vanish  for  large  x. 
This  is  confirmed  by  numerical  evaluation  shown  in  figure  3.1.  Visible  is  the  typical 
behaviour  of  a  relatively  constant  density  close  to  the  origin.  When  the  coordinate 
X  approaches  N,  the  integral  displays  a  relatively  sharp  drop  to  zero.  The  finite  case 
is  discussed  in  more  detail  later  in  this  chapter.  Hov^^ever,  the  value  of  the  density  at 
the  origin  varies  as  a  function  of  q.  Apparently  the  limit  <?  — >  1  cannot  be  obtained 
with  these  polynomials  (as  the  density  diverges),  but  rescaling  of  the  polynomials  /i„ 
(and  consequently  the  potential  w)  into  the  /i„  provides  a  well  defined  density  limit. 

7"/^M^;?)     =     hr.{^x;q)  (3.22) 


w[x)    =    ■w{y/^x)  (3.23) 

Alternatively  one  can  use  the  modified  recursion  relation 

-  -  1  -  q"- 

hn+i{x;  q)  =  2x/i„(x;  q) — /i„_i(x;  q)  (3.24) 

79 
The  results  for  the  density  of  the  scaled  polynomials  w{y^x)  J2n  ^n(^i  *?)  ^^^  displayed 
for  finite  N  in  figure  3.2.    Here  the  smooth  crossover  from  the  q-Hermite  ensemble 
[q  <  1)  to  the  Hermite  ensemble  is  obvious,  showing  Wigners  semi-circle  law  for  finite 
N... 

The  slower  decay  of  the  density  can  be  attributed  to  the  weaker  potential.  A 
strongly  confining  potential  (for  example  the  usual  square  potential)  keeps  the  eigen- 
values close  to  the  origin.  However,  in  a  weakly  confining  potential  3.7  the  eigenvalue 
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0.25 


Figure  3.1:  Eigenvalue  density  for  the  unsealed  q-Hermite  ensemble  for  N=10  at 
different  7  Here  the  solid  line  represents  the  density  for  7  =  0,  the  dashed  line 
corresponds  to  7  =  0.04,  the  dashed-dotted  line  to  7  =  0.1  and  the  dotted  line  to 
7  =  0.69 
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Figure  3.2:  Eigenvalue  density  for  the  scaled  q-Hermite  model  for  N=10  at  different 
7.  Here  the  solid  line  represents  the  density  for  7  =  0,  the  dashed  line  corresponds 
to  7  =  0.04,  the  dashed-dotted  line  to  7  =  0.1  and  the  dotted  line  to  7  =  0.69 
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repulsion  is  considerably  higher  and  therefore  the  eigenvalues  are  much  further  apart. 
Hence  the  level  density  is  much  smaller  at  the  origin.  This  imphes,  as  the  integral 
over  the  density  is  constant,  that  the  dropoff  of  the  curve  is  less  steep. 

The  density  at  the  origin,  as  it  is  fairly  constant,  provides  thus  a  measure  for  the 
mean  level  spacing  of  the  ensemble  and  further  rescaling.  It  can  be  determined  with 
the  help  of  the  Christoffel-Darboux  formula  applied  to  q-Hermite  polynomials. 

^w{x)w{y)  2_^  Hn[x]  q)Hn[y]  q)  = (3.25) 

n=l  '^n  X         y 

by  taking  the  limit  x  ^  y.  Here  the  //„  represent  the  normalized  q-Hermite  polynomi- 
als. Recognizing  that  odd  and  even  q-Hermite  polynomials  have  a  different  behavior 
at  the  origin  we  define 

An  =  limh2nix;q)  (3.26) 


and 


Calculating  A„ 


n        ,.      -^^271+1(3;;  q)  ,        - 

Bn  =  Urn (3.27) 


x->-0 


1  _  Q  1  _  o3         I  _      2n-l 

^"- -^■-       2!-!  (3-28) 

q  q-^  q^^    ^ 


and  similarly  Bn 


1  -q2" 

h2n+i{x^Q]q)    =    2xh2n{x;q) — h2n-i{x;q)  (3.29) 

qzn 

=     -2x -^-— /i2n-2(0w)  -  2a: ^h2n-2{0]q)         (3.30) 

+     -^2 ^^--3(0W)  (3-31) 

=    ...  (3.32) 

1  -  f?^"-!  1  -  g^"-^       1  -  (7^"  1  -  (72"-3      l-q 

=   -2x{^-^^ r-^  +  — r r-, ^3.33 

^        Qin  —  l  Qin  —  i  gin  Qin  —  .i  q       ^  ' 

+  ■•  + :f 7^ ^ J-)  3.34 

q2n         q2n~2  q^         q^ 

=     ~2xAn  (3.35) 
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we  get 

a-2iv(0;<?)    = ^ nin-i)i2i ^^4„5„  (3.36) 


V^    ^  TT  " 


This  expression  can  be  well  approximated  by 


E^^n^T^  (3-37) 

k=\    3=1  ^    y 

1        ^        ,  '    l-exp(-7(2j-l) 


cx..(0;,)«|,^     ;j^^^  (3.39) 


A  smooth  approximation  is  given  by 


a,N{Q-q)^^-^  (3.40) 


'? 


where  Nq  —  ^fr"-  For  7A^  <<  1  the  expression  3.40  goes  to  the  known  Wigner 
semicircle  law  (except  for  the  scahng  parameter  -^).  Also  the  Hmit  7A^  >>  1  is  re- 
covered as  required  by  equation  3.20  The  density  will  be  important  for  understanding 
the  scaling  properties  discussed  later. 

3.2.2     The  Kernel 

In  the  method  of  orthorgonal  polynomials  two  level  properties  as  the  number  vari- 
ance and  the  As-statistic  are  usually  calculated  with  the  help  of  the  kernel.  The 
kernel  depends  in  a  simple  form  on  the  eigenfunctions  3.12. 

K{x, y]q)  =  Yl  ^^(^'^ Q)4>n{y'^ "?)  (3-41) 

n-l 

The  kernel  for  the  q-Hermite  polynomials  can  be  exactly  calculated  for  A^  =  oo,  which 
was  first  shown  in  Ismail  [76].  For  this  calculation  one  uses  the  following  q-generalized 
linearization  formula  usually  used  to  calculate  a  product  of  two  polynomials  as  it 
appears  in  the  kernel. 

^m(m-l)/2^m(j;9)^n(n-l)/2^n(x:(?)    ^  ^^^^^ 

{q]q)m  {q;q)n 
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minin^m)      -n-m  +  k{k+l)/2+{m-k+l)(m-k)/2+(n-k+].)(n-k)/2 

Yl     7 TT ^ ? ^ hm+n-2k{x;q)     (3.43) 

k=o  {q;q)k[q;q)m-k{q]q)n-k 

which  expresses  the  product  of  two  q-Hermite  polynomials  as  a  linear  combination  of 
q-Hermite  polynomials.  This  linear  combination  is  a  generalization  of  the  generating 
function 


oo      ...^ 


^n    n(n-l)/2 


W 


E  -7—^ hn{x;  q)  =  {-t{x  +  Vx^  +  l);  t{Vx^  +  1  -  X-  q)^  (3.44) 

n=o   v"?'  q)°° 

hich  leads  to  (the  details  are  provided  in  [91] 

»  n(n-l)/2 

J2  hn{smh{x);q)hn{smh{y);  q)- ^BJ"  (3.45) 

n=Q  [q\q)n 

-  {-qRe'^+y-  -qRe-^'-y-  qRe^-y-  qRe^-^-  qU I {qR^ ]  q) 00  (3.46) 

This  formula  gives  the  kernel  (take  R=l)  for  A^  =  00 

Kix,y;q)    =     — — —  ^ /i„(a;;  g)/i„(j/;  g)  (3.47) 

w(x)w(y) 
^    ^     ^    ^'     qe=^+y--qe-^-y-qe^-y-qe-^+y-q)^  (3.48) 


The  infinite  products  can  be  further  expressed  in  Jacobi  theta  functions.    Using  the 

following  relations 

00 

{qe^-y;qe-^^-y^-q)     =     Y[{1  -  2q- coshix  -  y)  +  q'^)  (3.49) 

1  9M^-y)/2-,^) 

2q^l\q-qU   ^  sinh((x  -  y)/2)  ^"'^^^ 


anc 


-qe^^y--qe-(-+y)-q)    ^     fj  (1  +  2g"  cosh(x  +  y)  +  9^")  (3.51) 


n-l 


1  e^{i{x  +  y)l2;^) 

2q^l%q-q)^   2  cosh((x  +  y)/2) 

the  kernel  can  be  expressed  as 


(3.52) 


(7r/7)i/2  cosh^/'(x)cosh^/'(y) 

^ln(l/(?)[((?;9)c«]3      cosh((x  +  y)/2) 


e,{i{x  +  y)l2-^)       9,{i{x-y)/2;^) 


(3.54) 
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Since  the  imaginary  arguments  are  difficult  to  evaluate,  Jacobi's  imaginary  transfor- 
mations are  applied 

9,{ix-y^)     -     (^/7)i/^e^'/-^^4(7r.x/7,e-2^'/"')  (3.55) 

9,{ix;^)     =    z(7r/7)i/2e^'/"^^i(^x/7;e-^^'/^)  (3.56) 

resulting  in 

(^/7)'/'  cosh^/'(x)cosh'/'(y) 

2ln{l/q)[{q;q}^\^       cosh{{x  +  y)  /  2) 

94{n{x  +  y)/2r,  e"'"'/-^)  dyJTTJx  -  y)/2r,  e'^^'h) 

e\'\-Kxh]  e-^-'h)e\'\^yh\  e-^-'h)         sinh((x  -  y)/2) 

In  the  limit  7  <<  2TT^,y  Ri  x  one  can  further  simplify  the  result  by  approximating 

the  theta  functions  by 

e,(z,q)  =  2q'^Um{z)  (3.58) 

e^{z,q)  =  l  (3.59) 

finally  obtaining  (and  using  r,  .  \   ,3  —  exp(7r^/47)(^)^/^) 

1  cosh^''^(a;)  cosh^/^(y)  sin(7r(x  -  v)/27) 

K{x,y;q)  =  - y  ^^^     .\) ^^  (3.60) 

Ztt      cosh[[x  +  y)/2)       smh((x  —  y)/2) 

This  formula  represents  a  fundamental  result  for  the  q-Hermite  model  [77].  In  order 

to  concentrate  on  the  bulk  behavior  the  assumption  a:-y<<l,x-}-y>>l  can  be 

imposed,  which  simplifies  the  kernel  to 

r^/  ,         1   sin[7r(x  —  v)/27l 

K{x,y;q)=~ r\}z %^  3.61) 

Ztt  smh((a;  —  y)/2) 

This  kernel  has  apparently  the  density 

K{x,x-q)^a{x;q)=  ~  (3.62) 

which  is  constant  but  not  unity.  Fortunately,  the  transformation  which  leads  to  unit 

density  and  therefore  an  unfolded  spectrum  is  obvious.   With  x  -^  x/-f,y  — )■  y/7  we 

arrive  then  at  the  scaled  bulk  kernel 

-— ,  ,         1       sin(7r(x  —  y)) 

K{x,y-q)^—--y ^j-  (3.63) 

ZTT  smh(7(x  —  yj/2j 
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In  order  to  interpret  this  result  the  reduced  two  level  probabihty  function  R2  =  I  —  K 
as  a  function  of  the  distance  z  =  x  -  y  \s  plotted  for  different  7  in  figure  3.3.   Note 
that  only  the  square  of  the  kernel  has  a  probability  meaning  and  is  hence  physically 
relevant. 

The  resulting  limiting  cases  are  easily  understood  by  noting  that  in  the  Wigner- 
Dyson  Hmit  the  level  repulsion  prohibits  levels  close  to  each  other.  Accordingly,  the 
two  level  probability  function  is  zero  when  the  distance  between  the  two  levels  is 
negHgible.  When  the  levels  are  far  apart  they  don't  influence  each  other,  and  the 
two  level  probabihty  function  is  nonzero  (scaled  to  one  in  the  plot).  On  the  other 
hand,  if  the  eigenvalues  behave  in  accordance  to  Poisson-statistics,  they  are  entirely 
independent  from  each  other,  leading  to  a  constant  two  level  probability  function  even 
at  small  argument. 

These  limiting  cases  are  well  known,  however,  the  qualitative  form  of  the  transi- 
tion between  Wigner-Dyson  and  Poisson-statistics  is  characteristic  for  the  q-Hermite 
model  and  can  be  compared  with  either  experimental  data  or  numerical  experiments. 
The  power  law  decay  of  the  Wigner-Dyson  statistics  (obvious  in  the  Kernel  Kw-d{^,  y) 
j^sin_(7r(x-y]2^  is  replaced  by  the  exponential  decay  of  the  q-Kernel().  Note  that  the 
exponential  decay  is  present  for  all  7  >  0,  indicating  a  fundamental  difference  between 
the  Wigner-Dyson  limit  and  the  extension  to  g  <  1.  This  behavior  represents  one 
example,  where  it  becomes  visible  that  the  extended  (q-deformed)  random  matrix 
models  belong  to  a  different  universality  class  of  random  matrix  theories  compared 
to  'classical'  random  matrix  theories.  In  fact,  it  can  be  expected  that  all  q-deformed 
models  give  the  same  universal  behavior,  thereby  indicating  the  same  universality 
class,  because  their  one  body  potential  has  always  a  lognormal  behavior  for  large 
eigenvalues.  This  assertion  can  be  supported  by  investigating  another  q-deformed 
ensemble,  and  showing  that  the  kernel  is  quantitatively  the  same  in  the  bulk.  This 
has  explicitly  been  done  for  the  q-Laguerre  ensemble  [45,  46]. 


48 


Figure  3.3:  Reduced  two  level  probability  function  for  A^  =  oo  and  different  7.  Here 
the  solid  line  represents  7  =  0,  the  dashed  line  represents  7  :^  5  and  the  dashed-dotted 
line  represents  7  =  20. 
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Because  of  the  importance  of  this  argument  and  the  brevity  of  the  calculation  in 
the  reference  mentioned,  we  will  outline  the  calculation  here  in  more  detail.  The 
q-Laguerrc  ensemble  is  defined  by  the  following  power  series: 

Ln[x;  g)  -  2^ 7—7^ (3.64) 

This  definition  requires,  as  before  for  the  q-Hermite  polynomials,  the  knowledge  of  the 
recursion  relation  as  well  as  the  normalization.  Additionally,  the  invariant  measure 
is  required. 

wl{x]  q)  =  —— ^ —  (3.65) 

(-(1  -q)x;q)^ 

It  was  shown   [46],   that  this  invariant  measure  fulfills  the  requirements  for  a  q- 

generalized  model  (strong  confinement  at  the  origin  and  weak  confinement  at  large 

coordinate).  The  Christoffel-Darboux  formula  [4]  for  the  q-Laguerre  ensemble  is  used 

for  calculating  the  kernel 


N 


1 


KL^N{x,y;  q)  =  \Iwl{x]  q)wL{y;  q)  ^  7^I„(x;  (?)I„(y;  q)  (3.66) 

n=0  '*'" 


where  /c„  =  ,.J~^  „ .  Applying  the  formula  one  obtains 


(i-g)? 


KL,N{x,y;q)        = f — ^ — -j;^^/wL{x]q)wL{y]q)     (3.67) 

Liv+i(y;  q)LN{x;  q)  -  Ls{y]  <?)Iiv+i(a:;  <?) 


X  -y 

-^-^°°     —^—-■lwL{x-q)wL{y-q)  (3.69) 


7 
Looiy/q;  q)Loo{x;  q)  -  Loo(2/;  q)Looix/q;  q) 


(3.70) 


x-y 

The  occurrence  of  the  function  Loo  reflects  the  predictability  of  the  q-deformed  La- 
guerre  polynomials  at  large  N.  The  function  L^o  can  be  approximated  for  (1  —  q)x  '^  1 

by 


L^[x-q)     =     Ji'\2^j{\-q)x-q)  (3.71) 

^     —  cos(-ln(7(l  -  g)x))e^('"'(V('-')")-'^)  (3.72) 

7  7 
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The  potential  wi{x]  q)  is  expanded  by  use  of  Jacobi's  triple  product  formula: 

1/2  oo 

{q,qU{-zq'l'-,q)X-^-q)=     Y.     <?"'^'~^"  (3-73) 

n=  —  oa 

Choosing  z  —  ^^—p^  and  realizing  that  the  first  and  last  product  are  of  order  one,  one 
gets 


'q 

Approximating  the  sum  by  an  integral  the  potential  term  becomes 


wl{x-,  q)  =  {-Y"  exp(--[ln^((l  -  q)x)  +  7/2]^)  (3.75) 

7  47 

After  insfjrting  these  relations  into  the  kernel  one  obtains 

Fhr  \  —  a^~^^      r  IT  it 

KL{x,y;q)     =     ^j^^[{-y/'cosi-ln{{l-q)y))srn{-ln{{l-q)x)) 

-     (^)^/^in(|^ln((l-g)y))cos(|-ln((l-g)x))] 

Applying  the  similarity  transformation 

e    =     ^ln((l-<?)x)  (3.76) 

r]    =     ^ln((l-g)y)  (3.77) 

the  integration  element  w{x)dx  has  to  be  kept  constant  ,  i.e  w{x)dx  =  w{^)d^.   The 

result  is  then 

1         1         e^cos(^r/)sin(^n -e^sin(^ncos(^7?) 
27re«  +  e''  ginh(i^) 

In  the  bulk,  when  ?],  (  -C  1, 77  w  {,  the  exponential  functions  will  vanish  leading  then  to 
the  q-Hermite  kernel  3.61  This  proves  the  assumption  that  in  the  bulk  the  q-Hermite 
and  the  q-Laguerre  ensemble  have  the  same  behavior.  This  result  also  strongly  sug- 
gests that  every  q-deformed  model  will  have  the  same  and  therefore  universal  behavior 
in  the  bulk.  Only  close  to  the  origin  the  ensembles  will  have  a  markedly  different  be- 
havior. As  expected,  this  property  has  an  analogue  for  Wigner-Dyson  ensembles. 
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3.2.3     The  Two- Level  Form  Factor 

Based  on  the  kernel  calculated  before,  physically  relevant  quantities  can  be  cal- 
culated. The  two-level  form  factor,  which  is  simply  the  Fourier  transform  of  the 
two-level  Cluster  function,  is  one  important  example.  It  reads 

b{k)    =     r  e'^("-"'V(x  -  x')d{x  -  x')  (3.79) 

J  —oo 

=     r  e'^^'~''U<\x  -  x')d{x  -  x')  (3.80) 

J  —  oo 

''=(^-"')fl!(.x  -  ,x')  (3.81) 


7^     /■°°    sin(7r(2;  —  x')) 


47r^  J-oo  sinh(|(a;  —  x')) 
This  integral  can  be  calculated  exactly  (compare  with  Appendix  B)  with  the  result 

b{k)     =     —{{\k\  +  27r)coth[(|fc|  +  27r)^/7]  +  {\k\  -  2^)coth[(|fc|  -  2TT)Tr/-fl3.82) 
47r 

-2|fc|coth(|A;i7r/7))  (3.83) 

By  expanding  the  cotangent  hyperbolic  in  the  limit  7  — >•  0,  the  known  Wigner-Dyson 
limit  is  recovered 

m  =  {l      ■"     \l\tl  (3.84) 

At  k=0,  which  is  just  the  integral  over  the  Cluster  function,  the  two  level  form  factor 
takes  the  value 

6(0)  -  coth(— )  -  -^  (3.85) 

Note  the  nonanalytic  dependence  on  7,  which  excludes  the  possibility  to  obtain  the 
above  result  by  a  perturbative  approach  from  the  Wigner-Dyson  limit. 

The  Fourier  transform  of  the  two  point  correlation  function  or  two-level  form  fac- 
tor has  recently  been  of  interest  for  experimentalists  and  theorists  alike  [80,  78,  79] 
The  properties  of  a  spectrum  are  quite  difficult  to  measure  experimentally  due  to  its 
highly  oscillating  nature.  However,  the  Fourier  transform  can  reveal  'used'  or  'unused' 
frequencies  and  should  therefore  represent  a  more  robust  measure  of  the  eigenvalue 
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correlations.  Indeed,  in  the  case  of  random  matrices  where  the  level  repulsion  is  the 
fundamental  phenomenon,  the  Fourier  transformation  should  lead  to  the  occurrence 
of  a  correlation  hole.  The  two-level  form  factor  and  the  occuring  correlation  hole  are 
displayed  in  figure  3.4.  Smaller  frequencies  in  the  spectrum  are  suppressed,  because 
small  energy  differences  are  statistically  disfavored. 

3.3     Statistical  Analysis 
3.3.1     Nearest  Neighbour  Spacing  Distribution 

An  important  quantity  for  transition  models  is  the  nearest  neighbour  spacing  dis- 
tribution. The  spacing  distribution  cannot  be  accessed  with  the  correlation  functions 
alone,  but  additionally  requires  to  solve  an  integral  equation  [4].  This  difficulty  makes 
it  appropriate  to  calculate  the  spacing  distribution  numerically.  Figure  3.5  shows  a 
clear  discrepancy  between  the  Wigner-Dyson  statistics  and  Poisson  statistics  and  is 
fairly  easy  to  evaluate  for  numerical  eigenvalue  distributions.  For  the  transition  model 
it  is  particularly  interesting  to  see  how  the  Wigner  'surmise'  crosses  over  into  the  Pois- 
son neigbour  distribution.  This  family  of  curves  shows  a  remarkable  similarity  with 
numerous  investigations  of  transitional  statistics  from  many  sources  (as  a  few  exam- 
ples [8,  67,  48].  One  interesting  feature  visible  from  the  nearest  neighbour  spacing 
distribution  is  the  still  active  level  repulsion  even  at  finite  7.  Due  to  the  standard 
logarithmic  repulsion  term  the  curves  start  always  at  zero.  However,  the  slope  at 
the  origin  clearly  increases,  which  shows  the  qualitatively  different  nature  of  the  q- 
Hermite  ensemble,  whereas  the  slope  is  constant  for  fairly  general  circumstances  (see 
equation  1.32)  in  standard  random  matrix  theory  in  the  same  symmetry  class. 
The  maximum  of  the  curve  slowly  moves  towards  the  origin,  as  widely  observed 
[8,  81,  36].  This  feature  is  not  reproduced  in  other  purely  phenomenological  or 
semiphenomenological  theories  [15,  16]. 

Of  even  more  speculation  is  the  obvious  'scaling'  point,  where  all  curves  intercept. 
The  long  range  decay  of  the  nearest  neighbour  spacing  distribution  is  in  case  of  a 


53 


_Q 


O 
CM 


LO 


O 


in 
o 


o 
o 


^  =  0    (GUE) 

-  -    -  /?  =  5 

-  —  /?  =  20 

initial    drop    for   finite    n 


0.0 


0.5 


1.0 

k/2*7T 


1.5 


2.0 


Figure  3.4:  Two  level  form  factor  for  yV  =  oo  at  different  7.  The  minima  at  small 
energies  can  be  interpreted  as  a  correlation  hole.  The  figure  was  taken  from  reference 
[52]  where  the  parameter  /3  is  identical  to  the  parameter  7  used  in  here. 
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Figure  3.5:  Nearest  neighbour  spacing  distribution  for  N  =  oo  at  different  7.  The 
variable  t  is  given  by  t  =  x/D,  where  x  is  the  coordinate  and  D  the  mean  level 
spacing.  The  figure  was  taken  from  reference  [77]  where  the  parameter  f3  is  identical 
to  the  parameter  7  used  in  here. 
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Poisson  distribution  apparently  exp(  — x)  in  contrast  to  the  exp(  — x^)  decay  in  the 
Wigner-Dyson  case.  This  large  difference  in  the  decay  behavior  has  prompted  addi- 
tional interest,  because  it  could  lead  to  the  actual  'sharp'  transition  features,  which 
are  so  far  hidden.  One  idea  is  to  define  the  integral  from  the  scaling  point  to  infinity 
(the  tail  of  the  distribution)  as  the  order  parameter  of  the  transition  [53].  Another 
approach  is  to  calculate  the  tail  of  the  distribution  for  the  q-Hermite  model  in  a  more 
analytical  fashion  [82]. 

3.3.2     Number  Variance  and  As-statistics 

In  contrast  to  the  nearest  neighbour  spacing  distribution,  which  manifestly  mea- 
sures the  eigenvalue  correlations  for  close  range  eigenvalues,  the  number  variance 
and  the  As-statistics  are  long  range  measures  of  the  eigenvalue  correlations.  This  is 
achieved  by  integrating  a  correlation  function  over  the  entire  spectrum.  Therefore  it 
includes  also  the  effect  of  eigenvalues  which  are  far  apart.  The  number  variance  VJv 
in  the  interval  (-N/2,N/2)  is  given  by  [4] 
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with  the  two  level  form  factor  b{k)  covered  before.    To  an  excellent  approximation, 
b{k)  contributes  neghgibly  to  the  integral  in  equation  3.86  for  k  >  27r.  We  find  [52] 
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where  Ci  is  the  cosine  integral,  and  p  defined  as 
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goes  to  1(0)  for  7  —^  0(oo).  We  can  obtain  an  explicit  expression  for  the  number 
variance  by  approximating  the  factor  [coth(-^)  —  ~]  under  the  integral  by  kir/S'yN 
up  to  /c  =  3'jNp/Tr  and  by  p  for  larger  k.  The  result 


K(iV)    =    [l-p]N  +  ^ 


sin{3jpN/Tt) 


+4 


'2^2 


•"'i^J^«ra- 


(3.92) 
(3.93) 


agrees  with  the  numerical  evaluation  of  the  integral  up  to  a  small  7  dependent  constant 
for  large  N  as  well  as  corrections  of  order  1/A^.  Equation  3.92  clearly  exhibits  the 
crucial  dependence  of  the  number  variance  on  the  parameter  7.  As  expected,  7  =  0 
corresponds  to  the  exact  GUE  result,  with  a  logarithmic  dependence  on  A'';  increasing 
7  corresponds  to  a  transition  towards  a  Poioson  result,  which  is  linear  in  N .  Note 
that  the  variance  of  an  arbitrary  linear  statistic  /  =  En  f{xn)-,  of  which  K,  is  a  special 
case,  will  depend  on  the  form  factor  3.79  and  is  clearly  no  longer  universal 

Var(/)  =  -/     dk[l-h{k)]\J{k)\'  (3.94) 

where  f{k)  is  the  Fourier  transform  of  f{x).  This  result  generahzes  the  Dyson-Mehta 
theorem[7]  on  the  variance  of  arbitrary  Hnear  statistics  and  reduces  to  it  when  7^0. 
The  As-statistics  is  a  quite  popular  measure  of  long  range  eigenvalue  correlations. 
Already  introduced  in  the  first  chapter,  the  Aa-statistics  is  a  good  measure  of  the 
rigidity  of  a  given  eigenvalue  spectrum.  It  exemplifies  the  large  diiference  between  a 
Poisson  spectrum  and  Wigner-Dyson  spectrum  by  showing  a  linear  and  therefore  large 
deviation  from  a  perfectly  rigid  spectrum  in  the  Poisson  case  and  a  logarithmic  and 
much  more  weaker  deviation  from  the  Wigner-Dyson  case.  This  qualitatively  different 
behavior  motivates  the  calculation  of  the  As-statistics  for  the  transition  model.  The 
q-Hermite  model  cannot  be  evaluated  in  closed  form  but  the  interesting  limiting  cases 
are  analytically  accessible.  The  Aa-statistics  can  be  written  in  the  following  form: 
]S[  I        [^  sin^fTTul 
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For  small  7  (close  to  Wigner-Dyson)  we  can  calculate  the  first  7  dependent  correction, 
by  expanding  the  kernel  i^-^^^S^  *°  ^''^*  °^*^^^  ^^  ^-  ^he  result  is  then  [52] 


„2  AT2 


If  5l       'y  N 

A3(iV,  7)  =  —  [ln(27r,s)  +  7  -  ^ J  +  ^  +  O  [{^N)')  (3.96) 

for  7  <C  1.  This  result  reduces  obviously  to  the  Wigner-Dyson  limit  (7  =  0)  for 
7  ^  U.  On  the  other  hand  for  large  7  the  large  N  limit  is  appropriate,  where  the 
finite  integral  can  be  expanded  around  the  A^  =  00-case. 
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_27r^       exp(i^)-l 

for  7  >  1/yV,  which  reduces  to  the  Poisson  result  for  7  -)►  cx).  Here  C  is  a  numerical 
constant  independent  of  N  and  very  weakly  dependent  on  7.  Note  that  since  the 
ratio  A3/N  is  finite,  A3  always  has  a  linear  dependence  for  large  enough  N,  with  a 
slope  that  increases  with  increasing  7,  approaching  the  Poisson  limit  for  7  —>  00.  For 
values  in  between  a  numerical  evaluation  is  necessary.  For  a  small  set  of  eigenvalues 
(N  small)  the  family  of  curves  shows  still  a  logarithmic  increase  of  the  As-statistics. 
However,  for  larger  N  the  curve  smoothly  crosses  over  into  a  linear  increase.  The 
slope  of  the  large  N,  linear  increase  is  uniquely  dependent  on  7  and  thus  indicates  how 
far  the  statistical  transition  has  progressed.  Figure  3.6  shows  the  complete  solution 
(obtained  from  numerical  evaluation  of  the  integrals  in  equation  3.95)  for  various 
values  of  7.  The  deviation  from  the  GUE  result  as  a  function  of  some  parameter  is 
quahtatively  similar  to  the  deviations  seen  numerically  in  both  the  transfer  matrix  [50] 
and  energy  eigenvalues  [48]  corresponding  to  the  tight  binding  Anderson  Hamiltonian 
for  disordered  conductors,  as  well  as  for  the  eigenvalues  of  the  evolution  operator 
corresponding  to  the  Fermi-acceleration  model  [83]. 
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Figure  3.6:  Aa-statistics  for  N  —  oo  -dX  different  7. 
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Also  note  that  the  slope  exhibits  a  manifestly  nonanalytic  dependence  on  7,  show- 
ing the  nonperturbative  nature  of  the  transition  model.  This  also  indicates  the  subtle 
nature  of  the  behavior  from  7  =  0  to  7  — >■  0  in  the  present  formulation.  A  closer 
analysis  reveals  that  this  behavior  comes  from  the  fact  that  the  results  are  all  eval- 
uated for  an  infinite  set  of  eigenvalues.  Arguments  presented  later  in  this  chapter  as 
well  as  in  the  dynamic  model  formulation  show  that  the  limit  7  — >  0  is  only  properly 
understood  for  a  finite  number  of  eigenvalues.  This  originates  from  the  fact  that  the 
limits  7  ->  0  and  A^  ^  00  do  not  interchange. 

The  behavior  of  the  number  variance  and  the  Aa-statistics  are  apparently  quite 
similar,  even  though  they  contain  slightly  different  information.  The  reason  for  this 
similarity  is  the  analytical  agreement  at  large-N,  where  the  Aa-statistics  retains  only 
the  first  integral  in  3.95  as  a  correction  to  the  linear  term,  i.e. 

1  25  7^    /■«  sin^(7rii) 


Asis  <  1)  =  — s '-  /    du ^ — '—  (3.99) 

^  ^       15         15  47rio        sinh^(7w/2)  ^        ^ 

which  is  precisely,  apart  from  constants,  the  terms  contained  in  the  number  variance 

3.92.    Consequently,  the  large  N  behavior  is  indeed  similar.    However,  the  small  N 

characteristics  are  certainly  different  as  now  all  integral  terms  in  expression  3.95 

contribute  to  the  As-statistics. 

The  small  N  behavior  of  the  As-statistics  is  relevant  for  some  models.  Numerical 
investigations  with  the  tight  binding  Anderson  Hamiltonian  [48]  have  shown  that  for 
small  N  the  A3  curve  can  be  fitted  with  a  power  law  A3(A)  oc  N"" ,  where  the  coeffi- 
cient a  is  a  characteristic  of  the  physical  system  in  question.  But  the  transition  model 
discussed  herein  is  only  able  to  describe  the  gross  features  of  a  transitional  system 
and  is  therefore  not  able  to  recover  the  model  dependent  features  at  small  N.  A  simple 
fitting  procedure  for  the  coefficient  gives  an  alpha  of  0.63  for  the  transitional  model, 
comparing  well  with  the  alpha  of  0.70  for  the  analysis  of  the  eigenvalue  statistics  of 
the  tight  binding  Andersson  model. 

To  better  understand  the  range  of  the  transition,  one  can  plot  the  As-statistics 
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as  a  function  of  7  at  a  fixed  but  sufficiently  large  N.  The  result  is  depicted  in  fig-  ^ 

[■ 
ure  3.7  where  the  ordinate  is  normahzed  by  the  leading  linear  term.     Apparently 

the  transition  is  very  much  hnearly  dependent  on  the  parameter  7.   The  saturation 

at  larger  7  is  not  relevant,  as  the  kernel  3.61  used  in  this  calculation  is  no  longer  valid. 

3.4     Long-Range  Off-Diagonal  Oscillations 

So  far  the  investigations  of  the  kernel  have  concentrated  on  the  bulk  behavior,  ' 

which  is  close  to  the  main  diagonal  (x  ^  y)  and  far  away  from  the  origin,  which  is 
certainly  very  important  for  the  characterization  of  q-models.  The  complete  similar- 
ity in  the  kernels  at  the  bulk  for  both  q-Hermite  and  q-Laguerre  models  (cmp.  eq()), 
supports  this  statement.  However,  close  to  the  origin  the  q-models  have  a  different 
behavior  and  show  therefore  a  system  dependent  behavior.  This  model  dependency 
can  be  used  to  construct  a  specific  theory,  the  most  well  known  case  is  the  transfer 
matrix  formulation  of  conductance  ([21]  or  for  references  see  [85]),  where  the  cor- 
relation function  values  close  to  the  origin  determine  measurable  quantities  like  the ,  t 
conductance. 

The  q-Hermite  model  behaviour  at  the  origin  is  more  difficult  to  connect  with  physi- 
cally measurable  quantities,  because  the  eigenvalue  spectrum  does  not  lead  straight- 
forward  to  a  physical  quantity.  However,  the  inclusion  of  an  external  parameter  into 
the  q-Hermite  model  gives  again  reason  to  raise  the  question  of  the  behavior  close  to 
the  origin  and  is  discussed  in  the  next  chapter.  j- 

In  the  case  of  the  q-Hermite  model  another  phenomenon  occurs  which  is  based  on 
the  reflection  symmetry  towards  the  origin  of  the  underlying  polynomials.  The  two 
point  correlation  function  shows  long  range  oscillations  in  the  off-diagonal  direction 
(i.e.  the  (x,-x)-direction  on  a  x-y-plane).  These  oscillations  do  not  occur,  if  the  system 
is  in  the  Wigner-Dyson  limit,  because  the  correlation  function  depends  in  this  case 
only  on  the  energy  difference  variable.  This  translational  symmetry  is  spontaneously 
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Figure  3.7:  Normalized  Aa-statistics  as  function  of  7  for  N  =  00. 
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broken  in  the  case  that  the  disorder  parameter  7  is  nonzero.  It  has  to  be  emphasized 
that  the  long  range  oscillations  will  only  appear  in  the  q-Hermite  model,  and  are 
therefore  model  dependent. 
With  the  help  of  the  kernel  we  can  characterize  these  oscillations.  Using  the  equation 


'  3.60  with  the  scaled  coordinates  {x  — >■  7X,  y  — )■  71/)  yields 


^    ^  ^cosh(7x)  cosh(7y)      5in(7r(.x  -  y) 
A[x,y,q)       2^    cosh(7(x  +  y)/2)     smh(7(x  -  y)/2)  ^'^^^^' 

Concentrating  on  the  off-diagonal  direction  y  =  —x  the  result  is 

7 
K{x,  —X,  g)  =  -—  coth(7x)  sin(27rx)  (3.101) 

ZTT 

The  cotangent  hyperbolic  factor  will  approach  one  for  large  enough  x.  Then  it  is  easy 
to  see  that  the  long  range  oscillations  have  a  frequency  of  27r  and  grow  linearly  with 
7  as  long  as  7  <<  27r-^.  The  whole  kernel,  based  on  equation  3.47  is  plotted  in  figure 
3.8,  clearly  exhibiting  the  oscillations  also  for  7  ~  1. 

Note  that  the  value  for  7  is  already  relatively  large,  mainly  to  illustrate  the  off- 
diagonal  oscillations.  A  comparison  with  the  analytical  expression  3.101  would  predict 
that  the  height  of  the  peak  is  close  to  {-^Y  ~  0.5  which  is  a  small  underestimation 
of  the  peak  height  from  the  plot  (about  0.4).  This  shows  already  that  the  regime, 
where  formula  3.101  is  valid,  is  starting  to  give  imprecise  results. 

It  is  instructive  to  go  to  even  larger  7  and  plot  the  full  kernel  again.  Figure  3.9 
displays  the  kernel  for  7  =  8.1.  The  off-diagonal  oscillations  are  at  this  value  at 
full  height  and  saturated.  Interesting,  however,  is  the  formation  of  similar  oscillation 
peak  on  the  main  diagonal.  This  phenomenon  is  certainly  not  described  by  the  kernel 
3.61  as  7  >  1,  but  can  be  understood  by  considering  properties  of  the  q-Hermite 
polynomials. 

The  zeroes  of  q  deformed  polynomial  are  much  more  likely  to  be  found  at  one 
position,  compared  to  the  undeformed  case,  where  the  zeroes  cannot  be  fixed  with 
incresing  N.  This  'ergodic'  behavior  of  the  zeroes  in  the  Wigner-Dyson  case  leads  to 
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Figure  3.8:  Full  A''  =  oo-kernel  K^  for  7  =  4.6  for  N  =  00. 


64 


Figure  3.9:  Full  A^  =  cx5-kernel  A'^  for  7  =  8.1  for  N  =  00. 
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a  very  smooth  and  featureless  correlation  functions,  visible  in  translational  invari- 
ance.  However,  for  q  <<  1  the  eigenvalues  are  more  likely  to  stay  at  one  place  with 
I  increasing  N,  as  the  zeroes  of  the  polynomials  give  the  most  likely  position  of  the 

;  eigenvalues.  This  fundamental  difference  between  the  Wigner-Dyson  case  and  the  de- 

j  formed  (  to  be  more  precise  any  q  deformed)  ensemble  is  illustrated  by  the  existence 

of  a  A"  =  cxD  limiting  function  for  the  polynomials,  as  already  used  in  equation  3.71 
ij  for  the  Laguerre  ensemble. 

]  The  oscillations  in  the  (x,-x)-direction  may  have  some  physical  relevance.  Consider 

an  energy  spectrum  of  a  physical  system  in  transition  between  metalhc  and  insulating 
behavior,  which  preferably  exhibits  a  band  structure.  For  comparing  the  results  from 
random  matrix  theory  with  the  energy  eigenvalue  spectrum  one  limits  them  to  be 
from  the  'bulk'  (meaning  far  from  the  band  edges).  As  long  as  the  eigenvalues  are  not 
influenced  by  the  band  edge  behavior,  their  correlations  should  be  uniform.  Therefore 
the  system  in  the  bulk  exhibits  clearly  translational  symmetry.  The  fact  that  in  the 
case  of  the  q-Hermite  model  the  translational  symmetry  around  a  specific  point  (the 
origin)  is  broken  could  indicate  the  potential  relevance  of  the  q-Hermite  model  in 
understanding  the  effect  of  the  band  edge. 

An  important  test  for  this  assumption  is  to  calculate  an  integral  quantity  such  as 
the  number  variance,  which  uses  an  integration  region  including  the  origin.  As  the 
origin  plays  in  the  above  speculation  the  role  of  the  band  edge,  it  is  appropriate  to 
integrate  over  a  region  symmetric  towards  the  main  diagonal,  but  far  away  from  the 
origin  for  small  N.  Only  at  larger  N  the  origin  would  be  included,  and  the  effect  of 
the  long  range  oscillations  becomes  visible.  The  result  of  this  number  variance  calcu- 
lation is  shown  in  figure  3.10.  As  expected,  at  small  N  the  origin  is  not  included,  and 
the  number  variance  resembles  the  results  from  the  known  expression  3.92.  However, 
when  N  is  large,  and  the  integration  includes  the  origin,  the  behavior  of  the  number 
variance  shows  leveling,  constant  number  variance. 
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Figure  3.10:  Number  variance  over  the  full  kernel.  The  upper  line  represents  7  =  4. 
and  the  lower  line  corresponds  to  7  =  1.2. 
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3.5     Finite-N  Analysis 

The  discussion  so  far  has  focussed  on  the  investigation  for  infinite  many  eigenvalues. 
The  results  show  a  qualitatively  different  behavior  of  the  q-Hermite  model,  compared 
to  the  standard  random  matrix  Hermite  model.  However,  the  crossover  from  one 
model  to  the  other  is  difficult  to  apprehend  for  N  —  oo  results.  This  problem  is 
apparent  in  all  the  above  mentioned  results,  which  are  based  on  the  (A^  =  cxo)-kernel 
3.61.  The  kernel  already  has  a  qualitatively  different  behavior,  where  the  power  law 
decay  is  replaced  by  an  exponential  one.  This  difference  is  even  more  pronounced  in 
quantities  which  integrate  over  the  spectrum,  including  eigenvalues  far  apart,  such 
as  the  number  variance  and  the  As-statistics.  Figure  3.7  indicates  the  existence  of  a 
noncontinous  crossover  from  7  <§;  1  to  7  =  0.  No  matter  how  small  7  is,  the  integral 
measures  will  contain  a  linear  term  which  will  dominate  for  sufficiently  large  N.  This 
behavior  prohibits  a  smooth  crossover. 

In  order  to  compare  the  finite  N  results  with  the  N  =  00  calculations  one  has  to 
normalize  the  correlation  functions  correctly.  The  principle  relation  to  use  here  is  of 
course  the  normahzation  of  the  joint  probability  distribution  function  /  P({x„})(i{x„}  = 
1.  For  finite  N  this  leads  to  the  earlier  defined  reduced  probability  distribution  func- 
tions Rn  (in  particular  R2{x,y)  =  a{x)a{y)  —  K^{x^y))  and 


j  R2{x,y)dxdy  =  N{N  -I)  (3.102) 


one  nas 


a{x)dx^N  (3.103) 

I  K^{x,y)dxdy  =  N  (3.104) 

To  scale  the  energy  scale  out  of  the  problem,  i.e.    the  mean  level  spacing  D,  the 
energy  levels  are  rescaled  such  that  the  level  density  is  constant.  In  case  one  chooses 
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the  constant  to  be  unity  the  energy  scale  is  no  longer  present  in  the  problem.  The 
rescaling  is  represented  by  a  change  of  variables  from  {x„}  to  {x^}  and  the  correlation 
functions  in  the  new  variable  have  to  satisfy  the  same  integral  constraints  3.103  as  in 
the  old  variable  set.  For  the  density  this  leads  to  the  following  relation 

aM{x)dx  =  aN{x')dx'  =  [0(x'  +  N/2)  -  0(x'  -  N/2)]dx'  (3.105) 


w 


here  0(x)  =  1  for  z  >  0  and  zero  else.  The  last  equation  comes  from  the  choice  of 
unit  density  in  the  new  variables.  However,  the  choice  of  density  o-(x')  is  not  unique. 
Similarly  the  two  point  correlation  measure  obeys 


R2,N{x,y)dxdy  =  R2,!\[{x'^y')dx'dy  (3.106) 


leading  to 


Using  equation 


Kjf{x,y)dxdy  —  Klf{x',y')dx'dy'  (3.107) 


KU^'^y)   -    4fe^  (3-108) 


dx    dy 

Kjfix^y) 


(3.109) 


a{x)(j{y) 

Inserting  the  formulas  for  the  orthogonal  polynomial  expressions  2.25,  2.26,  one  ar- 
rives at 

For  all  q-models  the  two  point  correlation  function  decays  in  the  bulk  (x,y  >>  1)  as 
a  characteristic  for  the  class  of  q-models  investigated. 

To  compare  the  kernel  for  finite  N  with  the  results  already  obtained  one  can  focus 
on  the  bulk  behavior.  For  this  purpose  we  can  compare  the  sum  3.110  for  increasing 
N  and  compare  with  the  N  —  oo  kernel.  One  should  concentrate  on  the  scaled  poly- 
nomials in  3.24,  by  comparing  the  function  (x'  =  /^  a{x)dx^  Kn{x  +  N/2,  ~x  +  N/2)) 
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and  compare  with  the  kernel  3.61.  The  shift  in  the  kernel  arguments  is  necessary,  we 
would  like  to  observe  explicitly  bulk  behavior.  The  results  have  to  be  scaled  to  result 
in  unit  density. 

This  scaHng  becomes  clear,  once  being  reminded  of  the  fact  that  the  zeroes  of  the 
Hermite  Hn  polynomials  scale  with  y/n.  Dividing  the  argument  by  this  number  results 
in  a  curve  which  can  be  compared  to  the  analytical  kernel.  However,  for  stronger  q 
deformation  the  movement  of  the  zeroes  saturates  at  ;:;.  As  also  used  later,  a  good 
intermediate  approximation  is  ^/ng- 

This  approach  leads  to  a  good  description  of  the  scaled  q-Hermite  polynomials  and 
the  kernel  associated  with  them.  However,  the  procedure  for  the  unsealed  polynomials 
results  in  an  unsatisfactory  description. 


CHAPTER  4 
THE  DYNAMIC  DEFORMED  GAUSSIAN  ENSEMBLE 


The  transitional  model  discussed  in  here  can  be  extended  to  include  an  external 
perturbation.  This  perturbation  is  usually  an  external  field  acting  upon  the  system. 
As  this  dependence  allows  to  predict  new  phenomena  and  to  explain  perturbation- 
dependent  data,  the  physical  relevance  is  immediately  obvious. 

The  parameter-dependent  theory  can  be  treated  in  an  approximate  and  an  exact 
fashion.  The  approximation  provides  easier  access  to  some  important  results  as  well 
as  better  understanding  of  regions  where  the  exact  treatment  is  quite  complicated  to 
evaluate.  In  the  following  both  approaches  are  discussed  in  detail,  and  the  results  are 
put  into  the  context  of  parametric  correlations  for  chaotic  spectra. 

4.1     The  Hydrodynamic  Approximation 
4.1.1     Brownian  Motion  Model 

The  standard  random  matrix  model  describes  a  system  of  eigenvalues  in  static  equi- 
librium. Because  of  the  close  analogy  to  a  2D  electron  gas  the  model  was  also  called 
Coulomb  gas  model.  Dyson  introduced  the  idea  to  extend  the  modeling  efforts  to  in- 
clude a  dynamic  depence  on  an  external  parameter.  Dynamic  means  here  dependence 
of  the  eigenvalues  on  some  'fictitious'  time,  which  doesn't  have  to  be  the  physical  time 
but  any  other  external  parameter  on  which  the  eigenvalues  depend.  Dyson  realized 
that  a  description  in  form  of  a  conservative,  newtonian  formulation  cannot  be  done  in 
a  reasonable  physical  way  [6,  86],  i.e.  defining  velocities  and  inertia  of  the  eigenvalues 
(or  particles  in  the  Coulomb  gas  analogy)  is  impossible.  However,  a  nonconservative, 
dissipative  description  does  not  require  the  definition  of  the  velocity  of  a  single  eigen- 
value.   The  most  straightforward  ansatz  is  thus  to  assume  that  the  eigenvalues  will 
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evolve  according  to  Brownian  motion.  The  underlying  principle  of  the  Brownian  mo- 
tion model  are  the  frictional  forces  each  particle  is  exposed  to.  The  physical  picture 
implies  that  the  eigenvalues  are  experiencing  local  friction  forces  given  by  a  usual 
second  derivative  of  the  distribution  function  and  an  external  'electric'  field  which  is 
the  derivative  of  the  potential,  i.e. 

dW  \ 

£H^»})  =  -^=S[^^;^I-2^<  (4.1) 

The  fluctuating  (random)  forces  give  then  rise  to  a  Brownian  motion.  As  well  known 
from  statistical  mechanics  the  time  evolution  of  a  many-particle  Brownian  motion 
problem  is  given  by  a  Fokker-Planck  type  of  equation.  Mathematically  the  problem 
is  posed  by  assuming  a  set  of  eigenvalues  E^  depending  on  a  variable  r.  The  variable 
T  can  be  seen  as  a  'fictitious'  time,  where  it  is  important  to  use  physical  consideration 
to  identify  the  parameter  r  with  a  measurable  quantity.  At  r  =  0  the  eigenvalues 
have  the  initial  position 

P{{E^].Q)^Y[^{E,-E,,o)  (4.2) 

i 

For  r  >  0  the  system  evolves  according  to  the  following  Fokker-Planck  equation  [86] 

dP      ^    d      _,,^dW        dP  , 

57  =  ?  5^^     (^^  +  ^)  (4.3) 

Here  ^  represents  the  friction  coefficient.  The  term  W{{En])  is  given  by  the  known 
total  potential  of  the  respective  symmetry  class. 

N  N 

W{{Er.])  =  ^Y.^n\E.-  E,\+Y.y{E,)  (4.4) 

We  will  limit  our  discussion  to  the  unitary  ensemble  (/?  =  2),  because  the  transition 
model  is  also  confined  to  this  class.  This  equation  has  the  well  known  equilibrium 
distribution  function  at  i  — >  co 

N  N 

?({£,},  ^)=      n     {E,-E,Yexp{-f3j2V{E,))  (4.5) 
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4.1.2     External  Parameter  Model 

The  crucial  link  between  Dyson's  Brownian  motion  model  and  parametric  cor- 
relations consists  in  the  relation  between  the  'time'  r  and  the  external  parameter. 
Beenakker  [87,  88]  observed  that  the  initial  rate  of  change  for  energy  levels  depending 
on  an  external  parameter  at  small  X  is 

iE^-E.oY  =  X'^  (4.6) 

whereas  the  ensemble  average  of  the  Brownian  motion  model  requires  (again  at  small 
r  from  [86]) 

<  {E^t)  -  E,o{r)r  >=  -^  (4.7) 

This  implies  the  identification  between  the  fictitious  time  r  and  the  perturbation 
parameter  X: 

T  =  X^  (4.8) 

Dyson  converted  the  equation  above  [86]  into  an  evolution  equation  for  the  density 

N 

dE,..dEr.P{{En}, t)J2HE-  E,)  (4.9) 

-oo 


Wl 


th  the  solution  (A'^  =  oo) 

^^^^^'  ^^  ^  M^^^^'  ^^^^^^^)  "  II  dE' p{E\T)ln\E  -  E'\)\  (4.10) 

Linearizing  around  the  equihbrium  (  p{E,t)  —  peq{E)  +  Sp{E,T)  )  gives 

D{E,E')  =  -p-^Pe,ln\E-  E'\  (4.12) 

For  a  constant  density  of  states  p^g  =  po  the  integral  kernel  is  translational  invariant 
D{E,  E')  —  D{E  —  E  )  which  makes  a  Fourier  transform  suitable.   In  Fourier  space 
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this  equation  linearizes  and  results  in 

J^V(fc,  r)  -  -k^D{k)8p{k,  t)  (4.13) 

with  the  apparent  solution 

5p{k,  t)  =  6p{k,  0)  Qx^{-k'^ D{k)T)  (4.14) 

With  the  integral  kernel  as 

D{k)  =  -  r  dEe'^'^^-'p,,  \ii\E-E'\^  ^^^  (4.15) 

J  ~oo  p       I  /C I 


one  gets 


5p{k,  t)  =  5p{k,  0)  exp(-^|/c|r)  (4.16) 


The  decoupling  in  Fourier  space  is  due  to  the  linearization  procedure  and  allows  to 
solve  the  linearized  model  exactly. 

4.1.3     Density-Density  Correlation  Function 

The  principal  interest  is  focussed  on  the  parameter-dependent  density-density  cor- 
relation function 

SiE\E,T,T')=<  J2^KE-F4T))5{E'-E,ir'))>  (4.17) 

The  average  <>  here  is  the  statistical  average  over  the  eigenvalues.  With  the  nat- 
ural physical  assumption  that  the  system  is  translationally  invariant  in  the  external 
quantity,  i.e.  S{E,  E  ,t^t  )  =  S{E^  E  ^t  —  t  )  =■  S{E,  E  ,t).  we  can  write 

S{E.E\t)^<p{E,t)p{E\0)>  (4.18) 

=  <5p{E,T)5p{E\0)>+pl  (4.19) 

The  pair  correlation  function  is  at  r  =  0 

S{E,E\0)  =  pI-K{E,E')  (4.20) 
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Consecutively,  we  get 

S{k,T)    -    S{k,0)exp{-'^\k\T)  +  plS{k)  (4.21) 

=     -K{k)exp{-'^\k\T)  +  pl6{k)  (4.22) 

=    p^<5(fc)  +  (l-6(A;))exp(-^|fc|r)  (4.23) 

We  will  now  use  the  external  parameter  representation  t  =  X'^.   Then  the  velocity- 
velocity  correlation  function  can  be  expressed  as 

C{E,  E\  r,  r')  =<  £^  -^^S{E  -  E.{X))8{E'  -  E,{X'))  >  (4.24) 

leading  to  the  correlator  [  C{k,X)  =  -^■^S{k,X)] 

C{K  r)  =  h\  -  h{k))[2\k\  -  A\k\'X']  exp(-^l^lr)  (4.25) 


So  far,  the  result  is  valid  for  the  standard  random  matrix  ensembles.  In  order  to 
extend  these  considerations  to  the  transition  model,  the  idea  is  now  to  insert  the 
function  b{k)  =  6(/c,  7)  to  obtain  an  approximation  for  the  deformed  case. 

4.1.4     Decay  Behavior  at  Large  X 

A  simple  calculation  gives  for  b{k)  —  1  —  ^ 

5(0,^^«=)  =  ^(  — )^^  (4.26) 


and 


C(0,X->(X))  =  ---^  (4.27) 


Now,  for  6(A-)  =  6(/c,7)  =  ^(|A:|  +  27r)  coth[(|A;|  +  27r)7r/7]  +  (|/c|  -  27r)  coth[(|A;|  - 
27r)7r/7]  -  2|A;|  coth(|A;|7r/7)  one  can  approximate  6(A;,7)  by  the  behaviour  at  small  k 
(i.e.  A;  <<  27r),  which  will  dominate  the  behaviour  at  large  X. 

6(fc)  =  — |A;|coth(^)  (4.28) 

ZTT  7  ' 
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resulting  in  (for  simplicity  fi  =  I] 


5(0.X^oo)  =  ^i,  +  ^ij  (4.29) 


and 


C{0,X)^oc  (4.30) 

This  diverging  correlator  represents  a  problem  which  is  caused  by  the  derivation 
of  the  two  level  correlation  function  b(k).  This  function  was  calculated  (  see  [52] 
)  from  the  associated  kernel,  which  was  first  used  in  [77].  The  derivation  required 
a  A^  — >•  oo  calculation  of  the  sum  over  the  corresponding  orthorgonal  polynomials 
(compare  with  Ismail  [91]).  The  usage  of  the  N  =  oo  kernel  complicates  the  transition 
to  the  chaotic  case.  Despite  the  explicit  dependence  on  the  transitional  parameter 
7,  the  N  dependence  is  vital  for  a  smooth  crossover  to  7  =  0.  As  perceivable  from 
equation  (4.28),  6(/c,7)  is  qualitatively  different  from  b(k,0)  easily  noted  by  the  shape 
of  the  function  close  to  A;  —  0.  The  hmit  N  -^  oo  transforms  a  smooth  crossover  into 
a  sharp  transition.  This  discrepancy  is  probably  due  to  the  assertion  that  the  infinity 
kernel  is  valid  in  the  crossover  regime  only.  Still  for  smaller  E,  the  larger  k  regime 
provides  a  reasonable  approximation  as  visible  from  the  parametric  number  variance. 
We  now  aim  to  approximate  the  effects  for  very  small  7  by  phenomenologically  adding 
a  linear  term  to  b{k,^). 


^i±^lM^  k<k 


Here  is  kt  given  by 


(^-^™-i(i[i.i)coth(M.)  ,>,:         («i) 


'     ^        =     ^^'^°th(-^  4.32) 

Ztv  Ztt  7 


7 
27r2 


(4.33) 


for  7  <<  /?(.  Or 


*'  =  r^  (^■''" 
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Then  the  density-density  correlation  function  results  in 


S{k,  X)  -  p'a5{k)  +  (1  -  6(A;))  exp( ^|A:|X')  (4.35) 

Consequently,  integration  gives  then  for  S{E,X)  =  f  S{k,X)dk  for  X  large 

5(0,X7«1)     -     ^^{:^n^,il-^M-f^k,X'))  +  0{exp{-^k,X'))  J 

-     (l+7/2)5(0,X-^  00,7-0)  (4.36) 

i.e.  for  large  enough  X  the  expected  decay  law  is  recovered  and  provides  also  the  first  ( 

order  changes  for  small  7.  Furthermore,  the  scaling  law  proposed  by  Altshuler  et  al.  f 

[70]  is  recovered.  The  choice  for  the  slope  is  justified  later.  1 

The  range  of  validity  for  the  hydrodynamic  approximation  is  given  by  the  relations 
E  ^  D^X  ^  Xc,  where  the  characteristic  energy  scale  D  is  given  by  the  mean  level 
spacing  and  the  parameter  scale  Xc  represents  the  typical  oscillation  frequency  of  the 
parameter  evolution.  This  is  in  complete  agreement  with  the  picture  of  the  hydro- 
dynamic  approximation  as  a  continuum  approximation,  which  takes  the  correlations 
between  levels  with  large  energy  separation  into  account.  This  was  already  shown  in 
the  last  chapter  in  equation  3.61,  where  the  large  N  bulk  behaviour  was  identical  (and 
also  universal)  for  different  models  in  the  same  universality  class.  It  is  only  natural 
that  the  inclusion  of  an  external  parameter  (which  is  not  intrinsic  to  the  system)  will 
not  fundamentally  change  this  view.  As  long  as  the  eigenvalues  are  far  apart,  the 
hydrodynamic  (or  Coulomb  gas)  approximation  provides  reasonable  results. 

i 

4.1.5     The  Parametric  Number  Variance 

The  explicit  dependence  of  the  eigenvalues  on  the  parameter  makes  a  statistical  t 

measure  desirable,  which  connects  the  regularity  or  chaoticity  of  the  underlying  energy 
spectrum  with  the  external  parameter.  Such  a  quantity  was  introduced  by  Goldberg  I 

et  al.   [89].  As  usual  starting  point  serves  the  staircase  function 

N{E;X)  =  J2Q{E-E,{X))  (4.37) 
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However,  an  unfolding  of  the  spectrum  is  necessary  to  extract  the  general  statistical 
behavior.  Then  the  parametric  number  variance  is  defined 

PNV{X)  =<  [N{E,  Xo  +  X)-  N{E,  Xo  -  X)]^  >  (4.38) 

where  <>  is  the  average  over  intervals  of  length  X  at  fixed  E.  The  limiting  behavior 
is  PNV{X  ^  0)  — >  0.  For  large  X,  however,  one  has 

PNV{E,  X)  ^'^^'^  2ViE)  (4.39) 

where  V{E)  is  the  normal  number  variance.  The  PNV  can  simply  be  calculated  in 
the  hydrodynamic  approximation.  We  can  apparently  write 

PNV{A)  ^2<[N'{E,X)-  N{E,  A  +  '^)N{E,  A  -  ^)]  >  (4.40) 

By  inserting  N  =  J^ p{E',X)dE'  and  <  p{E^,X)p{E2,X)  >=  S{EuE2,0)  ,  taking 
into  account  that  S  is  only  a  function  of  the  energy  and  potential  differences,  at  least 
in  the  bulk  at  constant  density,  one  gets 

<N'{E,X)>=I    1^  dE,dE2<S{EuE2,0)>  (4.41) 

and  further 

<  N{E,  X  +  ^)N{E,X-^)>^  I''  l'^  dE.dE^SiEu  ^2,  A)  (4.42) 

From  4.21  we  know  S 

/OO  n 

dfce"=(^'-^^)(l  -6(A:,(?))e-V!'=i^      (4.43) 
-00 

giving 

f  f  dE,dE2S{E,,E2,A)  =  4 /  dk{^-'°'^^^f-^''^^){l  -  6(fc, q))e-^^'^'(AA4) 
Therefore  we  obtain 

PNV{X)    =    J    dE.J    dE^P  dke'^^^^-^'\l-b{k,q)){l-e-'TWx'^ 

=    4r^Mi^^^)(l-6(fc,,))(l-e-^l^l-^)  (4.45) 
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Note  that  for  X  -^  oo  we  get  PNV{X  ->  oo)  =  2T.^{E)  which  is  the  known  hmit. 

From  4.45  it  is  obvious  that  the  limiting  behavior  for  J'i'  ^-  0  is  PNV{X)  ex  X"^.  This 
clearly  differs  from  the  predicted  behavior  by  Faas  et  al.  [92,  71],  where  a  linear  slope 
was  calculated,  revealing  a  weakness  of  the  hydrodynamic  approximation,  which  is 
problematic  at  small  X  and  small  E  when  the  levels  are  too  close  represent  a  Coulomb 
fluid  anymore.  The  predicted  behavior  is  reasonable  and  agrees  also  with  numerical 
calculations  [71].  Note  that  the  integration  at  small  k  is  no  longer  relevant  here,  as 
the  PNV  is  a  statistic  over  the  entire  spectrum  and  quantifies  integral  behavior.  As 
apparent  from  figure  4.1,  the  actual  behaviour  deviates  from  the  Wigner-Dyson  case 
by  exhibiting  an  increased  value  for  the  parametric  number  variance.  Note  that  this 
deviation  can  be  nicely  explained  by  the  results  from  the  present  transition  model. 

4.2     Exact  Treatment  of  the  Parametric  Correlation  Model 
4.2.1     Ansatz  for  the  Density-Density  Correlation  Function 

As  discussed  before,  since  the  hydrodynamic  solution  has  its  limitations,  a  more 
exact  treatment  of  the  problem  is  useful.  Starting  out  from  the  fundamental  equation 
4.3  the  exact  solution  for  the  undeformed  case  (7  =  0)  was  calculated  by  Beenakker 
and  Rejaei  [88].  The  derivation  uses  a  mapping  from  the  N-level  Fokker-Planck  equa- 
tion to  a  decoupled  N-level  Schrodinger  equation  {p  =  2). 

^-^  -ioT.Jl^-^E  Er'  +  Nm{x^})  (4.46) 

The  single  level  eigenfunctions  of  the  static  N-level  system  can  be  calculated  exactly, 
and  a  natural  solution  can  be  found  by  expanding  the  multibody  dynamic  'wavefunc- 
tion'  in  terms  of  the  single  level  eigenfunctions.  As  a  well  established  procedure  this 
yields 

,V-1      00 

S{E,E\X)=  J2   E  HniE)H4E')Hm{E)Hm{E')exp{u;{n-m)X')        (4.47) 

n—O  m=N 

where  cj  is  the  oscillator  frequency.  The  generalization  can  be  justified  by  noting  that 
the  many-body  wave  function  can  be  written  as  a  Vandermonde  determinante  which 


79 


Figure  4.1:  Parametric  number  variance 


f 
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is  a  product  of  single  particle  wave  functions.  This  implies  that  the  density-density 
correlation  functions  can  be  formulated  in  analogy  to  equation  4.47.  Then  we  arrive 
at  the  following  density-density  function  for  7  >  0. 

N-l     00 

S{E,E\X)=  J2    E  h4E,q)K{E' ,q)hrr.{E,q)hm{E' ,q)exp{{n,  -  m,)uX') 

n=Q  m—N 

(4.48) 

It  follows  from  4.47  by  replacing  the  polynomials  /?,„(x)  and  the  eigenvalues  n  +  |  by  I 

f. 
their  natural  q  extensions  /i„(x,  q)  and  riq  +  ^  [91].  This  expression  shows  the  desired  | 

behavior: 

In  the  limit  X  =  0  the  known  result  from  the  theory  of  orthorgonal  polynomials  is 

recovered,  i.e. 

7V-1      CO  I 

S{E,E')     =     R2{E,E')^^    J2  hn{E,q)K{E',q)h^{E,q)h^{E',q)w{E)w{E')  I 

n—0  m=N 

N  N 


'£hl{E,q)w{E)5{E-E')-{^h4E,q)K{E',q)^w{E)w{E')Y 
a{E)5{E-E')-K'{E,E')  (4.49) 


Here  R2{E,E')  represents  the  two  level  (static)  correlation  function,  a{E)  the  level 
density  and  K{E,E  )  the  kernel.  The  eigenvalue  n,  =  yz-  appears  in  the  formal 
single  level  q-Greensfunction 

'diss 

Gs{E,  E\X)  =  Y.  ^n{E,  q)K{E\q)  eM-un.X'')  (4.50) 

which  is  the  solution  of  the  one  level  eigenvalue  equation.  Therefore  these  two  hmits 
for  the  density-density  correlation  function  are  exact. 

This  illustrates  that  equation  4.48  denotes  a  possible  q-generalization  of  the  density- 
density  correlation  function  in  4.47. 

This  q-generalization  will  not  be  unique,  though,  as  it  is  well  known  that  many  pos- 
sible q-Hermite  polynomials  exist.  A  more  fundamental  derivation  of  equation  4.48 
based  on  a  generalized  Fokker-Planck  equation  has  eluded  us  so  far. 


4.2.2     Large  X  Behavior  of  S 

In  order  to  understand  the  dynamic  properties  of  the  transition  model,  the  expres- 
sion 4.48  has  to  be  analyzed  for  large  X.  Starting  from  the  density-density  correlation 
function 

N-l     CO 
m=0  71=^ 

(4.51) 

we  can  choose  E  =  E    =0.    Then  one  can  calculate  the  values  for  the  q-Hermite 
polynomials  by  simply  applying  h2n{x.q)  =  2xh2n-i[x,q)  -  ^~fZi'  h2n-2{x,q). 

/^2n(0)  =   ^^^'^       J         '^  (4.52) 

The  Hermite  polynomials  with  odd  numbers  vanish  at  zero  argument.  Including  the 
normalization  one  can  write  then 

/^2n(0)/.2n(0)  =       „    ^^     rri  f  '^n  (4-53) 

q-^l{q,q)ooVulo{l-q'+^)  ^         ' 

The  potential  can  be  written  as 

{~qeyi^{2x):q)^{-qex^{-2x);q)^ 


V  TT  64(7rii,/7;exp(-7r'^/7)) 


(4.55) 


which  gives  for  E=0 

2  I 

^(0)  ^  •^'?^'^'(9;g)ooV7  ,    ,.„ -, ^^7-  (4.56) 

Vtt  1 +  2E„=i(-l)"exp(-n27r2/7) 

It  follows 

sixn)  =  aztV      [nr(i-^/^'-^)P  [nr(i-g^'-')]^ 

^tTvi'o  g-"7(«?,  9)00  n?:o(i  -  g'+M  9-7(9,9)00  n?ro(i  -  g'+O 

exp(cj(m,  -  nq)X'^)  (4.57) 

Here  the  constant  c  will  be  chosen  so  that  the  density  is  given  by  a  constant  inde- 
pendent of  7  and  N.  The  frequency  u:  is  given  by  cd  =  cj ji  (see  appendix  B)  and  the 
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potential  term  A  is 


^_2g2^ 1 


/OO 


The  constant  k  {A  ^  2/7r),  the  first  sum  and  the  second  sum  can  be  evaluated  for 
small  7,  i.e.  ^-^N  <<  1.  The  term  inside  the  sums  can  be  estimated  for  small  7 

„   _  [nr(i-^^-)i' 


3\        n  ^2n-r 


:i-g)(l-a--(l-9 
(l-g2)(i-54)..(i_g2„) 


Wl 


ith 


r^^ 


2n 


and 


The  logarithm  of  g  can  be  evaluated 


resulting  in 


_n-Y<< 


1  exp(-4^) 


722 


(4.60) 


(4.61; 


1  -  g'"-^  _  e-^(^"-^)(l-e-^) 

1   —  g2n  ~  "i^    _  g-2n7 

=.«!  i_  ^(1  +  ^/2)-^^  (4.62) 

"-^«^  (1+7)(1-^^)  (4.63) 


9n    =     (l+7)"n(l-^-^)  (4.64) 


ln(^„)     =    nln(l+7)  +  Eln(l-i-^)  (4.65) 

~    -^-ln(n2(i-2))  (4.66) 


^"  -^^^^  -fc;^  (4.67) 


By  neglecting  terms  of  smaller  order  in  717,  one  arrives  finally  at 

S{X)  =  —  ^  g"^^e-™'^^  ^  <?"^„e-"'^^  (4.68) 

^7  m=0  n=jV 
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We  are  able  to  determine  the  constant  c  by  demanding  a  constant  density  a.  Since  we 
can  calculate  the  density  only  at  E  =  0  exactly,  we  will  choose  c  so  that  the  density 
is  unity  at  that  point.  This  assumption  is  clearly  justified  by  looking  at  the  figures 
3.1  and  3.2  ,  which  illustrate  that  the  density  does  not  significantly  change  close  to 
the  origin. 
A  short  calculation  gives  (after  approximating  Qk) 


n 

The  condition  (Jn[^)  =  const.  =  1  leads  to 


'^^(o)   =   ^E-T^r  (4.69) 

v/2c7r 


{N\  (4.70) 


77r^ 


u 


By  replacing  the  sums  by  integrals 


(4.71) 
(4.72) 


^  g'^^^e-'"'^'  ~  /     A^e"^'-^'  (4.73) 


0     ,/mn 


m=0  ■'^        V    "■? 

and  substituting  rriq  ^  m  we  arrive  at 


^(^)  =  W  L      ^^^^  L    -^T-^"  (4-74) 


/Vg  Jo        -Jm         Jnc      \/n 


After  a  further  substitution  step  (5^  =  m)  and  insertion  of  the  relations  for  c  and  lo 
one  obtains 

S{X)  =  y^dsl^^dte^^'-^'^^'/-  (4.75) 

This  equation  represents  the  generalization  to  the  standard  expressions  in  [88]  and 
[70]  for  disordered  systems. 
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The  equation  4.75  can  be  used  to  determine  the  large  X  behaviour.     Starting  by 
calculating  the  two  integrals  as 


/ 

Jo 


ds     —        ,  I  -^ds'  = 


=e 


N^ujX 


V^JJO  Jo            V^  ^N,uX^ 
^^     —     /             ^+'  _       V  '   „-wX  h   I  _ 

IqL 

where  the  large  argument  expansion  for  the  incomplete  Gamma  function  was  used  ( 
7(a,x)  =  4°°  e~H°-'~^dt  ~  x''~^e~^).  Multiplying  the  two  integrals  as  well  as  inserting 
c  and  id  results  in 

S{X)  =  ^^(1  -  v'l-exp(-7/V)e-^^-P(--^)-=/^-)  (4.76) 

This  expression  has  for  X'^  »  /ie"^^  the  known  limit  S{X)  —  ^,  which  was  already 
obtained  from  the  hydrodynamic  approximation  4.29.  However,  for  X"^  <^  y.t'^  the 
correlation  function  decays  with  the  markedly  different  behavior 

^W  =  ^e-^  (4.77) 

This  result  demonstrates  the  breakdown  of  the  universal  decay  law  for  for  sufficiently 

large,  but  not  too  large,  X  (/i  <C  A'^  <  e~^^).    Equation  4.76  contains  not  only 

the  expected  general  decay  term  a   ^,  but  also  an  exponential  correction.    This 

additional  term  is  quite  small  at  large  X  and  small  -^N ^  but  for  larger  7A^  it  will  lead 

to  weaker  decay  behavior  at  large  X  compared  to  the  universal  result. 

For  small  X,  which  is  defined  by  the  relation  X^  <<  /i,  the  density-density  correlation 

reads 


^(^'  =  ^"^'^ 


in  terms  of  the  error  function  erf{x)  =  -^  J^  due~'^' .  By  distinguishing  the  two  cases 
X^TT^  «  2fj.jNg  and  X^n^  >  2/j.'yNq,X'^iT'^  <<  2/i  one  can  easily  understand  the 
small  X  behavior.  In  the  case  X'^tt'^  <<  2fi-yNq  the  expansion  of  the  error  functions 
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leads  to 


S{X)  =  ^[~-l]  +  0{X)  (4.79) 


On  the  other  hand  we  have  for  X  tt    >  2n^Ng,  X  tt    <<  2/i 


This  result  points  out  the  special  nature  of  the  value  X=0.  For  small,  but  not  too 
small  values  of  X  the  density  correlator  has  a  large  term  of  the  order  of  -^.  At 
very  small  X,  however,  the  density  correlator  converges  to  the  value  i —  —  1.  This 
clearly  describes  a  peak  at  the  origin.  In  case  of  the  undeformed  model  we  have 
-^Nq  =  1  —  e"^^  —  0  and  therefore  a  divergent  value  at  S{X  —^0).  A  variation  of  the 
above  calculation  obtains  indeed 


indicating  the  diverging  result.  The  calculation  demonstrates  that  the  divergence  in 
the  case  of  vanishing  7  broadens  into  a  wider  and  finite  peak  for  increasing  disorder. 
The  interpretation  of  this  phenomenon  is  connected  to  the  presence  of  level  self- 
correlations.  The  level  repulsion  excludes  levels  at  the  same  energy,  but  the  statistical 
fluctuation  causes  interference  of  the  same  level  at  the  same  energy.  This  contribution 
appears  naturally  in  the  parameter-independent  density-density  correlation  function 
S{E,  E')  a,sa  function  (equation  4.21  and  also  4.57).  The  finite  distance  between  the 
levels  in  the  transition  model  prohibits  on  the  other  hand  the  strong  self-correlations, 
so  that  the  correlation  function  converges  to  a  large,  but  finite  value  (Here  we  assume 
jN  «  1). 

The  phenomenon  of  the  broadening  of  a  correlation  function  at  small  arguments  ap- 
pears also  in  the  velocity- velocity  correlation  function  which  we  will  discuss  next. 
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4.2.3     The  Velocity  Correlator 

In  the  Fokker-Planck  formulation  of  equation  4.3  the  density-density  correlation 
function  is  the  fundamental  quantity.  However,  from  physical  consideration  the  ve- 
locity correlator  is  of  equal  importance. 

The  velocity  correlator  can  directly  be  derived  from  the  density-density  correlation 
function  as 

CiE,  E\  X,  X')  =  ^^  jdEJ  dE'S{E,  E',  X,  X')  (4.82) 

The  evaluation  of  the  complete  function  C{E,  E\  X,  X')  is  technically  difficult,  since 
a  large  N  approximation  for  the  q-Hermite  polynomials,  which  partially  constitute  the 
eigenfunctions,  for  all  arguments  is  presently  not  available.  However,  in  analogy  to  the 
density-density  correlation  function,  it  is  possible  to  calculate  the  velocity  correlator 
at  vanishing  X.  The  technical  procedure  is  outHned  in  appendix  B  with  the  result 

[/V/2] 


27^£'o„.[^,]+i(2n  +  l),!(2m),!((2n-l),-(2m),)2     ^^'^'^ 

(4.84) 


/92 

Mm)^-[2n+l)q)^(X-X  )2 


dXdX' 

t'^f-^      ^         g^"+^-((2n-l),!.TO2m  +  l),!!)^ 


^i    n=[m+i  (2"),!(2m  +  l),!((2n  -  3),  -  (2m; 


'd' 


9  ,{{2m),-(2n-3),MX-X')^^ 


for  N  >  2.   Since  the  sum  representation  is  difficult  to  evaluate,  the  approximation 
by  an  integral  is  desired.  As  outlined  in  appendix  B,  we  get 

yLi    */  0  i/  1  S     ~~  Z 


7T 


{x-x'Y  /■! 


2^ 


ds  /  ^^  dt{e  +  3^)e(='-*')(^-^  )^-V2,.    (4^ 


From  here  on  we  use  X  =  X  ~  X' .  This  formula  can  be  evaluated  for  small  and  large 
X.  For  large  X  we  write  the  first  integral  as 

C{X)     =     -J  dC-^)j^dsj^^^^dt{s'  +  f)e~^^'-'')^'-'l'^    {Am) 

(4.90) 
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In  both  integrals  appears  now  the  same  expression.    Concentrating  only  on  terms 
which  are  relevant  for  4p  >>  1  we  obtain  as  an  intermediate  result  [ 


/o 


dsl'^^dt{f  +  s')e^^'-''^^'-y''^  (4.9r 


,2 


-~^[l-iy/W,  +  ^=)e~^^'^'^-''-^    '^  +  0{1/X')  (4.92) 


Inserting  this  result  in  equation  4.87  and  evaluating  the  parameter  integration  in 
leads  then  to 

1  (/^ 


^W--;^i^  +  — ^:^ — e    ^'^  (--.-^0(1/;^^)  (4.93) 

This  result  recovers  the  universal  velocity  correlator  in  the  accordingly  called  'uni- 
versal' regime,  i.e.  CiX)  =  -^.  Our  result  also  conveys  that  even  at  7A^  ^  0  we 
reestablish  the  universal  correlator  for  sufficiently  large  X,  in  particular  at  '^"^^  >>  1. 
It  is  interesting  to  note  here  that  this  universal  decay  behavior  can  be  also  recovered 
from  the  completely  chaotic  Sinai  billiard  by  use  of  the  Gutzwiller  trace  formula  [92]. 
Furthermore,  it  is  obvious  that  the  correction  term  in  jN  leads  to  a  smaller  velocity 
correlator.  This  is  in  perfect  agreement  with  the  result  for  the  density-density  corre- 
lation function  4.75  as  well  as  with  the  intuitive  picture  that  the  deformation  in  the 
transition  model  weakens  the  correlations.  Even  though  this  is  directly  apparent  only 
from  the  energy  level  point  of  view,  the  parametric  correlations  will  show  the  same 
tendency  as  they  are  tightly  coupled  to  the  energy  eigenlevel  behaviour.  In  particular, 
for  the  Poisson  case  one  would  expect  a  velocity  correlator  identically  zero. 
However,  the  weak  exponential  correction  term  is  unlikely  to  be  discovered  in  exper- 
imental or  numerical  work. 

For  small  X  the  velocity  correlator  has  a  nonuniversal  behavior.  As  the  transitional 
model  possesses  a  new  scale  7A^,  the  correlator  will  exhibit  new  features  at  small  X. 
In  writing  the  expression  4.75  as  (take  r]  =  —) 
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-2r^X's'e^(^'-''^'''  +  2rjXh'e^^^'~''^^']  (4.95) 

the  five  individual  terms  can  be  analyzed  separately.  Neglecting  for  •qX'^  <<  1  the  s 
integral,  the  first  and  last  terms  are 

i 

In  the  approximation  chosen  the  first  and  the  third  term  are  similar  and  combine 
with  the  fifth  term  to 

1        -^         -.X- 


,  e    ^-,  -e---  )  +  0{^X^ri)  (4.96) 

The  fourth  term  can  be  shown  as  small  in  0[y/Xh]),  but  the  third  term  obviously 
contributes  a  constant  term. 


"=l  +  ^[i(l-7iV,)ln(i^,-^|  (4.97) 


The  result  can  be  written  as 


nX" 


^W  =  >^--V^i^(-.^,.Ki±^)„         ,«3, 


The  entire  correlator  is  depicted  in  figure  4.2.  calculated  from  the  straigthforward 
integration  of  the  integrals  4.75. 

It  is  apparent  from  this  expression  that  the  velocity  correlator  exhibits  a  peak 
around  X=0.  The  width  of  the  peak  is  of  the  order  7A^,  which  approaches  zero  in  the 
weakly  disordered  case.  In  this  limit  the  peak  diverges  and  becomes  a  (^-function.  This 
can  be  seen  in  the  limit  of  7/V  ->  0,  where  the  only  term  remaining  is  the  nonanalytic 
term  C{X)  =  7;^^  '"'''  +  0{const.).  This  term  describes  the  divergence  as  well  as 
its  width.  The  peak  is  visible  in  figure  4.3  for  a  small  7yV. 

At  7A/'  =  0  the  limit  X  ^  0  is  not  well  defined,  as  it  is  clear  from  equation  4.75 
that  the  Hmiting  procedure  C{X  -^  0,77V  =  0)  will  lead  to  C{X  =  0,7A^  =  0)  =  0. 
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Figure  4.2:  Unsealed  Velocity  Correlator  for  five  different  values  of  7 
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Figure  4.3:  Peak  structure  in  level  sensitivity 
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This  is  an  explicit  violation  (quite  apparent  from  figure  4.3  without  peak)  of  the  sum 
rule 

y^     C{X)dX  =  0  (4.99) 

which  is  an  obvious  conclusion  from  equation  4.82.  The  ^-function  peak  at  X=0  has 
to  be  accounted  for  the  Wigner-Dyson  limit. 

The  problem  to  evaluate  the  velocity  correlator  for  q=l  at  X=0  was  already  rec- 
ognized by  Simons  and  Altshuler  [70].  They  suggested  a  level-broadening  procedure, 
where  the  statistical  averaging  is  not  performed  over  sharp  ((^-function)  energy  levels, 
but  over  slightly  smeared  levels  (a  gaussian  broadening  function  was  used  in  [70]). 
This  broadening  leads  to  a  well  defined  as  well  as  positive  C(0)  and  subsequently  to 
compHance  with  the  sum  rule.  This  procedure  is  physically  justified,  since  the  energy 
levels  show  nonzero  equilibrium  fluctuations. 

Note  that  the  effect  of  the  broadening  procedure  on  the  velocity  correlator  is  quite 
similar  to  the  inclusion  of  disorder  in  the  transition  model.  A  singularity  at  X=0 
is  smoothened  out  into  a  finite  peak  as  a  function  of  a  parameter  which  determines 
the  width  of  the  peak.  It  would  certainly  be  quite  interesting  to  consider  the  level 
broadening  procedure  also  for  the  transitional  model  in  order  to  include  the  universal 
(weakly  disordered)  regime.  However,  this  would  pose  strong  technical  problems,  as 
the  'smoothened'  ensemble  average  over  a  q-Hermite  ensemble  is  not  clear.  There- 
fore, we  will  restrict  here  the  physical  validity  to  7yV  >  1.  In  this  regime  the  peak 
is  determined  by  the  disorder  parameter  -/VV,  because  the  characteristic  energy  scale, 
the  mean  level  spacing,  is  set  to  one. 

4.2.4     Scaling  of  the  Velocity  Correlator 

The  universal  behavior  of  random  matrix  results  at  7  =  0  is  assumed  to  appear 
also  in  results  with  parametric  correlations.  However,  the  additional  variable  (the 
external  parameter  X)  will  result  in  the  occurrence  of  a  new  characteristic  scale  in  X, 
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which  can  be  thought  to  be  analogous  to  the  mean  level  spacing.  For  the  parametric 
correlations  the  velocity  correlator  at  zero  field  is  a  suitable  description  (see  [68]). 
Moreover,  after  rescaling  with  the  characteristic  quantity,  the  correlator  becomes 
independent  of  the  physical  properties  of  the  involved  system  and  displays  universal 
behavior  (see  [68,  69]).  This  was  even  claimed  for  weakly  disordered  or  not  completely 
chaotic  systems,  so  that  a  comparison  to  the  presented  theory  is  appropriate.  We  will 
limit  our  discussion  here  to  disordered  systems  as  the  theory  is  more  developed  on 
this.  For  this  investigation  one  can  first  observe  that  the  rescaling 


X    -^    x^^C{0)X  (4.100) 

c(x)     =     C{x)/C{0)  (4.101) 

is  performed  by  a  constant  appropriately  chosen  by  the  zero  argument  correlator. 
The  significance  of  the  zero  velocity  correlator  resides  also  in  the  fact  that  it  can  be 
identified  with  the  conductance  in  the  case  7  =  0.  The  assumption  consists  in  the 
applicability  of  a  fluctuation-dissipation  theorem  for  the  random  matrix  ensembles. 
One  indicator  is  already  the  presence  of  a  Fokker-Planck  equation  which  is  a  funda- 
mentally nonconservative  equation.  The  occurrence  of  a  dissipative  system  indicates 
already  that  usual  transport  relations  can  be  applied.  The  validity  of  the  theorem 
requires  that  the  external  perturbation  remains  small,  which  will  only  affect  the  par- 
ticles in  the  medium  considered  but  not  the  medium  itself.  The  system  reacts  to 
the  external  perturbation  in  sole  dependence  on  the  equilibrium  fluctuations.  This 
necessitates  to  formulate  the  dissipation-fluctuation  theorem  in  the  statistical  context 
used  so  far.  The  averaged  energy  is  e{X)  =<  p{X)  >  where  the  brackets  indicate  a 
statistical  average  over  the  distribution  P{{En},X),  and  p  is  given  by  equation  4.9. 
One  can  then  assume  ^^  =  0  without  loss  of  generality.  The  fluctuation-dissipation 
theorem  can  be  written  as  [69] 

de        ,  de  ,,,5X,, 
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where  (J|)^  ^  lim^,,^^  <  E,  |f  ff  ^(i?  -  E,{XME  -  E,[X'))  >^^^,^„,  and 
q;  is  a  constant.  A  similar  relation  was  originally  proposed  by  Wilkinson  [93,  94]. 
The  derivative  in  the  external  parameter  of  the  energy  corresponds  to  the  velocity 
correlator  C(X),  and  the  equation  above  can  be  written  as 


-  =  .hCm^?  (4.103) 

This  allows  to  identify  the  velocity  correlator  at  X  =  0  with  the  conductance  de- 
scribing the  sensitivity  of  the  spectrum  to  changes  in  the  external  parameter  X.  This 
identification  was  in  another  context  also  derived  in  reference  [95]  where  the  differ- 
ent nature  between  the  dissipative  conductance  given  by  the  Kubo  formula  and  the 
conductance  due  to  curvature  of  the  level  at  the  origin  C(0)  was  pointed  out.  The 
equality  between  those  two  quantities  is  only  given  for  system  sufficiently  close  to  the 
the  universal  (metallic)  regime. 

Even  though  the  Thouless  argument  is  questionable  at  larger  disorder,  we  can  still 
demonstrate  the  breakdown  of  the  universality  by  scaling  the  velocity  correlator  by 
the  value  C(0).  The  result  is  depicted  in  figure  4.4. 

From  the  expression  4.98  and  4.93  the  scaling  procedure  yields  for  the  transition 
model  at  small  arguments,  i.e.  -^^  «  C(0) 


c{x)     =     I  \A^(l  +  7^.-2(7yV,)3/^)  ^^ 

(1  -  7iVJ/n[(l  +  y7^v7)/(l  -  ^/^)]^ 


i.104) 


whereas  for  large  argument  {~^jf  >>  C(0)) 

1       \liK+ 


I  cix)     — 1 ^ — -e    '"[(i+v^N^J/li-vTA^)]' 


(4.105) 


Note  that  for  all  7/V  the  universality  still  holds  true  for  large  enough  x,  where  the 
correction  term  becomes  exponentially  small.  For  smaller  x,  however,  the  strong 
disorder  parameter  can  not  be  scaled  away,  and  the  universality  in  this  regime  breaks 
down.  This  breakdown  is  associated  with  the  breakdown  of  the  universahty  of  spectral 
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Figure  4.4:  Scaled  Velocity  Correlator  for  three  different  values  of 


95 

correlations  at  strong  disorder  [52].  The  numerical  results  from  figure  4.4  support  that 
conclusion.  The  scaled  velocity  correlator  for  the  largest  value  of  7  =  2.4  is  clearly  not 
fitting  the  known  asymtotics  anymore.  In  particular,  for  small  X,  the  7  =  2.4  curve 
hes  below  the  asymptotic  behavior  even  though  the  third  correction  in  the  asymptotic 
I  curve  is  positive. 

These  predictions  should  be  possible  to  verify,  for  example  by  comparing  the  behav- 
ior with  numerical  results  for  the  tight  binding  Anderson  model  for  strong  disorder 
near  the  critical  regime. 


CHAPTER  5 
DISCUSSION 


This  chapter  recapitulates  the  basic  transitional  statistical  model  at  the  core  of  this 
work  and  restates  the  central  results  of  the  approach.  These  results  are  then  compared 
with  actual  expectations  from  physical  systems,  and  several  predictions  are  drawn. 
The  beneficial  contributions  of  the  model  will  become  clear  in  a  detailed  comparison 
with  other  approaches  to  the  transitional  statistics  problem.  Finally  an  outlook  will 
be  presented  focussing  on  the  future  of  this  physical  area  and  on  additional  possible 
model  developments. 

5.1     The  Transitional  Statistics  Model 

Complex  quantum  systems  exhibit  certain  common  characteristics.  The  analysis 
of  their  energy  levels  often  reveals  a  similar  universal  behavior. 

Standard  random  matrix  theory  has  been  instrumental  in  explaining  general  level 
statistic  behavior  in  chaotic  and  metaUic  systems.  It  is  based  on  the  underlying 
assumption  that  the  energy  levels  of  a  complex  system,  about  which  we  have  limited 
knowledge,  are  better  described  by  a  Hamiltonian  matrix,  which  is  as  random  as 
possible  and  subject  to  certain  global  symmetry  conditions  hke  hermiticity. 

The  great  success  of  random  matrix  theory  for  otherwise  difficult  to  approach  sys- 
tems and  the  possibility  to  connect  phenomena  from  seemingly  different  areas  in 
physics  has  made  it  a  widely  accepted  and  helpful  interdiscipfinary  theoretical  tool  in 
physics. 

Many  of  the  areas  implied  contain  also  physical  systems,  which  are  neither  absolutely 
chaotic  nor  completely  ordered  and,  therefore,  show  a  different  level  behavior.    The 
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characteristics  of  these  systems  are  in  between  the  random  matrix  theory  and  a  theory 
of  complete  level  randomness  (Poisson  statistics).  Of  most  interest  are  systems,  whose 
characteristics  change  as  function  of  a  certain  intrinsic  parameter,  which  determines 
its  degree  of  chaos  or  disorder.  By  varying  this  parameter  the  system  crosses  over 
from  a  chaotic  spectrum  to  a  regular  (integrable)  spectrum. 

To  understand  these  transitional  systems  better,  it  is  a  well  motivated  idea  to 
develop  a  new  extended  random  matrix  theory,  which  can  account  for  the  change  in 
behavior.  The  model  described  in  this  work  represents  a  contribution  in  this  direction. 
Its  goal  is  not  only  to  explain  already  observed  spectra,  but  also  to  predict  general 
features  for  systems  in  spectral  transition. 

The  model  can  intuitively  most  easily  be  understood  in  terms  of  the  two-dimensional 
Coulomb  gas  model,  which  replaces  (in  an  exact  mapping)  the  rather  abstract  energy 
levels  by  two-dimensional  particles  (electrons).  The  particles  repel  each  other,  but 
they  do  not  escape  from  each  other  due  to  the  presence  of  a  confining  potential. 
The  repelling  force  has  been  identified  as  the  source  of  the  most  fundamental  effect, 
the  level  repulsion.  As  the  system  changes,  the  level  repulsion  is  still  present,  albeit 
significantly  less  pronounced. 

Therefore  it  is  natural  to  keep  the  two  particle  repelling  potential  and  to  generahze 
the  confining  potential.  Standard  random  matrix  theory  features  a  very  strongly 
confining  potential,  which  keeps  the  eigenvalues  close  together.  This  potential  has 
the  characteristic  behavior  to  confine  the  particles  to  a  finite  distance  from  the  origin, 
even  when  the  number  of  particles  diverges  in  the  thermodynamic  limit.  This  pro- 
perty leads  obviously  to  a  vanishing  distance  between  particles  for  iV  ->  oo,  which 
explains  the  picture  of  the  particle  system  as  compressible  fluid. 

The  system  does  not  exhibit  compressibihty  in  case  of  a  weakly  or  softly  confining 
potential,  as  the  particles  in  the  thermodynamic  limit  move  towards  infinity  instead 
of  being  compressed.   The  system  is  therefore  incompressible,  and  one  can  expect  a 
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qualitatively  different  behavior. 

The  justification  for  changing  the  confining  potential  can  be  found  within  the  ma- 
ximum entropy  approach  framework.  It  states  that  the  potential  is  determined  by 
using  the  physical  eigenvalue  density  as  input.  As  the  density  in  both  random  matrix 
ensembles  are  constant,  the  transitional  model  complies  with  the  general  approach. 
The  model  can  be  exactly  solved  in  terms  of  orthorgonal  polynomials,  which  certainly 
represents  a  strength  of  this  formulation.  It  also  introduces  one  parameter  governing 
the  transition  from  a  chaotic  to  a  regular  spectrum. 

The  results  of  any  transitional  model  have  to  be  carefully  evaluated.  The  problem 
lies  in  the  fact  that  directly  observable  and  predicted  quantities  such  as  the  unsealed 
correlation  functions  are  difficult  to  compare.  Therefore,  the  usual  way  of  comparison 
is  to  establish  a  statistical  analysis,  in  which  the  quantities  are  rescaled  so  that  the 
level  density  becomes  constant.  Consequently,  the  fluctuations  gain  universal  meaning 
and  can  be  interpreted.  These  statistical  results  can  then  be  compared  with  either 
numerical  or  experimental  data.  Any  transitional  model  should  be  able  to  recover 
the  two  fundamental  statistical  limits,  the  Poisson  statistics  on  the  one  hand  and  the 
Wigner-Dyson  statistics  on  the  other  hand.  Furthermore  a  model  should  be  capable 
of  describing  the  crossover  between  these  two  limits  as  a  function  of  the  transition 
parameter. 

The  transitional  model  evaluated  herein  possesses  all  those  features.  In  particular, 
the  two-level  correlation  function  R2  ^  I  -  A'^  has  an  explicit  crossover  into  the 
well  known  Wigner-Dyson  correlation  function  [4]  for  the  parameter  limit  7-^0. 
The  other  limit  is  obtained  by  imposing  7  ^  cx).  The  two-level  correlation  function 
becomes  then  i?2  =  1,  as  expected  from  Poisson  statistics.  However,  this  limit  is 
discontinuous,  as  it  can  be  seen  from  the  fact  that  for  all  7  <  00  we  have  /?2(0)  =  0, 
while  for  7  =  00  the  value  is  /?2(0)  =  1.  This  is  a  remnant  of  the  fact,  that  the 
transitional  model  contains  still  the  two  body  repulsive  force,  i.e.     two  levels  are 
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disfavored  to  be  at  the  same  energy.  In  addition,  it  is  quite  interesting  to  note  the 
fundamentally  different  decay  behavior  of  i?2,  indicated  by  the  change  from  a  power 
law  to  exponential. 

The  nearest  neighbour  spacing  distribution  is  a  measure  which  shows  a  distinct 
difference  between  Wigner-Dyson  and  Poisson  statistics.  The  result  of  the  model  in 
figure  3.5  interpolates  smoothly  between  these  two  extreme  cases.  The  maximum 
present  in  the  Wigner-Dyson  spacing  moves  for  larger  transition  parameter  towards 
the  origin  and  joins  in  the  Hmit  7  -)►  cxd  into  the  exponential  decay  of  a  Poisson 
spacing.  Two  important  observations  are:  First,  the  slope  of  the  nearest  neighbour 
spacing  distribution  at  the  origin  can  be  shown  to  directly  depend  on  the  symmetry 
class  for  standard  random  matrices  1.32.  The  obvious  deviation  from  this  general 
behavior  emphasizes  the  presence  of  a  new  (universality)  class  of  random  matrices 
in  the  transitional  model.  Secondly,  all  the  different  curves  intercept  at  one  point 
in  figure  3.5.  This  striking  feature  could  imply  a  scaling  point  and  therefore  a  new 
statistics  accompanying  the  model. 

Of  vital  interest  are  also  long-range  integrated  statistical  properties.  The  usual 
measures  of  this  behavior  are  the  number  variance  and  the  Ag-statistics.  The  differ- 
ence between  the  two  limiting  statistics  is  substantial,  as  the  Wigner-Dyson  model 
implies  a  slow  logarithmic  increase  in  both  measures,  indicating  a  stiffness  of  the 
spectrum  due  to  the  level  repulsion  compared  to  the  linear  increase  in  the  Poisson 
statistics.  The  transitional  model  recovers  both  extremal  statistics  as  well  as  interpo- 
lates bet-ween  these  two  curves  for  intermediate  values  of  the  transitional  parameter. 

The  result  for  the  Ag-statistics  displayed  in  figure  3.6  supports  that  claim.  The 
analysis  of  these  measures  indicates  that  the  statistics  can  be  written  in  the  form 

A3(yV)  ^aln(27ryV)  +  6iV  (5.1) 

The  Unear  term  in  N  appears  only  in  the  transitional  model  and  will  dominate  the  A3- 
statistics  or  number  variance  for  sufficiently  large  N.  The  slope  of  the  linear  increase, 
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the  constant  b.  solely  and  nonanalytically  dependent  on  the  transitional  parameter  7, 
can  be  used  as  a  measure  of  how  far  the  transition  from  regular  to  chaotic  behavior 
has  progressed  [47,  51]. 

Furthermore  of  interest  is  the  Fourier  transform  of  the  two-level  correlation  func- 
tion, the  two  level  form  factor  6(A;,7).  It  can  exactly  be  determined  from  R2  and 
displays  a  robust  and  measurable  feature  of  a  chaotic  spectrum,  known  as  correlation 
hole.  The  transitional  model  provides  a  precise  description  of  the  correlation  hole  and 
how  it  will  disappear  when  the  system  crosses  over  to  a  regular  spectrum. 
This  characteristic  is  already  observable  in  atomic  as  well  as  molecular  systems  and 
should  be  furthermore  measurable  in  small  mesoscopic  systems  at  small  temperatures. 

The  system  properties  mentioned  so  far  concentrate  on  the  translational  invariant 
part  of  the  model  in  the  bulk  of  the  spectrum.  However,  as  this  model  breaks  spon- 
taneously the  unitary  symmetry  associated  with  the  translational  invariance  [51],  the 
validity  could  be  also  extended  to  include  the  energy  level  edges.  The  theoretical 
predictions  from  this  full  model  cover  the  appearance  of  long-range  order  oscillations 
in  the  correlation  functions,  which  will  appear  as  a  qualitative  change  in  long-range 
statistical  quantities. 

An  important  generalization  of  this  model  is  the  inclusion  of  an  external  parameter 
and  its  interaction  with  physical  systems  exhibiting  transitional  level  statistics.  As 
many  systems  (for  example  disordered  conductors)  are  usually  measured  with  an  ex- 
ternal field  present,  the  importance  of  this  extension  becomes  immediately  obvious. 
This  extension  also  allows  to  study  new  directions,  as  correlation  functions  solely 
depending  on  the  external  parameter  can  be  analyzed.  One  fundamental  result  is 
the  occurrence  of  a  characteristic  quantity  in  the  external  perturbation  describing  the 
sensitivity  of  the  eigenvalues  to  the  external  perturbation.  In  case  of  a  pure  Wigner- 
Dyson  level  statistic  system  a  well  established  relation  exists  using  a  fluctuation- 
dissipation  type  of  equation  which  interprets  the  level  sensitivity  as  conductance, 
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as  originally  pointed  out  by  Thouless  [96,  97].  The  scaling  with  the  characteristic 
quantity  results  in  universal  correlation  functions.  However,  the  generalization  to 
transitional  statistics  will  result  in  a  breakdown  of  this  universality  at  larger  disorder 
(soft  chaotic  spectra).  Furthermore,  the  relation  of  the  level  sensitivity  to  measurable 
quantities  becomes  unclear  but  some  relation  to  the  conductance  should  remain. 

5.2     Relevance  of  the  Model  for  Physical  Systems 

The  questions  approached  by  the  transitional  model  are  relevant  for  a  large  variety 
of  physical  systems.  One  very  important  example  are  disordered  conductors.  This 
topic  has  been  widely  investigated  both  experimentally  and  theoretically,  because 
the  disorder-induced  metal-insulator  transition  is  important  to  our  understanding  of 
disordered  materials.  Furthermore,  the  phenomenon  of  localization  is  fundamental 
in  its  own  rights.  In  particular,  the  connection  between  disorder  and  dimensionality, 
which  has  been  shown  to  be  nontrivial  in  Abrahams  et  al.  [98]  by  exhibiting  a  phase 
transition  only  for  c/  >  2,  has  not  been  understood  on  a  fundamental  level. 

The  general  belief  in  the  community  is  that  the  metal-insulator  transition  at  d=3 
is  sufficiently  described  by  the  3-dimensional  tight-binding  Anderson  Hamiltonian 
[99].  While  the  measurable  physical  properties  like  the  critical  exponents  have  been 
extensively  investigated,  a  clear  picture  about  the  underlying  energy  level  structure 
is  missing  so  far.  However,  there  have  been  several  studies  of  the  eigenlevel  statistics 
of  the  3D  tight  binding  Anderson  Hamiltonian  [48,  100,  101]  which  clearly  suggest 
that  the  transitional  model  describes  the  eigenvalue  spectrum  of  the  metal-insulator 
transition  quite  well.  In  particular,  the  Ag-statistics  in  figure  3.6  exhibits  the  same 
behavior  as  in  reference  [48].  It  should  be  noted  that  not  only  the  large  N  linear  term 
is  present,  but  that  also  for  small  N  the  agreement  is  remarkable. 

Reference  [53]  investigates  the  nearest  neighbour  spacing  distribution,  where  the 
results  show  a  very  similar  'maximum  shifting'  behavior  as  observed  in  the  present 
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model  in  figure  3.5.  Also  visible  from  figure  3.5  is  the  intersection  of  all  curves  in  one 
point,  whicli  is  also  observed  in  other  efforts.  The  appearance  of  a  common  'scaHng' 
point  raises  the  intriguing  possibihty  of  a  universal  behavior  in  the  critical  regime 
and  therefore  an  independent  statistical  distribution  [63,  64].  This  is  also  clearly 
supported  by  numerical  investigations  which  indicate  scale  invariance  at  the  critical 
point.  The  scaling  hypothesis  in  [53]  focusses  on  the  decay  exponent  of  the  spacing 
distribution  for  large  arguments.  A  similar  approach  has  also  been  undertaken  for 
the  present  transitional  model  [82].  At  the  moment  this  is  still  an  undecided  and, 
thus,  strongly  researched  topic. 

In  summary,  the  agreement  both  in  qualitative  as  well  as  in  quantitative  respect 
between  the  energy  level  statistics  in  the  presented  transitional  random  matrix  model 
and  the  3D  tight  binding  Anderson  model  is  quite  good.  However,  there  are  recent 
indications  that  the  spectral  statistics  at  the  transition  point  are  universal  and  that 
they  follow  a  unique  statistical  distribution  [63,  53]. 

Another  area  of  physics,  where  a  transitional  model  is  highly  desirable,  is  quan- 
tum chaos  theory.  Usually,  chaotic  systems  can  only  theoretically  be  approached  in 
semiclassical  terms  based  on  the  Gutzwiller  trace  formula.  The  Gutzwiller  trace  for- 
mula allows  to  calculate  the  density  of  states  as  a  sum  over  closed  orbits.  This  results 
therefore  in  direct  information  about  the  energy  spectrum,  which  makes  a  comparison 
possible.  Furthermore,  many  quantum  chaos  systems  permit  numerical  treatment,  so 
that  statistical  analysis  is  the  preferred  way  of  comparison. 

Even  though  a  direct  link  between  standard  random  matrix  theory  and  complete 
(hard)  quantum  chaos  is  so  far  not  established,  there  is  no  serious  doubt  in  the 
scientific  community  that  the  Wigner-Dyson  statistics  govern  the  energy  levels  of 
complete  chaotic  quantum  systems.  Since  many  systems  also  exhibit  transitional 
spectral  features,  the  interest  in  a  suitable  random  matrix  theory  is  natural.  The 
hydrogen  atom  in  a  strong  magnetic  field  represents  an  important  and  well  studied 
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problem  [102,  103],  and  the  statistical  energy  level  properties  are  very  similar  to  the 
herein  studied  transitional  model,  when  the  strength  of  the  apphed  magnetic  field  is 
increased. 

Both  the  nearest  neighbour  spacing  distribution  and  the  Ag-statistics  have  the  same 
qualitative  behavior  as  the  transitional  model.  A  comparison  with  quantum  billiards, 
however,  gives  a  different  picture.  The  As-statistics  show  for  integrable  systems  a 
cutoff  point,  from  which  on  the  function  remains  a  constant.  This  cutoff  point  can  be 
explained  by  the  presence  of  a  smallest  periodic  orbit  in  the  billiard.  As  long  as  the 
level  number  is  below  the  one  of  the  smallest  periodic  orbit  number,  the  behavior  is 
Poisson-like,  but  when  it  is  crossed,  the  As-statistics  level  off. 

This  picture  is  more  involved  in  the  case  of  chaotic  systems,  even  though  the 
indications  point  into  the  direction  of  a  logarithmic  increase  of  the  Ag-statistics  and 
a  cutoff  at  a  specific  point  [104].  It  is  interesting  to  note  that  this  behavior  is  also 
observed  for  the  transitional  model  in  case  of  the  integration  over  the  entire  kernel. 
The  cutoff  appears,  as  soon  as  the  origin  is  included  in  the  integration  limits. 

The  field  of  mesoscopic  physics  knows  many  applications,  where  transitional  statis- 
tics are  of  importance.  Well  known  examples  are  Aharonov-Bohm  rings  and  qua- 
sicrystals.  One  important  system  are  small  metallic  particles.  The  interest  in  these 
systems  started  with  the  investigations  of  Gor'kov  and  Ehashberg  [32]  who  already 
suggested  that  random  matrix  theory  applies  to  this  type  of  system. 
The  correlation  hole  proposed  in  chapter  3  should  already  be  an  observable  and  ro- 
bust feature.  Measurements  have  so  far  been  limited  to  atomic  or  molecular  systems. 
However,  this  property  should  be  also  observable  in  systems  of  mesoscopic  size.  It 
is  quite  important  to  choose  temperatures  and  particle  size  small  enough  in  order  to 
avoid  continuously  broadened  levels,  which  would  not  exhibit  level  correlations.  For 
a  metallic  system  the  assumption  of  an  electron  gas  would  lead  to  reasonably  sharp 
levels  for  about  50  nanometer-particles  with  the  temperature  in  the  miUikelvin  range. 
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Of  further  interest  would  be  the  choice  of  metallic  particles,  which  upon  doping  with 
another  material  become  insulating.   The  correlation  hole  would  become  smaller,  as 
predicted  by  equation  3.79. 

5.3     Comparison  with  Other  Approaches 

Tho  generalization  of  random  matrix  models  to  include  transitional  statistics  has 
also  been  pursued  with  other  approaches.  It  is  quite  important  for  the  understanding 
of  these  systems  to  explain  similarities  and  differences  between  these  approaches  and 
the  presented  model. 

Banded  random  matrices  have  been  a  well  studied  approach  for  expanding  standard 
random  matrix  theory.  The  leading  idea  in  this  model  is  to  weight  the  matrix  elements 
depending  on  their  distance  from  the  main  diagonal  of  the  random  matrix.  In  case 
this  decay  is  strong,  only  the  eigenvalues  on  the  main  diagonal  will  remain  recovering 
Poissonian  statistics.  In  case  that  the  drop-off  is  extremely  small,  the  matrix  elements 
will  have  approximately  the  same  size  recovering  Wigner's  ansatz.  This  transition 
occurs  as  a  function  of  the  parameter  b/N,  where  b  is  the  bandwidth  and  N  is  the 
number  of  eigenvalues  [35,  61]. 

The  model  can  be  mapped  onto  a  nonlinear  cr-model  and  solved  exactly  for  /V  =  oo 
[59,  60].  The  solution  exhibits  a  smooth  crossover  from  Wigner-Dyson  to  Poisson,  and 
transitional  features  seem  to  be  absent.  The  analogous  behavior  of  the  parameter  ^/L 
in  quasi  one-dimensional  systems,  where  <f  is  the  localization  length  and  L  the  system 
size,  might  suggest  a  suitable  applicability  for  these  systems. 

The  preferential  basis  model  [55],  already  mentioned  in  the  second  chapter,  is  also 
capable  of  describing  the  statistical  crossover  from  Wigner-Dyson  to  Poisson.  It 
introduces  a  unitarity  breaking  term  into  the  well  known  Wigner  probability  density. 

P{H)  =  Cn  exp{-TrH'  -  h'N'Tr{[A,  H][A,  H]^))  (5.2) 

The  matrix  A  determines  a  preferential  direction  in  matrix  space  and  will  be  aver- 
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aged  out.  After  averaging  the  model  will  regain  unitary  symmetry,  even  though  basis 
preference  is  implicitly  contained  in  the  model.  This  procedure  is  quite  reminiscent  of 
the  transitional  model  treated  here.  The  construction  of  the  theory  conserves  unitary 
symmetry  as  seen  from  the  JPD  in  equation  1.18,  but  the  transition  from  a  strongly 
to  a  weakly  confining  potential  suggests  also  here  a  spontaneously  broken  symmetry. 
Other  results  also  underline  the  strong  correspondence  between  these  two  models. 
The  most  striking  similarity  consists  of  the  fact  that  in  the  critical  regime  of  the  pref- 
erential basis  model  the  two  point  correlation  functions,  as  given  by  equation  3.61, 
are  identical.  The  preferential  basis  model  shows  three  different  regimes,  which  are 
separated  by  a  product  of  the  number  of  eigenvalues  times  a  nontrivial  parameter  de- 
termined by  the  coefficient  of  the  original  symmetry  breaking  term.  This  is  analogous 
to  the  consideration  of  the  product  7yV  in  the  weakly  confining  potential  model,  as 
suggested  in  chapter  3. 

In  summary,  the  strong  similarity  with  the  preferential  basis  model  supports  the 
picture  of  the  /V  =  cx3  model  to  characterize  the  critical  regime  and  the  parameter 
7A^  for  describing  the  whole  transition. 

As  in  the  transitional  model,  the  approach  by  Kravtsov  and  Canah  [51]  is  also  based 
on  the  global  maximum  entropy  approach  and  involves  the  solution  of  the  mean  field 
equation  2.19.  It  concentrates  on  the  aspects  of  the  unitary  symmetry  breaking  and 
uses  a  similar  weakly  confining  potential  as  in  3.1  to  derive  the  correlation  functions. 
It  recovers  the  bulk  kernel  3.61  of  the  transitional  model  and  also  the  off-diagonal 
oscillations  3.101.  However,  the  authors  concentrate  on  the  aspect  of  the  spontaneous 
breaking  of  the  unitary  symmetry  and  associate  the  linear  term  in  the  number  variance 
or  As-statistics  with  models  which  explicitly  break  the  symmetry. 

As  order  parameter  of  this  transition  the  coefficient  of  the  linear  term  was  conse- 
quently suggested.  The  similarity  in  the  model  construction  makes  the  model  treated 
herein  and  the  mean  field  model  in  most  aspects  identical.  Therefore,  the  numerical 
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investigation  of  the  mean  field  model  [47]  can  be  used  without  hesitation  to  addition- 
ally support  claims  in  here. 

5.4     Conclusion  and  Outlook 

In  summary,  the  transitional  model  with  weakly  confining  potential  has  provided  us 
with  new  insight  in  the  cross-discipline  area  of  transitional  statistics.  The  two  point 
correlation  function  has  exactly  been  determined,  an  expression  which  has  appeared 
since  from  other  approaches  [55],  and  the  statistical  analysis  contains  many  features 
of  a  transition  in  spectral  statistics  visible  in  many  physical  systems. 

In  particular,  the  analysis  was  able  to  exactly  predict  a  linear  term  in  the  number 
variance  and  the  Aa-statistics.  This  linear  term  has  been  clearly  observed  [48]  and 
presents  a  sign  of  spontaneously  broken  symmetry.  Even  though  it  lacks  so  far  micro- 
scopic foundation,  one  can  seriously  believe  that  the  model  contributes  to  unsolved 
questions  concerning  the  spectral  aspects  of  the  metal-insulator  transition. 

The  scope  of  the  original  model  has  been  extended  to  include  parametric  corre- 
lations and  it  has  been  shown  that  universality  present  in  the  weakly  localized  or 
chaotic  regime  will  break  down  in  the  appropriately  scaled  velocity  correlator,  when 
the  disorder  is  strong  enough  [105].  The  level  sensitivity,  partially  connected  to  the 
conductance,  depends  on  the  height  of  a  peak  smeared  out  by  disorder. 

These  results  demonstrate  two  separate  aspects.  First  of  all,  the  weakly  confining 
transitional  model,  which  would  already  have  been  of  interest  from  the  mathematical 
physics  point  of  view  alone,  shows  unique  and  interesting  phenomena  compared  to 
standard  random  matrices.  Secondly,  the  actual  physical  model  about  energy  level 
distributions  is  quite  successful  in  explaining  real  level  distribution  obtained  either 
from  numerical  or  experimental  work. 

While  the  success  of  this  model  is  encouraging,  further  investigations  are  desirable. 
The  full  two-level  properties  certainly  warrant  additional  studies.  Of  particular  inter- 
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est  are  the  long-range  oscillations  which  are  the  manifestations  of  the  spontaneously 
broken  unitary  symmetry  [51].  The  long-range  properties  show  a  constant  behavior, 
once  the  origin  is  included  in  the  integration  range.  One  can  speculate  that  the  special 
nature  of  the  origin  is  correlated  to  special  points  in  the  spectrum.  The  band  edge 
of  a  disordered  solid  state  system  would  be  one  example,  and  a  statistical  study  of 
numerical  eigenvalues  of  a  typical  system  (tight  binding  Anderson  model)  including 
the  band  edges  could  support  this  point. 

Another  important  case  would  be  the  smallest  periodic  orbit  of  a  closed  quantum 
biUiard  system  (chaotic  or  regular).  These  systems  show  a  clear  cutoff  in  their  A3- 
statistics  at  a  number  N  which  is  associated  with  the  smallest  periodic  orbit.  This 
behavior  can  be  achieved  from  the  current  model  by  using  a  shifted  complete  two 
level  correlation  function. 

An  additional  topic  is  to  expand  the  parametric  correlation  calculations.  The  ex- 
pansion of  the  number  operator  4.44  for  finite  N  to  nonzero  values  of  the  energy,  as 
discussed  in  appendix  B,  would  lead  to  the  straightforward  investigation  of  the  exact 
parametric  number  variance  and  an  improved  treatment  of  the  level  sensitivity.  How- 
ever, their  is  no  clear  way  towards  this  goal,  as  the  large  but  finite  N  nonzero  energy 
asymptotics  of  the  q-Hermite  polynomials  are  presently  not  known. 

Furthermore,  the  connection  between  the  level  sensitivity  and  the  conductance, 
usually  known  as  Thouless  relation,  might  be  extended  to  stronger  disorder  as  sug- 
gested by  the  presented  model. 

A  long  outstanding  topic  is  the  connection  between  this  transitional  statistics  model 
with  microscopic,  i.e.  Hamiltonian,  formulations.  The  encouraging  resemblance  bet- 
ween this  model  and  the  preferential  basis  model  could  eventually  help  to  clarify  this 
point. 


APPENDIX  A 
THE  TWO  LEVEL  FORM  FACTOR 


The  calculation  of  the  two  level  form  factor  in  3.79  is  quite  central  to  the  discussion 
of  the  correlation  hole  and  the  associated  fourier  transform  we  were  not  able  to  find 
in  standard  mathematical  references.  Therefore  a  more  detailed  calculation  has  been 
added  in  this  appendix. 
The  integral  to  calculate  is 

y-oosinh^(6x)  ^^-^^ 

,2/ 


i    f°°  sm  (ax  b)     (.ns 
=     a/  .   \,\^e'^'"^^dx  (A.2) 

0  i-co    smh  (x)  ^       ' 

The  sinus  function  can  be  rewritten  with  the  complex  exponential  function  resulting 


m 


1       rex,    gi[k+2a)ylb  ^      rcx)        ^iky/b  ^      foo    Qi{k-2a)y/b 

=    /1  +  /2+/3  (A.4) 

These  integrands  have  now  a  singularity  at  y=0,  in  contrast  to  the  original  integrand. 

They  have  furthermore  singularities  at  y  -  nrn,  so  that  Cauchy's  theorem  can  be 

applied. 

Studying  the  integral  12  in  more  detail  one  can  see  that  the  standard  half  circle 

contour  is  adequate.    For  k^O  the  contour  can  be  closed  over  the  upper  half  plane 

resulting  in 


Using  the  formula 


^2=^2.zX:/?e._.„(^^^  (A.5) 


Res.^Jiz)  =  y^^j^^i^  -  ariiz)  (A.6) 
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one  gets  [ 

D  ^  ;•  ro/  ■        \       ^  /-,         /  •       X  cosh  2:     ,  .      ^ 

smh  [z]  smh\z)  smh(2)  ^        ' 

+     {z-i7rnY  .  (A.8) 

smh  (z)  ^       ^ 

Expanding  the  hyperbolic  functions  around  the  singularity  2nd  order  gives  then 

smh''(2)  b  ^       ' 

.A:,        1  L 

and  (taking  the  principal  value) 

;2(i  >  0)  =  ^(M  _i_)  (A.„, 

For  kjO  the  singularities  are  at  y  =  -inn  and  the  contour  closes  over  the  lower  half 
plane.  The  above  argument  leads  then  to  the  same  residue  and  the  integral  is 

^^  =  1^4-]nW)  (A-12) 

The  integrals  II  and  13  are  quite  analogous  to  12  with  the  results 

/I  =  -^j2a  +  k\i-—±---  -  1)  (A.13) 


and 


This  results  in 


262'  1^1  _  e-|fc+2a|7r/6         2' 


^^-^l^-^^l^l-e-Uk/^-^)  (^-14) 


^    =     26^l^  +  ^^l(l-e-lUw6-^)  (A.15) 

Trlfcl    1  1 


'     62    '2       1  -e-IM-/*^  ^^-^^^ 

+  2^'^-^^'(l-e-UK/^-^)  (A-17) 

(A.18) 


no 

from  which  we  finally  obtains 


^(^)     =     ^[(|A;i  +  27r)coth[(|fc|  +  27r)^//3]  (A.19) 

+     (|/c|-27r)coth[(|fc|-27T)^//?]  (A.20) 

-     2|A;|coth(|/c|^//?)]  (A. 21) 


This  is  the  expression  for  the  two  level  form  factor. 


APPENDIX  B 
THE  VELOCITY  CORRELATOR 


To  calculate  the  velocity  correlator  from  the  standard  form  of  the  density-density 
correlation  function  a  rather  technical  procedure  was  followed,  which  is  the  reason  it 
appears  in  this  appendix. 
The  correlator  is  defined  as 

C{E,E'.X^X')=  r  r  ^L^S{E,E',X,X')  (B.l) 


dXdX 
where 

N-l       CO 

S{E,E',X,X')    =    cw{yfcE)w{^cE')J2  12^n{smh{y^E))hn{smh{^cE')) 

'l  m=0  n=N 

I  hm{smh{^cE))hm{smh{V^E'))e^^''-^'>^^^''-^'y  (B.2) 

The  constant  c  and  u  are  scaled  in  a  way  so  that  S{E,E\0,0)  =  a{E)  =  1.  From 
equation  ()  they  are  c  =  ■jTr^/2Ng  and  tu  =  c/ii.  The  calculation  of  C{E,E\X,X') 
has  to  be  limited  to  E  —  E  =0,  because  for  finite  E  the  q-Hermite  polynomials 
are  very  difficult  to  estimate.  The  reason  is  that  the  polynomials  h^  exhibit  different 
behaviour  depending  on  whether  the  argument  sinh(v^£^)  is  larger  or  smaller  than 
■yn.  In  case  it  is  smaller  the  curve  is  very  close  to  the  Hermite  polynomials,  compared 
to  a  large  q-deformation  for  larger  arguments.  An  approximation  procedure  capturing 
this  crossover  in  the  polynomials  has  not  been  found  so  far,  but  it  would  be  certainly 
of  interest  in  order  to  expand  the  results  obtained  here. 
For  E  =  E   =  0  the  velocity  correlator  has  the  following  form 

jV  — 1      oo  g 

C{X,X')  =  cJ2J2i         h4smh{V^E))h^{smh{^E))w{^E)dE)  (B.3) 

(  r    /i„(sinh(y^£'))/^„(sinh(^^£:')).^(v^EV£;')— ^e(™'-"'M^-^')^ 

•^-oo  oXuX 
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The  derivatives  in  X  are  fairly  simple  to  calculate,  so  the  problem  reduces  to  evaluation 
of 

/     hn{smh{y^E))hm{smh{y^E))w{y/iE)dE  (B.4) 

J  —CO 

For  q  =  1  these  integrals  can  be  approximated  by  using  the  large  N  expansion  for  the 
normalized  Hermite  polynomials  if„ 


H2n{V2nE)    =    (-1)"^^= -^cos(^2c(4n  +  l)E)e^^  (B.5) 


TT     n 


H2n+i{V2nE)    =     (-ir24= ^sin(y2c(4n  +  3)£;)e'^^'  (B.6) 

Note  that  the  factor  y/c  is  replaced  here  by  \/2c,  which  is  necessary  because  of  relation 
3.24.  This  results  in  the  integral  representation  of  the  density-density  correlation 
function  [88]  (c  ~  ^) 

S{E,  E  ,X-X)  =  jds  j  dt  cos{tt{^E')s)  cos{7T{E-E')t)e^'-''^^^-^  ''^'/^'^  (B.7) 
The  velocity  correlator  is  then 

,  ,  /"l        /■oo  71-2  2        y2 

C{E,  E  ,  X-;^')  =     ds  /  ^tJ,,(2^(i^-5^)  -  4(^)^  — (^^-.2)^)e(^^-*^)^V2.  (g.S) 
Jo     Ji  Zfi  2^X11  ^       ' 

where  the  expression  Js^t  stands  for 

Js,t^  f      f  \os{ir{E- E')s)cos{tt{E- E')t)dEdE'  (B.9) 

The  term  Js,t  can  be  evaluated  by  integration  and  with  some  straightforward  trigono- 
metric manipulations  to  be 

2    cos(7r(E- £')(.- 0)       cos(7r(£:-i^')(.  +  t)), 

^^'^-^^        (]^7p        + (TTJy ^  ^^-^^^ 

As  expected,  this  depends  only  on  the  difference  E  —  E' .  This  expression  is  needed 
to  recover  equation  4.75  in  the  text. 

The  integral  over  s  and  t  is  dominated  by  the  region  s  ■^  t  ^  I.    Hence  one  can 
transform  the  integral  with  the  substitutions  s^  ~  t'^  —  k;  s"^  +  t"^  —  m  approximately 


113 

into  the  expression  (the  integral  over  s"^  +  t^  is  roughly  estimated  with  s"^  -\-  t'^  ^  2) 

a  r  il^.5S±m_lMe^KX-x')^.VUdk  (B.ll) 

By  noting  that  the  extension  of  the  integration  region  to  negative  numbers  (note 
that  always  t  >  s)  will  result  in  the  absolute  value  for  k,  and  that  the  extension  to 
k  — >  ±oo  is  only  appropriate  for  X'^  >>  %  the  integral  results  in 

J —oo  \k^ 

Apart  from  the  constants,  this  is  the  expression  for  q=l  needed  for  the  calculation  of 
the  parametric  number  variance  in  4.44.  Hence  this  argument  can  with  appropriate 
restrictions  on  X,  as  seen  above,  recover  the  hydrodynamic  approximation.  It  is  also 
easily  possible  to  recover  the  well  known  expression  4.29  from  here. 
We  are  mostly  interested  In  E  =  E'  =  0,  because  it  gives  an  X-dependent  quantity, 
without  energy  complications.  Then  equation  B.9  transforms  into 

-    4     s^  +  t^ 

and  therefore 

CiX-X')     =     -it  ds  r  dt^-^e^^'-^'^i^-^'^'-'/'-  (B.14) 

jj,  Jo        Ji         t^  —  s^  ^ 


I  dsj^dt[t'  +  s')e^^'~''^^^-^'^'^'l'^]  (B.15) 


2^ 

The  conjunction  is  now  that  the  same  expression  is  valid  iox  q  <  1,  the  only  difference 
being  the  upper  integration  limit.  However,  this  has  to  be  shown  exphcitly. 
The  integral  4.38  can  exactly  be  evaluated  for  E  =  E'  ={).  First  the  evaluation  will 
concentrate  on  the  case  q=l,  which  leads  to  the  result  4.45  in  the  text.  This  does  not 
only  show  the  procedure,  but  it  also  leads  to  the  relation  which  is  a  result  previously 
not  obtained.  As  before,  we  will  have  to  evaluate 

f    Hn{V2cE)H^{V2cE)e-^''^' dE  (B.16) 

J -oo 


114 


Except  for  the  different  integration  limit  this  expression  represents  the  orthogonality 
relation  for  Hermite  polynomials.  Then  we  see  immediately  that  in  case  of  both  n 
and  m  even  the  integral  vanishes  because  of  the  symmetric  integrand. 


J  —  oo  Z  J  —  ca 


2me-'^"dE 


_nym 


0 


(B.17) 


RecaUing  the  well  known  'raising'  and  'lowering'  operators  for  Hermite  polynomials 
we  have 


dEHn{V2'cE)    =    2^cnH^_,{V2cE 


-cE- 


e=^  dE{Hr.{V2cE)e-'^  }    =    -V2cHn+i{V2cE) 


(B.18) 


Taking  the  case  that  n  is  even  and  m  odd  first,  the  use  of  B.18  and  partial  integration 
gives 

r-o      /-o 


With 


and 


we  get 


/      /     H2n+i{V2cE)H2^{V2'cE)e-'^"  dE 

J  —  oo  J  —  oo 


0        /'O 


(B.19) 
=    2(2n  +  l)/       /      H2n{V2~cE)H2„.-x{V2^E)e-'^"dE  (B.20) 

J—oo  J  —  oo 

=    2\2n  +  l){2m-l)  r    f    H2n-i{^cE)H2m-2{V2cE)e-'^' dE  {^.21) 

J—oo  J  —  oo 


:B.221 


-    2'"^(2n  +  l)(2n-l)---(2n-2m  +  l)(2m-l)(2m-3)---l        (B.23) 
/     HoH2n-2m+i{V2~cE)e-'^' dE  (B.24) 

J—oo 


f        H2n-2m+i{V2^E)e-'''^"dE  =  H2n-2m-2{0) 

J—OO 


^2n(0)  =  (-ir2"(2n-l)!! 


;B.25) 


(B.26) 


H2n{V2'cE)H2miV2'cE)e-'^"dE  = 


2"+"-^(2n  +  l)!!(2m-  1)!! 
2(2n  -2m-  1) 


(B.27) 
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Taking  the  normalization  into  account 

1 


7r22n+2m+l(2n  +  l)!(2m;:    .-oo 

1        [(2n  +  l)!!]2[(2m-l)!!]2 


r    H2n+i{V2~cE) H2m{V2~cE)e-'^" dE]'  [B. 28) 

J  —CO 

(B.29) 


2tt  (2n  +  l)!(2m)!(2n  -  2m  -  1)2 
The  case  m  even  and  n  odd  works  out  in  an  analogous  fashion  (note  that  n  >  m  -\-2) 

.22^.2..H2Ul)!(2m)![/->-(^^)^-^-(^^)--""^-^^l^(^-^Q) 

1  [(2n)!!]^[(2m  +  l)!!]^ 

27r(2n)!(2m  +  l)!(2n-2m-3)2  ^^''^^^ 

This  term  can  be  expressed  with  F  functions  as 

[(2n  4-  l)!f  [(2m  -  l)!!]^  1  r(n  +  3/2)r(m  +  1) 


(2n  +  l)!(2m)!(2n  -  2m  -  1)2    '  tt  r(n  +  l)r(m  +  3/2) 


;b.32) 


anc 


[(2n)!!]2[(2m+l)!!]2  ^  1  r(m  +  3/2)r(n  +  1) 

(2n)!(2m  +  l)!(2n-2m-3)2       tt  r(m  +  l)r(n  +  3/2)  ^       "^^ 

The  velocity  correlator  is  then 

C(X-X')     =     -Lv    f    r(n  +  3/2)r(m  +  l)  1 

'  27r2£'„^^^/^r(n  +  l)r(m  +  3/2)(2n-2m-l)2         ^'^'^^^ 


5^ 


dxdx' 


,(m-n)w(X-X')^ 


(B.35) 


Af/2-l         oo 


1    "i^'     ^      r(m  +  3/2)r(n  +  l)  1 


2^'    i'l  „=^2+i  r(m  +  l)r(n  +  3/2)  (2n  -  2m  -  3)2  ^^'^'^^ 


_^ (7n-n)w(J>!:-X')2 


;b.37) 


Using  an  integral  instead  of  the  sum  and  estimating  ph^?(^!^^">Vl    -   ^^^  (the 

^     1  (m  +  l)l  (n+3/2)  n+n/i      ^ 

approximation  w  ^  leads  to  the  same)  one  obtains 


1         r^^  '^  /"°°  \/m  +  7r/4  1 

^'i^-^)     =     ^-iTT/        d^  dn^/  To ^ ^     (B-38) 

/TT'^yV,  Jm=o        Jn=N/2         L^  j^  7^/4  (2n  -  2m  -  1)2     ^         ' 


g^i■m;-^iT^-|-i;;w^yY-A   T  I'D   OQN 


5^-5;.' 


,((2m)-(2n+l))u;(;C-X  ) 
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1         /-iv/2-i       |-oo  Jn  +  7r/4 


1 


(im  /  (Jn  '^  ^— (B.40) 


2'K'^Nq  Jm  =  \  Jn^N/2  +  l  /^ttT+TtM  (2?^  —  2m  —  3)^ 


<9        ^({2m+l)-2nMX-X-)^  ,^^^y 


\2 


dxdx' 

Using  s^  =  2Nm,t'^  =  2Nn,c  =  ^^uj  =  -  one  gets  (Neglecting  terms  of  order  1/N) 

This  represents  the  integral  approximation  for  the  velocity  correlator  for  7  =  0. 
Now  we  can  switch  to  g  <  1.    Similar  to  the  previous  discussion  the  integrals  with 
both  n  and  m  even  vanish 

/O  1       /*oo 

h2nh2me-'^  dE  =  -  h2nh2me-'^' clE  ^^>^  0  (B.43) 

-co  Z    J —  CO 

The  'raising'  and  'lowering'  operations  are  [91] 

D,K{s\n\i{^cE))     =    2^g-("-i)/'^-^/i„_i(sinh(^£;)) 

K^,{sm\i{^fcE))     =     — -^^D,{K[%mh{^cE))w{^cE)]         (B.44) 

1  '       V  c 

The  partial  integration  rule  can  be  expanded  to  include  the  q-Hermite  polynomials 
[76]  resulting  in 

j     h2n+iismh{^E))h2m{smh{^E))w{^/^E)dE  (B.45) 

=    -     L+i    g'^^'^'^T"^"^    /°    h2m-i{smhiy/~cE))h2n{smhiV~cE))w{V~cE)dE 

Q  i   —  Q'        J— CO 

2^(1 -o)2  .0 

=     Jm+2n+2  (2^  +  1), (2m  -  1 ),  /    /i2m-2 (sinh( v^£;))/i2„-i (sinh( VH^))u;( Vc£:)dis 
y  J —CO 


22'"(1   _  o)2m 

g2.^^^(2^+l)gn(2n+l)(2^  +  l),(2n  -  1),-  ■  •  (2n  -  2m  +  l),(2m  -  l),(2m  -  3),-  •  ■  1 

rO 

ho{smh{^/cE))h2n-2m+l{smh{^/cE))w{^/cE)dE 
Now  with 

rO 

/       ho{smh{^/^E))h2n-2m+i{smh{^E))w{^E)dE  (B.46) 

J  —oo 

l-q 


q2n-2m+\^ 


/i2n-2m  +  l(0)u;(0)  (B.47; 


and 


one  arrives  at 
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1-T2n-2m+l/'i 

h^^iVcE  -^  0)  =  ^^-^(,„_,1,,).;,      ^  (B.48) 


roo 

J     h2n+i{smh{y^E))h2m{smh{^/^E))w{y/^E)dE  (B.49) 

(1  -  g)(2n+2m+.)/2       ^2n  +  l),!!)((2m  -  1),!!) 


q2m+m{2m  +  l)+{2n+l){n+l)  ^2n  —  2m  —  1) 


(B.50) 


here  (2n  +  1),!!  =  ^^~''^^y~.t"  '^'"'  (for  simplicity  (-1),!  =  1)  was  used.  Also  tk 


w..^.^  V-'-   I   -yg--  -  (1-,) 

natural  extension  of  the  factorial  n,!  =  ^^"'^"l^^'C    ^'^  will  be  of  help.  Including  the 

normalization 

„(2n+l)(2n+2)/2+m(2m+l) 
{— ■ -)''" 

liQ]  (i)ooiq;  q)2n+l{q]  q)2m 
and  the  relation 


:B.5r 


(2n  -  2m  -  1),  =  -^{{2n  -  1),  -  (2m),)  (B.52) 


the  integral  becomes 


h2n+l{sm^^E))h2m{sm^^/iE))^v{y^E)dE]^  (B.53) 

u;^(Q)        g^"+^"^((2n  +  l),!!)^((2m  -  1),!!)^ 
7^(9;  g)oo  (2n  +  l),!(2m),!((2n  -  1),  -  (2m),)2 


(B.54) 


The  case  m  odd  and  n  even  {n  >  m  +  2)  is  treated  analogously  with  the  result 

/■° 

[/      h2m+i{smh{^/^E))h2n{smh{^E))wi^E)dE]''  (B.55) 

-'— oo 

^        u;^(0)         g^"+^-((2n-l),!!)^((2m  +  l),!!)^ 

j'{q;  qU  (2n),!(2m  +  l),!((2m  -  3),  -  (2m),)2  ^^"^^^^ 

The  final  expression  for  C(X,X')  is  then 

r(XX')    =    ^\T      f         g^"+^-((2.+  l),!!)^((2m-l),!!)^ 

^     '       ^  27L^o„=4],,  (2n  +  l),!(2m),!((2n  -  1),  -  (2m),)^^^"+^'^-(^-^  ^ 

,      '""f"      V  c/^"^^-((2n-l),!!)^((2m  +  l),!!)^ 

±^1    MWi(2^)-(2m  +  l),!((2n-3),-(2m),)^^^"-^-+^^^-'  ^^ 
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for  A^  >  2.  The  function  f  is  given  by  /„,^(X,X')  =  ^-9!^e((2"^)'-(2"-3)?M^-^')^ 

Analogously  to  the  undeformed  case  it  is  desirable  to  find  an  integral  approximation 
to  this  sum  for  large  N. 
Realizing  that  the  expression  g^,^  =  ^^^S^f'^    ^^^  ^he  two  distinct  limits  a„,„  =~'^° 

^^'^}tnp    ~  ^  °^  ^^^  °^^  ^^^^  and  an,m  ='^^°  ;^  on  the  other  hand,  the  following 
approximation  is  natural 


(B.57; 


TT\/n 


Then  we  can  write  for  the  sum  the  following  approximate  integral  expression  (Neg- 
lecting terms  on  the  order  of  1/N) 


l-K^  jQ  Jn/2       y  m,  (2n,  -  2mg)2         ^ 


1         rN/2  roo  fj^  2n+2m 


Again  the  substitution  n,  =  2Nt,  rUg  -  2Ns  and  use  of  the  relations  c  =  n^j/2Ng ,  u;  = 
iT^/2iJ,Nq  results  then  in 


1 


-.2     I     42 


C{X,X')    =     ![/■  (/5  /-^<ii^^^^e(^'-*')(^-^')-'/2- 
/U  Jo         Jl  s'^  —  P 

7T\x-x'y  /•!     -^ 


j  ds  j -^^  dt[e  +  s'')e^'"-'"^^^-'''>'l^>^]  (B.59) 


2|/ 
This  represents  the  integral  approximation  for  the  velocity  correlator  for  7  >  0. 
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